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Meek’s criterion postulates that a streamer will prop- 
agate if the positive ion space charge field close to the 
anode is equal to K times the impressed field, where K 
should be a constant independent of gap configuration and 
pressure. Measurements in confocal paraboloid gaps, 
designed to substantiate this, revealed that K was not 
constant, but decreased with decreasing pressures. Morton 
showed that on a negatively charged smaller electrode 
Townsend's ionization functions do not represent the 
jonization produced in divergent fields. Consequently it 
was believed that the calculations of the space charge field 
of the positive ions produced in such gaps must be inac- 
curate in view of this erroneous assumption. A series of 
measurements on the dark current i in the same gap showed 


that it was not represented by i=i) exp [/adx] at lower 
pressures. The trend and the apparent quantitative varia- 
tion of K with pressure is explained if the apparent io versus 
voltage curves are extrapolated to streamer onset for each 
pressure. Thus for practical purposes under the conditions 
studied K is sensibly constant in conformity with Meek’s 
theory. Per contra, the applications of Meek’s theory to 
gaps with very divergent fields where the currents can be 
expected to deviate markedly from those computed by the 
relation i=io exp [ f'adx] are not justified, thus placing a 
practical limitation on the use of Meek’s theory to streamer 
studies. Roughly this phenomenon appears when the 
change of the field exceeds two percent change over an 
average electron free path. 


INTRODUCTION 


R some time it has been clear that the 
classical Townsend theory ‘of sparking 
cannot account for the formation of sparks in 
air when the product of pressure and plate 
distance exceeds 200 mm of Hg Xcm. Summaries 
of the inadequacy of the Townsend theory have 
been given by Meek! and by Loeb and Meek.’ 
The discovery of the streamer in the positive 
point-to-plane corona and its explanation by 
Loeb, Kip, and Trichel* supplied a new mecha- 


* Now at the University of New Mexico. 
** Now at the University of California Hospital. 
/ Meek, Phys. Rev. 57, 722 (1940). 
. Loeb’ an J. M. Meek, The Mechanism of the 
Hectic. Spark (Stanford University Press, 1941). 
+A. F. Kip, Phys. ag hy , 139 (1938); G . W. Trichel, 
ibid. 55, 382 (1939); L B. Loeb and A. F. Kip, De App. 
Phys. 10, 142 (1939). 
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nism for spark breakdown. Almostsimultaneously 
the same mechanism was proposed by Raether* 
on the basis of cloud-track pictures. The streamer 
theory of sparking was developed qualitatively 
by Loeb, but it was not until Meek proposed his 
criterion for streamer formation that the theory 
was placed on a quantitative basis. 

Meek’s criterion states that a streamer will 
form when the positive ion space charge field 
close to the anode is equal to K times the im- 
pressed external field (providing the photo- 
ionization in the gas is adequate). The space 
charge field is computed on the assumption that 
most of the positive ions created by the first 
electron avalanches are contained in a hemi- 


*H. Raether, Zeits. f. 
(1938); 112, 464 (1939); H 


107, 91 (1937); 110, 611 
osta, ibid. 113, 531 (1939). 
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spherical boss which has a radius equal to the 
radius of electron diffusion at the head of the 
avalanche. The radius of electron diffusion 
depends upon the quantity Ao/\/f, where Apo is 
the electron free path at 760 mm of Hg, and f is 
the fraction of energy lost by an electron per 
impact. Both Xo and f are functions of the elec- 
tron energy. Since the electron energy distribu- 
tion function is not well known, especially at 
high pressures and in strong electric fields, there 
is some uncertainty in the value of \o/+/f which 
should be used. It has been shown, however, that 


the error introduced by this uncertainty is not 


serious.® 

Originally, Meek chose the empirical condition 
that K=1.0, and his calculations gave a fair 
agreement with published values of sparking 
potentials, but a more careful scrutiny of the 
data indicated that a value of K=0.1 was more 
nearly correct. 

If the streamer theory is to have physical sig- 
nificance, K should be a constant, independent of 
pressure and gap configuration. In order to study 
the behavior of K, one can use values of measured 
sparking potentials in plane parallel gaps. Analy- 
sis of reliable sparking potentials in the literature 
yields a factor of K which varies by a factor of 60 
at atmospheric pressure. More recent measure- 
ments® of sparking potentials yielded data for 
which K varies only by a factor of two. At lower 
pressures, however, there are no adequate spark- 
ing potential values and, as recent studies® have 
indicated, there are not likely to be any reliable 
values of sparking potentials at lower pressures 
for some time. 

Another method can be used to evaluate K. 
Fitzsimmons’ determined the onset potentials 
for streamers in air in confocal paraboloid gaps 
at atmospheric pressure. By knowing the field 
distribution over the axis of the gap at the 
measured potential, it is possible to evaluate the 
positive ion space charge field, and values of K 
may be calculated. The results of Fitzsimmons 
showed that K has the value 0.1 for air, and this 
conclusion seemed consistent with the value of 
K selected by Loeb and Meek from sparking 
potential data. But there was still no information 


5 See reference 2, page 52, last paragraph. 
*L. Fisher, Phys. Rev. 65, 153 (1944). 


7K. E. Fitzsimmons, Phys. Rev. 61, 175 (1942). 


about the behavior of K with pressure. It was 
therefore decided to ,extend Fitzsimmons’ ob. 
servations to lower pressures, and for the sake 
of generality to other gases. The results of these 
measurements, based on the assumptions given 
above, revealed that K is not a constant but 
decreased with decreasing pressure. A carefy| 
analysis of the physical processes involved in 
streamer formation indicated that no known 
factor could be responsible for the results. It was 
therefore concluded that either some unknown 
factor was being neglected, or that Meek’s 
criterion was in principle fundamentally wrong, 
In view of the fact that the criterion, rough as it 
was, appeared to be reasonable, and was based 
on physical processes believed to operate in the 
streamer, it was hoped that a reasonable ex- 
planation could be found. At this time, Morton 
was investigating the validity of applying Town- 
send’s ionization coefficients to ionization pro- 
duced in divergent fields and discovered that 
they were not applicable where the field-change 
per electron mean free path exceeded 2 percent. 
It occurred to one of the writers (L. H. F.) to 
inquire whether the space charge field of the 
positive ions in applying Meek’s criterion to 
divergent gaps was thus not being calculated 
upon erroneous assumptions. Experimental data 
bearing out this hypothesis were obtained in the 
very gap in which the variation of K was ob- 
served. The two sets of experiments will be 
presented separately, while a combined treat- 
ment of both will be given in the discussion. 


THE DEPENDENCE OF K ON PRESSURE 


The experimental arrangements were the same 
as those used by Fitzsimmons’ except for the gap 
itself. A hollow paraboloid electrode of focal 
length four centimeters was used as a plate in 
conjunction with a paraboloid point of focal 
length either 0.02 or 0.005 cm. Both point para- 
boloids were considerably sharper than those 
used by Fitzsimmons. The surfaces of the elec- 
trodes were nickel. These gaps were enclosed in 
a large Pyrex bulb and were connected to a 
mercury diffusion pump and gas _ purification 
system. Mercury, however, was carefully ex- 
cluded from the chamber. Provisions were made 


8 P. L. Morton, Thesis, University of California Library, 
Berkeley, California. 
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to circulate air through the chamber, but this 
procedure was found to be unnecessary if the 
currents passing through the gap were smaller 
than one microampere. If the currents exceeded 
this value by increasing the potential above that 
for streamer onset, oxides of nitrogen were 
formed which changed the streamer onset poten- 
tial. The chamber was baked out under vacuum 
at 400°C for several hours and filled with either 
air or hydrogen. However, after hydrogen had 
been admitted to the chamber, baking alone was 
not adequate to remove the hydrogen. It was 
found necessary to subject the gap to a pro- 
longed glow discharge with nitrogen at low 
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Fic. 1. Streamer onset potentials versus pressure. Curves 
A and B were obtained with a point paraboloid of 0.02-cm 
focal length, curves C and D with a point paraboloid of 
0.005-cm focal length. The arrows indicate the lowest 
pressure where streamers were observed in dry air. 


pressure. Only in this way could the hydrogen 
dissolved in the nickel be eliminated. The air 
was dried by passing it over three traps sur- 
rounded by a mixture of alcohol and dry ice, 
while the hydrogen was cleaned over hot copper 
filings in order to remove oxygen. 

In order to facilitate the pulse formation pre- 
ceding corona onset, a sample of thorianite was 
placed about eight inches from the gap, thus 
increasing the dark current, and making the 
observations more reproducible. 

The chamber was initially filled with the gas 
to approximately atmospheric pressure and the 
potential difference across the gap was increased 
until pre-corona onset streamers appeared.’ The 


pressure was then decreased in steps of 20 mm 
of Hg, and corresponding streamer onset poten- 
tials were observed at each pressure, until a 
pressure was reached at which streamers were 
no longer formed. The pressure at which 
streamers no longer appear is in no case clearly 
defined. On decreasing the pressure of a given 
gas, and maintaining a suitable potential dif- 
ference across the gap, fewer and fewer streamers 
occur, while the number of burst pulses increases. 
This observation indicates that at lower pres- 
sures, the discharge tends to spread laterally 
over the point forming a burst pulse rather than 
to grow out of the high field region in the form 
of a streamer. This preferential formation of 
burst pulses at low pressures is probably due to 
a decrease in the density of photo-ionization in 
the gas. The curves for pre-corona streamer onset 
potentials versus pressure showed a linear rela- 
tionship in all cases (Fig. 1). Within the common 
pressure range for streamer formation, the 
curves for dry air coincided with those for room 
air within one percent. No attempt was made to 
determine the relative humidity of the room air 
used. Hydrogen would not show any pre-corona 
onset streamers if the gas had been purified by 
negative corona from an auxiliary point in the 
chamber® in addition to the hot copper treat- 
ment. However, since all measurements of the 
first Townsend coefficient a in hydrogen were 
taken without electrical purification of the gas, 
this step was omitted in the present experiments. 
Consequently, streamers were observed in the 
potential region preceding corona onset. . 

Since the constant K is the ratio of the field 
X, of the positive ion space charge to the im- 
pressed field X at the tip of the point parabola, 
it was necessary to calculate both X; and X. The 
surface field corresponding to a measured 
streamer onset potential was calculated from the 
relation 


V 1 
X= 
loge (f/F) x+f 


This expression is the solution of Laplace's 
equation for the field on the axis between two 
confocal paraboloids with focal lengths, f for the 
point parabola, and F for the plate parabola, 


(1) 


volt/cm. 


*G. L, Weissler, Phys. Rev. 63, 96 (1943), 
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TABLE I. Values of K and associated factors. wi 
= of 
‘ Pres- Streamer onset N 
sure Point Focal charms 
Gas He v/cm (cm) Xe (xe/p)* v/em Xi/X a ads ait’ pr 
Dry air 748 11950v 1.13X10° 0.02 0.034 7.5X10* 0.67 1210 7.9X10® 9.6108 13.6 
Dry air 654 10600 1.00X10° 0.02 0.034 7.2K10-* 3.210 0.32 1100 4.0x10® 44 129 oc 
Room air 469 8575 8.1010‘ 0.02 0.034 85x10 1.5010 0.18 990 2.4xX10®° 2.4x10® 124 It 
oomair 370 7260 6.85X10* 0.02 0.034 9.6X10-? 3.5310 0.05 880 64x10" 11.2 pl 
H: 766 9000 85X10 0.02 0.09 10.8X10- 39X10! 045 613 1.2X10* 7.35x10* 14.05 sli 
550 7400 7.0X10* 0.02 0.09 12.8x10-* 74X10 0.11 533 3.2X10° 1.7X108 12.7 d 
2 450 6470 6.1110 0.02 0.09 14.110 2.110% 0.04 481 1.1K10° 5.310" 116 
H: 350 5500 0.02 0.09 68X10? 0.013 427 4.6xX10* 2.010" 10,7 ge 
Dry air 746 7070 2.12X10® 0.005 0.017 2.33x10 1.1 3000 7.16108 2.15X10* 13.5 im 
Dry air 547 5810 0.005 0.018 5.74X10-? 0.69 2460 13.2 str 
Dry air 429 5040 1.50X10° 0.005 0.018 048 2250 3.9X10®° 8.8 x108 12.9 
Room air 399 4780 1.43X10° 0.005 0.018 6.7K10-% 3.810 0.27. 2120 2.24108 4.75x10* 12.3 tie 
Room air 252 3700 1.1110® 0.005 0.019 8.7x10-% 7.2108 0.07 1620 7.3X10* 1.18X10* 11,2 str 
| H: 738 5700 1.71105 0.005 0.0475 8x10-% 2.2X10 0.13 1796 1.78108 3.2xK10® 121 
2 478 4370 1.31X10® 0.005 0.044 4.0xX10° 0.03 1400 4.9X10* 6.8610" 10.8 
2 399 3850 1.16X10® 0.005 0.040 10x10-% 8.3310 0.007 1250 1.2xK10* 1.5107 94 
respectively. The quantity x is the distance from hydrogen. Since the integral occurs as an ex- chi 
the point parabola to the point in question. ponent in Eq. (2), its value will be the dominant lat 
For the calculation of the field X; the equation factor. Even if the value of a, is replaced by the ex] 
developed by Loeb and Meek" was employed. value of @ corresponding to a distance of one by 
m3 ionizing free path away from the surface, and - 
| a, exp f adr] even if one allows x, to take on different values, of | 
} ° the final value of X, will be changed only slightly. ust 
X,=5.27X1077 : volt/em, (2) The numerical constant in (2) refers to air fiel 
(x./P) only. In order to calculate the constant for fre 
where the factor 5.27X10-’ contains the nu- hydrogen, it is assumed that the electrons have ra 
merical value of \o/1/f, as is that value of the 4M average energy of about two volts. For this ae 
first Townsend coefficient which contributes the ©4S€ rye gives f= os and Brose and os 
largest number of positive ions to the formation Saayman™* give cm, and thus 
of the space charge boss (and is consequently 0/Vf=1.74X10~ cm. (If we assume an average a 
found at or near the surface of the point), p is ¢lectron energy of seven ev we find do/V/f=48 
| X10-‘ cm.) Thus for hydrogen, Eq. (2) becomes the 
the gas pressure in mmm of Hg, and x, is that cm. or nycrogen, =q. 1 
| distance from the point surface in cm where, on a lati 
: the average, one ion pair is formed by an electron 5.55X 10-7 a, exp f adt_| anal 
advancing toward the point, that is, /o"* adx =1. 0 velt/em. the 
Since the field varies rapidly over the gap, the = j en . 
} calculation of X,; involves the graphical integra- aes 
tion of adx over the axis of the gap. For this The calculated results of the pre-onset streamer nd 
purpose, a must be known as a function of studies at lower pressures are compiled in Table fielc 
distance from the point. Such values are obtained __I. In the case of room air, K varies by a factor of va 
{ from the measurements of a/p as a function of 12 in either gap, while for hydrogen the factor is cyli 
X/p in uniform fields from the measurements of 20 for the more divergent gap and even larger for lh 
Sanders" for air and from those of Hale” for the other gap. In dry air, and in both gaps, XK ‘ale 
apt varies only by a factor of two. Not too much fals; 
1 See reference 2, 45, Eq. (25). ——— _ 
uF, H. Sanders, Phys. Rev. 41, 667 (1932); 44, 1020 4H. Ramien, Zeits. f. Physik 70, 351 (1931). witl 
(1933). “H. L. Brose and Saayman, Ann. d. Physik 5, 797 aes 


2 D. H. Hale, Phys. Rev. 54, 241 (1938); 56, 815 (1939). 


(1930). 
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weight can be attributed to the absolute value 
of K since both x, and a, are not uniquely defined. 
Nevertheless, the decrease in K with decreasing 
pressure appears to be real. 

The range of pressure in which onset streamers 
occur is much larger in room air than in dry air. 
It has been pointed out® that the efficiency of 
photo-ionization is very sensitive to extremely 
slight changes in the composition of a gas. This 
dependency of photo-ionization on purity, to- 
gether with the knowledge that photo-ionization 
in the gas is the primary mechanism involved in 
streamer formation, provides a probable explana- 
tion for the difference in pressure ranges for 
streamer formation in dry and in room air. 


ION CURRENTS IN DIVERGENCE GAPS 


In order to determine K from the onset poten- 
tials of streamers, it has been seen that the space 
charge field of the positive ions must be calcu- 
lated, and that this field is proportional to 
exp [fadx]. It has been recently pointed out 
by Morton® that in the glow discharge, the 
ionization does not take place where the values 
of the field strength are high, and that one cannot 
use the Townsend ionization functions when the 
field strength varies appreciably over a mean 
free path. Such considerations induced Morton 
to study Townsend currents in concentric cylin- 


ders in hydrogen and have led to the most inter- 


esting results. In much less divergent fields than 
were used in the present study, and at pressures 
of the order of a millimeter, he finds that the 
ionization is not given properly by integrating 
the ionization functions for uniform fields. 

The question arose as to whether the calcu- 
lations of the space charge field X, given in the 
earlier part of the paper could be in error due to 
the same phenomenon, thus explaining the vari- 
ations in K. True, the pressures employed in the 
streamer studies were 300 to 700 times as large 
as the pressures used by Morton, but then the 
fields of the confocal paraboloids were very much 
more divergent than the fields of Morton’s 
cylinders. 

In order to decide whether the space charge 
calculations in the confocal paraboloid gap were 
falsified, the more divergent gap was provided 
with an external source of radiation, and the 
currents well below streamer onset were measured 


as a function of potential and pressure. From 
these measurements it is possible to determine 
the importance of the divergent field in calcu- 
lating the space charge field X,. 

The sharpest point parabola (focal length 
0.005 cm) was used and the gap length was again 
four centimeters. Ultraviolet light from a quartz 
mercury arc was focused by means of a quartz 
lens through a graded quartz seal on the center 
of the plate paraboloid. The radius of the circular 
image was estimated as two millimeters. The 
chamber was thoroughly baked out and was 
filled with hydrogen purified over hot copper. 
The current below onset was measured by means 
of a Dolezalek quadrant electrometer (sensitivity 
1000 mm/volt) which was used to determine the 
potential drop across a high ohmage resistor. By 
keeping the ultraviolet arc on continuously 
during the run, the intensity of illumination was 
kept constant. The entire set of readings was 
taken during one continuous run, and the 
current through the arc changed by less than one 
percent. 

Table II gives the different pressures and 


Taste II. Dark currents and related quantities as a func- 
tion of voltage and pressure. 


S adz 
13.6 
12.9 
12.4 
11.2 
14.05 
12.7 
11.6 
10.7 
13.5 
13.2 
12.9 
12.3 
11.2 
12.1 
10.8 
9.4 
nant 
y the 
and 
lues, 
ntly, 
) air 
have Voltage Sads cS eds (ampere) (ampere) 
“a 676 3550 49 
676 3770 67 
$76 3910 
4 
rage 676 4150 41 
4270 
ames 676 4390 
616 3350 
616 3560 
616 3660 
616 3950 
(2’) 616 4060 
616 4150 
548 3110 
$48 3310 
3750 
able 548 3840 
4 
477 3200 
for 477 3480 
s, K 408 2860 
auch 408 2970 
408 3200. 
408 3350 


potentials at which currents were measured. 
Also given are values of fadx and of exp[ fadx ] 
calculated from the field distribution by graph- 
ical integration over the axis of the gap. The 
current versus potential curves are given by Fig. 2 
for various pressures. 

If the current is correctly determined by the 
use of Townsend's coefficient alpha, the measured 
current 7 is given by i exp [ f'adx ], where io is 
the initial photoelectric current. By dividing the 
measured current by exp [fadx], one obtains 
an apparent io which should be constant. The 
departure of 7» from constancy is a direct measure 
of the discrepancy of the true ionization from the 
Townsend value predicted from io exp [ fadx]. 
A plot of the apparent i» against potential is 
given in Fig. 3 for the various pressures. 

Experiments of the nature of those reported 
above were also performed in air, but the currents 
below streamer onset were too small to be 
measured conveniently. The currents in air were 
certainly less than one-hundredth of the currents 
measured in hydrogen. Comparison of values of 
Townsend’s alpha for air and hydrogen leads 
one to expect currents of the same order of mag- 


150 


16 


2500 3000 3500 V = 4000 4500 
Fic. 2. Dark current versus potential curves for various 


pressures. 


nitude in the two gases. The large currents in 
hydrogen must be ascribed to the great trans- 
parency of the gas to electrons and photons, and 
to the small energy losses of electrons in hydrogen 
as shown by the smaller threshold potentials for 
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streamers in hydrogen as compared to air. Kip" 
has reported measurable currents in air with 
positive points below streamer onset only under 
strong radioactive irradiation. 


CONCLUSIONS CONCERNING THE TWO EXPERI. 
MENTS 


It is clear from Fig. 3 that the ionization 
predicted from Townsend’s function is far from 
correct, and that the apparent value of ip is not 
constant even at a constant pressure at varying 
potentials. At the highest pressure the value of 
the apparent % is sensibly constant for a con- 
siderable voltage range and Townsend’s equation 
is not far from correct. This pressure region cor- 
responds to the pressure studied by Fitzsimmons, 
But as the pressure decreases, the deviation Of the 
apparent to from constancy becomes more and 
more marked. At high pressures the mean free 
path of the electrons are short; at low voltages 
all of the apparent io versus voltage curves seem 
to be approaching roughly the same limiting 
value of io. From these considerations it seems 
reasonable to assume that the ionization at high 
pressure and low voltage may be taken as being 
accuractely represented by i=i exp [fadx], 
but that there is more ionization at lower pres- 
sures than is accounted for by Townsend's coef- 
ficient as applied to this equation. Such increased 
ionization at low pressures effectively introduces 
an apparent decrease in the value of K, not 
because of the failure of Meek’s criterion, but 
because the theoretical calculation of the space 
charge field X, with Townsend's alpha does not 
correspond to reality. 

The experiments on the onset of streamers 
actually determined the factor by which the 
actual space charge field in the gap differed from 
the space charge field calculated by means of 
exp [fadx]. Therefore the apparent variation 
of K is not due to a failure of Meek’s criterion. 
The current measurements below onset are not 
only in the right direction to explain the trend in 
K, but the quantiative variation of K with 
pressure is explained in order of magnitude if the 
apparent 1% versus voltage curves are extrapolated 
to streamer onset for each pressure. Thus the 
apparent failure of a constant value of K is to 


4% A. F. Kip, Phys. Rev. 54, 139 (1938). 
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be ascribed to the use of the first Townsend coef- 
ficient in calculations which Morton showed 
cannot be applied when the change of field 
strength per electron free path is approximately 
two percent or Jarger, which applies to these 
investigations. 


AN ATTEMPT TO EVALUATE THE FAILURE OF 
TOWNSEND’S IONIZATION FUNCTIONS 


Between impacts, an electron in a gas in an 
electric field acquires the full energy of the work 
done by the field over the length of the com- 
ponent of the mean free path directed in the 
field direction. Suppose electrons exist in a gas 
in a divergent field such that the field strength 
varies appreciably over the distance of an electron 
free path, and such that the electrons travel 
from strong to weak fields. In such cases, where 
the field changes very rapidly over a few mean 
free paths, as they do about sharp points, the 
electrons gain a large share of the energy of the 
total potential drop (hundreds of electron volts) 
in a very few free paths. Under these conditions, 
the excitation functions and free paths are such 
that energy loss due to excitation and inelastic 
impacts is very small. The electrons gain a ter- 
minal energy nearly equal to the potential drop 
and thereafter can ionize more efficiently than 
they excite. Such ionization is not materially dif- 
ferent from cathode-ray or slow 8-ray ionization, 
where ionization reaches the order of an ion pair 
for 20 to 30 ev of energy expended. 

In a uniform field of the same length as the 
gap, the energy gain of an electron free path near 
sparking is of the order of one electron volt. The 
electrons gain their terminal energy slowly, 
passing through regions in which excitation and 
energy losing impacts are quite frequent before 
and after gaining the ionizing energy. The elec- 
trons rarely gain much more than twice or three 
times the ionizing energy, with the result that 
the average energy expenditure per ion pair is 
much more than the 20 to 30 ev where electrons 
obtain an initial high energy at once and then 
proceed to ionize in relatively weak field regions. 

Thus, under nearly uniform field conditions, 
where Townsend's equation holds, the ion pro- 
duction is small compared to the case where, as 
in glow discharges or about negative points, the 
energy gain is very rapid and is succeeded by a 
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low field region where the energy is expended in 
ionization. 

In the present experiments, however, the point 
is positive, and it would appear as if the ioniza- 
tion should not be greater than for the equi- 


e 

2 an H, 


Fic. 3. The variation of the apparent io with different 
potentials for various pressures. 


librium case where Townsend’s equation applies. 
In fact, it might even be expected to be Jess. In 
the equilibrium Townsend case the value of 
JS adx is constant, irrespective of whether the 
electrons travel from high to low fields, or vice 
versa ; that is, the function fadx is independent 
of the sign of the field. However, where the free 
paths are longer and the field changes appreciably 
over a free path it was seen that as in Morton's 
case, where the high field was negative, the 
ionization could be expected to be, and was, in 
fact, found to be greater than for the uniform 
Townsend case. 

If now we consider the case where the high 
field region is positive with large gradients over 
the free paths near the anode, the situation is 
different. In the cathode region, electrons gain 
energy and may fritter away even more energy in 
excitation and in elastic impacts than in the 
higher Townsend gradients. However, near the 
anode, the sharp gradients over a mean free path 
will again rapidly raise the electron energy to 
possibly hundreds of electron volts over the last 
free paths with very little energy loss to excita- 
tion and elastic impacts. Now if pressures are so 
low that the electrons make very few collisions 
before they strike the anode, then in fact their 
whole energy will be lost and the ionization will 
be Jess than in Townsend’s case. If, on the other 
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hand, the pressure is not too low, the high energy 


of the electrons and their random motions, de- 
spite the high anode field, will cause again a 
greater ionization than in the Townsend case. It 
is not inconceivable that owing to the attraction 
of the anode, these high energy electrons will not 
have as many collisions with gas molecules before 
capture by the anode as they did in moving 
through the low field region in Morton’s case. 
However, both the electrons and their progeny 
generated in the high field region can be most 
effective before capture. This is indicated in 
Bradbury’s'* study of photoelectric currents in 
gases, where back diffusion of electrons is large 
even in the presence of fields, especially when the 
electronic energy is large. The magnitude of the 
excess ionization will obviously be critically 
pressure dependent, and will depend on Ram- 
sauer free paths, excitation, and ionization func- 
tions. The effect will vary in different gases. 
Thus the excess ionization observed in the 
present experiment will be expected, as indicated 


_ above, provided the electrons, before being 


drawn in, spend more time than anticipated in 
the neighborhood of the point. Actually, high 
energy electrons of 100 ev near the anode will do 
just this at higher pressures, as has been recog- 
nized in Bradbury’s study of photoelectric 
currents in gases. Thus it is not surprising that, 
under the conditions of the present experiments, 
the currents were found to be in excess of the 
value of i9 exp [fadx] calculated from Town- 
send’s alpha, though one would expect the 
currents to be relatively less in excess than would 


1%N. E. Bradbury, Phys. Rev. 40, 980 (1932); L. B. 
Loeb, Fundamental Processes of Electrical Discharges in 
Gases (John Wiley & Sons, Inc., New York, 1939), 
Chapter 7. 
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be the case with the high field region at the 
cathode as in Morton’s experiments. 

In the present experiment, the field strength 
changes by 50 percent over a distance of 2.5 
X10-* cm from the point, and the electrons just 
capable of ionization have approximately a mean 
free path of 10-* cm at the lowest pressure 
studied. Hence, in the most extreme case, the 
field strength varies by approximately two per- 
cent over a mean free path. This is the amount of 
variation in the field strength over mean free 
path at which Morton finds abnormalities to 
begin in his currents. 

The effect is undoubtedly a very complicated 
one, and its explanation must lie in the behavior 
of the energy distribution functions of the elec- 
trons as a function of voltage and pressure in a 
divergent field, as well as in the spatial distribu- 
tion of collisions of higher energy electrons with 
atoms. 

The present results in air may be explained by 
the attachment of electrons to oxygen molecules 
in the weak portion of the field. Loeb'’ has shown 
that such ions do not shed their electrons until 
they are in a field of X/p=90 volts/cm/mm. In 
these extremely divergent fields, perhaps, the 
electron is released even closer to the point than 
corresponds to X/p=90. If this explanation is 
correct, the calculation of Townsend currents in 
a divergent gap in air should be approached with 
even more caution than a similar calculation in 


hydrogen. 


The authors wish to express their sincere 


thanks to Professor L. B. Loeb for his many: 


helpful suggestions and for his valuable assistance 
in the interpretation of this work. 


17L. B. Loeb, Phys. Rev. 48, 684 (1935). 
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On the Binding Energy of Deuteron and the Neutron-Proton 
Scattering by a New Potential 


Mu-HsiEN WANG 
Department of Physics, National University of Chekiang, Meitan, Kweichow, China 
(Received June 2, 1944) 


The binding energy of deuteron and scattering cross section of proton by fast neutron are 
calculated by using new forms of nuclear potential suggested by Wang. The results obtained 


are found to be in good agreement with the experimental values, when “zero cut 


potential is employed. 


" of the 


T is pointed out by Wang! that the force 

between two nuclear particles may be related 

to the gravitational force. He takes two alter- 
native forms of the nuclear potential, 


= —Aekir, 
V=—(B/r)eX'', 
K=h/mc=3.84X10- cm, 


(1a) 


and (1b) 


with 


where the constants A and B, as determined by 
the gravitational constant, are 4.78X10-* and 
1.84 10-®5, respectively. 

The purpose of the present work is to deter- 
mine whether the potential (1a) and (1b) can 
give correct results about the binding energy of 
a deuteron and the scattering cross section of the 
neutron by the proton. The calculations follow 
closely those of Bethe and Bacher.? 


I. THE BINDING ENERGY OF 
THE DEUTERON 


The wave equation for the relative motion of 
the two nuclear particles is 


Ay + (M/h?)(E—V)y=0. (2) 


The potential is spherically symmetric; (2) can 
thus be separated in polar coordinates r, 6, ¢, by 
putting 

V(r, 8, &) 


where P;, is a spherical harmonic. The wave 
equation for u; is then 


h? 


1K. C. Wang and H. L. Tsao, Phys. Rev. 66, 155 (1944). 
936)” Bethe and R. E. Bacher, Rev. Mod. Phys. 8, 82 


For the ground state /=0, we have 


M 
V)uo=0. (3) 


dr* 


To solve Eq. (3), two ways of cut-off of the 
potential given by (1a) and (1b) are employed; 
viz., 
(a) zero cut-off: 
V=0, for r<a; 
(b) straight cut-off : 
V=V(a), for r<a; 


(a) Zero Cut-Off 
For r<a, V=0; Eq. (3) takes the form, 
Me 
——-—%™=0, 


where e= —E, the binding energy of the deu- 
teron. Its solution is 


uo = D(e* —e-*r), 
B=(Me/h*)}, 
duo/dr = DB(e*" +e-*"). 

1 duo 
a r>a, V= V(r); the asymptotic solution of 
3) is 


V= V(r), for r>a. 


V= V(r), for r>a. 


where 
and 
Therefore 


coth a. (4) 


Uo 


Now we treat c as a slowly varying quantity 

whose second derivative with respect to r may 

be set equal to zero. In so doing, we have 
duy/dr =c'e—** — Bce~**, 


d*uo/dr? = — 2Bc'e~®" +- 
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TABLE I. Numerical values for 
forr>a and V=0 for r<a. 


«in10-erg coth [1.24 X10"a(Vo 
3.50 1.33 1.38 
3.60 1.32 1.32 
3.70 1.31 1.26 


TABLE II, Numerical values for Eq. (7) when 
V=—(B/r)e*'" for r>a and V=0 for r<a. 


. (7) when 


«in 10 erg coth [1.24 2e 
3.50 1.29 1.38 
3.60 1.29- 1,32 
3.70 1.27 1.26 
3.80 1.26* 1.19 


Substituting this last expression in (3), we get 
the differential equation for c, 


dc/dr=c' = —MVc/2gh’, 


the solution of which is 


Thus 
( MV 
uo= F exp Br r 
uo a 
5 


In order that the wave function u» could be 
joined at r=a smoothly, the expressions (4) and 
(5)-must be equal; i.e., 


(6) 
2Bh? 
(b) Straight Cut-Off 


For r<a, V= V(a); Eq. (3) takes the form 


V(a)| 


sin yr, 


“ v(a)| 


Thus 


1 duo 
( =¥ cot ya. 

Uo dr r=a 

For r>a, (1/uo)(duo/dr) is same as given by 
(5). Therefore, instead of (6), we have 


M|V(a)| 


— (6a) 


Equations (6) and (6a) serve to determine the 
binding energy of the deuteron ¢ if we use the 
experimental value of — V(a)= Vo, the depth of 
the potential well, and calculate a, the range of 
the nuclear force, from (1a) or (1b). Conversely, 
if we use the experimental value of ¢, we can 
calculate the value of Vo, and a by the aid of 


Eq. (1a) or (1b). 


(c) Numerical Calculations 


(1) V=0, for r<a; V=—Ae*’’, for r>a. We 
take the recent experimental value V»=10.5 
Mev; then a calculated from (1a) is equal to 
4.21X10-" cm. Equation (6) becomes numeri- 
cally, 


coth [1.24 (1) 


From Table I we see that the root of (7), the 
binding energy of deuteron, is 3.60X10-* erg 
or 2.26 Mev. 

(2) V=0, for r<a; V=—(B/r)e*’’, for r>a. 
Also take Vo=10.5 Mev; then a is equal to 4.42 
X10-" cm from (1b). 


TABLE III. Numerical values for Eq. (7) when V= V(a) 
for r<aand V=—Ae*’' for r>a. 


Vo0o/2Ve—vVe 
«,in10-erg cot [1.24 X10"a(Vo—e)#] 
7.80 0.018 0.067 
8.00 0.038 0.037 
8.20 0.061 0.027 


TasB_e IV. Numerical values for Eq. (7) when V=Vi@) 
for r<a and V=—(B/r)e*'" for r>a. 


Vo0o/2Ve—ve 
10-erg cot [1.24 
8.20 —0.020 0.206 
8.30 —0.013 0.005 
8.40 — 0.004 — 0.083 
— — 
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BINDING ENERGY OF DEUTERON 


The binding energy of deuteron from Table II 
is 3.66 X 10~* erg or 2.30 Mev. The results given 
by (1) and (2) are in good agreement with the 
experimental value 2.17 Mev. 

(3) V=V(@), for r<a; V=—AeX!", for r>a. 
Also take Vo= 10.5 Mev, then a is 4.21 X10-" cm 
from (1a). Equation (6a) turns out numerically 
to be 
cot [1.24 Vo—e)*] 

=([(Vo/2V/ 6) — Ve]/(Vo—«)!. 

(4) V=V(a), for r<a; V=—(B/r)e*’’, for 
r>d. 

From Tables III and IV, the binding energy of 
the deuteron is 8.0 10~* and 8.3X10~ erg, re- 
_ spectively. These values are too large in com- 


ison with the experimental value. The correct 
value will be obtained if we take Vo=5.93 Mev. 


II. THE SCATTERING OF NEUTRON 
BY PROTON 


Let us denote E the kinetic energy of a proton 
and a neutron in a coordinate system in which 
the center of gravity of the two particles is at 
rest, which is equal to one-half of the kinetic 
energy of the incident neutron in a system at 
rest. The wave function u; will satisfy the 
equation, 

M\ dr 
Asymptotically for large r, the solution of (8) 

u,=c sin (Kr—4$lx+5)), (9) 
with 
K=(ME/h*)}. 

Then the cross section do, the number of the 
neutrons scattered per unit time through an 
angle 6 and 6+48, if there is one neutron crossing 
unit area per unit time in the incident beam, is 
given by the well-known formula, 


do = (x/2K") 
X | [exp (275,) —1]|* sin (10) 
It has been shown’ that if 1/K>>a, a being the 


range of the force, all phases 6; will be small 
except 59. Then 


do sin? 59 sin 6d@, 
*Cf. reference 2, p. 115. 


TABLE V. Comparison of values of « X 10". 


@(obs.) 


1.1-1.5) 


43 
86 


(11) 


For the ground state of deuteron, we have 
already shown that 


Now, in the present case E is positive ; we should 


have® 
1 dugt 
=a-— 
naM(E+e) 


= 
h? 
where uot and uo~ are the up» functions for r<a 
corresponding to the positive and negative values 
of E, respectively ; i.e., 


c= f do=4nK~ sin? do. 


M(E+e) 


(13) 


Uo =B sin Kr, uo = D(e®* ; 


therefore 


1 a 


8 coth Ba—K cot Ka 
B+ K? 


(14) 
For large r, 
1 dugt 


uot dr 

In order to join the wave function at r=a 

smoothly, the expressions given by (13) and (15) 
must be equal, thus 

K cot (Ka+6o) =A. (16) 


Under the assumption that Ka is small, we 
have 


) =K cot (Ka+6o). —_ (15) 


cot (17) 

and by (11), 
o=4n/K*+A?*. 

If » is small A ~a, and we have 


oo=4r/K?+ 02. 


(18) 


(19) 
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-The numerical results are given in Table V 
in taking V=—Ae*!', and the values of a@ are 
obtained from Table I at the point r=a. Cor- 
responding to (18) and (19), Bethe and Bacher 
have derived’ the formulae: oo(B) =4rh?/e+E, 
o(B)=30/2, and the values calculated from 
these are given for comparison in the fourth and 
fifth columns of Table V. 


COx 


It is seen that the present results are in better 
agreement with the experimental results than 
those given by Bethe and Bacher. If we take 
V=—(B/r)eX'", the results do not differ appre. 
ciably from that given above. 

In conclusion, the author wishes to express hig 
thanks to Dr. K. C. Wang for suggesting the 
calculation and for helpful discussions. 
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Momentum and Energy of Photon and Electron in the Cerenkov Radiation 


R. T. 


Cox 


Johns Hopkins University, Baltimore, Maryland 
(Received June 22, 1944) 


The conservation of momentum and energy between an emitted photon and an electron in 
the Cerenkov radiation gives an equation for the direction of emission. This equation differs 
from the one which Frank and Tamm derived from electromagnetic theory only by a negligible 
term involving the ratio of the wave-length of the electron to that of the photon. 


it has been observed by Cerenkov! that fast 
electrons traversing a transparent medium 
emit a radiation, continuous through the visible 
spectrum and beyond, which differs in respect to 
its polarization and angular distribution from 
radiation previously found with similar sources. 
Frank and Tamm? have given an electrodynamic 
theory of this radiation similar to the dynamic 
theory of the bow wave of a ship or the conical 
wave around the path of a bullet moving through 
the air at a speed greater than the speed of 
sound. According to this theory the emitted 
radiation has a component of every wave-length 
of which the speed in the medium is less than the 
speed of the electron. The energy per unit 
length of path radiated by the electron in the 
frequency range dy is (2me/c)?(1—1/6°n*)vdv, 
where e¢ is the electronic charge, c the speed of 
light in vacuum, 8 the ratio to this speed of the 
speed of the electron, and n the index of refrac- 
tion of the medium for light of frequency v. The 
emitted rays make an angle @ with the electron 
velocity given by cos@=1/n8, and they are 

1P. A. Cerenkov, Comptes rendus Acad. Sci. U.S.S.R. 8, 
451 (1934) ; 12, 413 (1936); 14, 102 (1937); 14, 105 (1937) ; 
Phys. Rev. 52, 378 (1937). 

21. Frank and Ig. Tamm, Comptes rendus Acad. Sci. 
U.S.S.R. 14, 109 (1937). The subject has also been treated 
by Fermi, Phys. Rev. 57, 485 (1940), in the development of 


a classical theory of the retardation of charged particles in 
gaseous and condensed media. 


polarized with the electric vector in the plane of 
this angle. The theory accords with the ob- 
servations of Cerenkov, and it has been con- 
firmed by observations under more favorable 
experimental conditions made by Collins and 
Reiling* and, over a considerable range of 
electron speeds, by Wyckoff and Henderson.‘ 
Although the phenomenon has thus an ade- 
quate classical explanation, it is interesting also 
to treat it by considering the conservation of 
momentum and energy between the electron and 
an emitted photon of the radiation. Let it be 
supposed that the electron, traversing the 
medium with a speed u, emits a photon of energy 
hv in a direction making an angle 6 with the 
initial velocity of the electron. After the emission, 
let the speed of the electron be v and let the angle 
between its final and initial velocities be ¢. By 
the conservation of momentum, the angles @ and 
¢ will be coplanar and on opposite sides of the 
initial direction of the velocity of the electron. 
Also, by the conservation of momentum, 


cos 6+ (h/X) cos 6 
=mu(1—u?/c*)-4, 
mv(1—v*/c?)-! sin ¢—(h/A) sin 6=0. 


a pb ha B. Collins and V. G. Reiling, Phys. Rev. 54, 499 
+H. O. Wyckoff and J. E. Henderson, Phys. Rev. 64, 1 
(1943). i 
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MOMENTUM AND ENERGY 


And by the conservation of energy 
me?(1 hy = mc*(1 — u?/c*)-4, 


where m is the rest-mass of the electron and A 
the wave-length of the photon in the medium. 


—That the wave-length in the medium, rather 


than in vacuum, is the divisor of hk in the 
momentum is required by the fact that in quan- 
tum mechanics it is the momentum which deter- 
mines any interference pattern in the medium 
which might be used for the measurement of the 
wave-length} It has also been pointed out else- 
where” that this value for the momentum is 
needed for consistency with the laws of refraction 
and radiation pressure. 

In the paper just cited, Michels and Patterson, 
considering relativity in a non-dispersive medium, 
have also proposed that c/n, the speed of light 
in the medium, should replace c in the expressions 
for the momentum and energy of a particle. If 
there were a non-dispersive medium, the cal- 
culation of Frank and Tamm indicates that in 


‘it a charged particle would have as its limiting 


speed c/n, rather than ¢ as in empty space. For 
the expression (2re/c)?(1—1/6n?)» would have 
to be integrated, with m constant, over all values 
of the frequency » to give the energy radiated per 
unit length of path. This integral being infinite, 
an infinite force would be required to maintain 
the particle at any speed above c/n. Thus the 
calculation of Frank and Tamm provides some 
support, at least of a formal kind, to the con- 
clusion of Michels and Patterson. At the same 
time, however, it shows that the assumption of 
a non-dispersive medium is inappropriate for the 
treatment of the Cerenkov radiation, since in 
this phenomenon the electrons actually traverse 
appreciable distances at speeds greater than the 
speed of the radiation in the medium. 

In the present treatment the use of c/n rather 
than c in the relativistic expressions for the 
momentum and energy of the electron would 
make them and hence also the momentum and 
energy of the photon imaginary at electron 
speeds at which radiation is actually emitted. 


5 L. Brillouin, rendus 178, 1696 (1924). 
*R. T. Cox and J. C. Hubbard, Proc. Nat. Acad. Sci. 
Ml, 1925). 
. Michels and A. L. Patterson, Phys. Rev. 60, 589 


(1944). 
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Justification for considering the modification 
which the medium makes, through its modifica- 
tion of the wave-length, upon the momentum of 
the photon, while ascribing to it no modification 
of the momentum of the electron, may probably 
be based on a comparison of the wave-lengths of 
photon and electron with the interatomic dis- 
tances in the medium. The wave-lengths of light 
for which the index of refraction is appreciably 
greater than unity are so much greater than the 
interatomic distances that it appears reasonable 
to treat the medium as continuous in considering 
the motion of the photons. On the other hand, 
the wave-lengths of fast electrons are so much 
smaller than the interatomic distances that in 
considering their motion it would seem necessary 
to regard the medium from the Lorentzian point 
of view as simply space with a distribution of 
charged particles. Thus the ordinary relativistic 
expression for the momentum of a particle in 
space would be required. 

The elimination of ¢ and v among the three 
equations above, together with the substitution 
v=c/nx, gives the result 


cos 
This may be written 
cos 6 = 1/nB+(A/d)(n?—1)/2n?, 


where 8 = u/c and A is the de Broglie wave-length 
of the electron before the emission of the photon. 
“The first term on the right is the expression for 
cos 6 derived classically by Frank and Tamm. 
The second, in which 4 is a factor, illustrates the 
rule that a quantum equation becomes identical 
with the corresponding classical one when h is 
replaced by zero. Or, to say much the same thing 
in another way, this term is a natural one to find 
as the difference between a quantum calculation 
and a classical one, since the classical calculation 
of course ignores any effect associated with the 
wave-length of the electron. | 
This term will be very small since, as was 
mentioned earlier, the wave-length of the electron 
is very much smaller than that of the photon. 
Hence the classical approximation is a very good 
one, and the direction of emission of the light 
does not appreciably depend upon its wave- 
length except as the wave-length determines the 
refractive index. | 
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Free Vibrations of Anisotropic Bodies 


H. 
Research Laboratory, Commercial Equipment Company, pe City, Missouri 
(Received April 17, 1944) 


Approximate solutions for free vibrations of a finite anisotropic body are derived by a per- 
turbation method. As an example, some extensional modes of thin crystal plates are calculated. 
Calculated frequencies and deformation patterns are compared with observations. 


1. SKETCH OF THE METHOD 


NLY very few problems of elastic vibration 
are susceptible of rigorous solution by the 
methods of the general theory of elasticity. 
During the last years stationary vibrations of 
crystalline bodies have become of practical im- 
portance because of the discovery of piezo- 
electric excitation of vibrations. Rigorous solu- 
tions for free vibrations of a crystalline body are 
known only for an infinitely extended plane 
parallel plate and an infinitely thin rod. It has 
become desirable to obtain even rough approxi- 
mations, as rigorous solutions for most problems 
are inaccessible. Some attempts in this direction 
have been made by Mason! and Bechmann.? 

It is attempted here to give a more systematic 
method of approach. The perturbation method 
has been used for a long time especially in quan- 
tum mechanics. To our knowledge, adaptation 
of this powerful method to the theory of elasticity 
has not been made. 

We first make a resume of well-known equa- 
tions of elastic vibrations. If U, V, W are the 
Cartesian components of the displacement and 
if we have stationary vibrations of angular fre- 
quency they take the form® 


,/dz, 
(1) 
,/dx+0Y,/dy+0Z,/dz. 

The strain components are defined by 
x,=0U/dx, x,=dU/dy+0V/dx, 


Y=9V/dy, (2) 


2,=0W/dz, x,=dU/dz+0W/odx. 


1W. P. Mason, Bell Sys. Tech. J. 13, 405 (1934). 

?R. Bechmann, Zeits. f. Physik 117, 180 (1941); 118, 
515 (1941); 120, 107 (1942). 

3W. Voigt, Lehrbuch der Kristallphysik (B. G. Teubner, 
Leipzig, 1910). 


The general linear relation between stress com- 
ponents and strain components is 


3) 


with 
Cik=Cki; 


and the inversion of (3) is 
—Xy= 2+See Y,+ +566X y. 

If we eliminate the strain and stress components 


from these equations, we obtain three equations 


of the form 
V, W)=0, 


p2’?V+L(U, V, W)=0, 
p2’W+M(U, V, W)=0, 


(4) 


with 


K(U, V, (5) 


Ox 


where K , L, M are linear operators of second 


order which act on three components U, V, W. 
The boundary conditions for the free body 
take the form: 
X,cos (n,x) +X, cos (nm, y) +X, cos (n, z) =0, 
X, cos (n, x) + Y, cos (n, y) + Y, cos (n, z) =0, (6) 
X, cos (n, x)+ Y, cos (n, y)+Z, cos (n, z) =0, 
where n signifies the normal to the surface, 
directed inwardly to the body. 
The fact is remembered that any two rigorous 
solutions of (5) and (6) are orthogonal to each 
other :4 


f (7) 


where the integration is extended over the body. 


4A. E. H. Love, The Mathematical Theory of Elasticity 
(Cambridge University Press, 1934). 
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We shall try to represent the solution (U, V, W) 
approximately as a sum of zero-order functions 
(u;, Vi, Wi) Which are chosen so that they are 
rough approximations of the actual solution. We 
write: 

W=CLAwi. (8) 


However, this expansion is possible only if the 
functions (u;, v;,w,) satisfy the boundary con- 
ditions (6). In the following, we shall not suppose 
this to be true, so that (8) does not hold in 
general. We will construct a set of auxiliary 
functions («;’, v;’, w;’) so that they are equal to 
the functions (u;,v;,w,;) everywhere but in a 
small domain close to the boundary. In ‘this 
domain, (u;’, v;’, w;’) will behave so as to satisfy 
the conditions (6). We will write: 

U= LA V= LA vi, W= XA (9) 

We now define strain and stress functions: 

4 = + (10) 
an 

— =CuXzit (11) 


The functions (u;,v;, w,;) will not satisfy Eq. 
(5), but will satisfy equations of the form: 


—K (ui, vi, Wi) 
+0X,;/d2= — fi, 
—L(uj, vi, wi) 
+0Y,,;/dz= —gi, 
—M (ui, vi, W,) =0X Y,;/dy 
+0Z,;/d2= —hy. 
Instead of the boundary conditions (6), the 


functions (u,;, v;, w;) will satisfy equations of the 
form: 


(12) 


X,; cos (n, x) +--+ =¢,(s), 
Xi cos (n, x)+-+- =yi(s), (13) 
X cos (m, x)+ +++ =x;(s). 

Similarly, the corresponding functions derived 


from the auxiliary functions (u,’, v;’, w;’) will be 
denoted by 


fi, g., and vi, xi’. 


It will be noticed that u; --- can be considered 
as solution of the following problem: Given a 
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volume force of components —f,/p, —g;/p, and - 
—h;/p, respectively, and surface tractions of 
components ¢i, Wi, xi, respectively, the equi- 
librium deformations of the body are to be 
determined. From Eqs. (12) and (13) it can be 
seen that the requirements for the solution of 
this problem are met by %,, v;, and w,. 

Let S’ be a surface internal to the surface S of 
the body and parallel to S, having a small 
distance « from S so that to each point s on S 
there corresponds a point s’ on S’. The cor- 
responding points are joined by the normal n 
which is directed toward the interior. Let the 
domain included between S and S’ be called D, 
while the rest of the body volume is V. The 
tractions across S’ which have the form of the 
first members of (13) will be found by integrating 
the functions —f;, —gi, —h; along the normal n: 


Ox oy dz 
= | X.; cos x) cos (n, y) 


cos (m,2)|nm0, (14) 


OX 
S( dn= cos | nmo- 
Ox 


If we use (12) and (13) we obtain 
{Xz cos (n, x) +Xyi cos (n, y) 


where the dots stand for the two analogous 
equations. 

The same calculation for the functions u,’- -- 
yields similarly : 


{ Xi; cos (m, x) +Xy; cos (n, y) 
+X% 08 (m, f (16) 


Returning to the definition of the auxiliary 
functions u,’---, we remember that they should 
be equal to u;--- within V. Therefore, they 
must be solutions of a problem involving the 


1944 
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| 
7 
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same forces —f;/p--- within V and the tractions 
across S’ given by (15). Thus, the first members 
of (15) and (16) are equal: 


By definition of the u,’---, the tractions 
must vanish: 


‘fidn= f 


f 


We substitute (9) into (5): 
pO? DA + LAK (u,’, w,’)=0, 
(19) 
pQ? w’) =0. 


We multiply these equations by u,’, 1%’, wx’, 
respectively, add them, and integrate over the 
crystal : 


{ pQ? + + w;'w,')dv 


M(ul, (20) 


As u,’--- tends toward u;---, the first integral 
in (20) becomes: 


f Nx. (21) 


In view of (7) the non-orthogonality integrals 
N« would vanish for 1+ if the functions (u;- - -) 
were exact solutions of the problem. 

The integrals N;; will be determined arbitrarily 
by normalization : 


1= Nx. (22) 


The volume V inside of S’ contributes to the 


second integral of (20) the term: 
The contribution of D to this integral is 


D D 


where u,’ is put equal to u% and considered as a 
constant along the normal n. Substituting (18) 
into (23) we find 


+0 J | | 


When e tends toward zero, the integrals 
Jo'‘fidn--+ vanish and we have 


=— 
Equation (20) now takes the form: 


(uj, ---)}dv 


- f 0, 
Ai(o2?Na — Hx) =0, (24) 
where 


+ f (25) 


If we consider, as above, two sets of functions 
Ui, Vi, W; and Uy, M%, We as solutions of two dif- 
ferent equilibrium problems involving the given 
volume forces and surface tractions, 


—fi/p, —gi/p, —hi/p; gis Wir Xi 


d 
—fi/p, — gr/P, —hy/ p; Pky Vey Xk 


or, 


= 
{ 
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respectively, then a theorem due to Betti* states 
that: 


- f (fiat 


+ f 
=— f (fats t+ gedit 


+ f (26) 


or, according to Eq. (25): 


Ax = (27) 


In other words, the matrix Hy, is symmetric. 
If the Eqs. (24) are to be compatible, the 
determinant 
p?— Hy 
p2?Ni2— 


pQ?Ny2— Hy 


(28) 


must vanish. This will be possible only when 2? 
assumes certain characteristic values which cor- 
respond to the approximate characteristic fre- 
quencies of the free vibrations. With a particular 
value 2,2 we can solve Eqs. (24) for the coef- 
ficients A; and by substitution into (9) we can 
obtain the approximate solutions U, V, W. For 
all practical purposes, the functions u;--- can 
then be substituted for the u,’; in other words, 
Eq. (8) can be used instead of (9). . 


2. APPLICATION: EXTENSIONAL VIBRATIONS OF 
THIN CRYSTAL PLATES 


Let a crystal plate have its two parallel larger 
surfaces normal to the Z axis of the reference 
system. The origin is situated in a plane midway 
between the layer surfaces. On the boundaries, 
for a free crystal 


X,=Y,=Z,=0. (29) 


For a thin crystal, we may assume that the 
variation of the stress is linear, in a first approxi- 
mation. As (29) is true for both surfaces, it must 
then hold everywhere. Substituting (29) into (4), 
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we obtain 


—Vy= y, 

Vy + SeeX y. 

And if we solve (30) for the stresses: 

—-X,= t Very, 

Very, 

—X y= Yoery- 

In this equation, X,, Y,, X, mean the average 

value of the stresses taken over the thickness of 

the plate. Equation (31) is the equation given by 

W. Voigt.* Explicit formulas for the y% are given 
in the same book. 

As explained by Voigt, Eq. (31) applies really 

only to points of the plane z=0. For these 

points, W=0 in the case of extensional vibra- 


tions, and the strains (more exactly, the average 
values of the strains) are given similarly to Eq. 


(30) 


(31) 


This problem has been solved for the special 
case of a circular isotropic plate by Love.‘ For a 
rectangular plate, no general solution is known. 
V. Petrzilka® has given an exact solution for a 
mode of an isotropic rectangular plate, but the 
solution exists only for certain dimensions of the 
rectangle. 

For this two-dimensional case, Eqs. (10)-(13), 
(21), (22), and (25) are simplified : 

Ou; Ov; 


OX OX yi OX yi 

cos (n, x) cos (n, y) = ¢:(s), 

Xi cos (n, x) + Y,; cos (nm, y) 

*V. Petrzilka, Zeits. f. Physik 97, 436 (1935). 


Ou; (10a) 


(11a) 


(12a) 


(13a) 
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f (uj4p+0,0,)dxdy = Nix, (21a) 
f =1, (22a) 


— f 
+ (om+van)as. (25a) 


(a) Low Frequency Longitudinal Modes 


A group of extensional modes has been called 
longitudinal because the strain components x, 
and y, are large in comparison to the shear xy. 

These vibrations have been studied experimen- 
tally and theoretically by several authors.'?%7 
Their analysis is, however, unsatisfactory in 
several respects, mainly because it fails to 
account for the fact that three low frequency 
longitudinal modes are observed, usually with 
nearby frequencies. | 

We shall obtain the zero-order solutions from 
consideration of a simplified problem. In (31), 
we put 


Yu = Y22; 
V12 = = = 0. (34) 
Yes = 


These are the equations for an isotropic plate 
with vanishing Poisson’s ratio. Furthermore, we 
suppose that the plate is square, the edge having 
the length a. The directions X and Y shall be 
parallel to the edges, and the origin is in one 
corner of the plate. One exact solution is the 
one given by Petrzilka: 


u, = C; cos kx sin ky, 
where k=x/a. (35) 
v= —C, sin kx cos ky, 


The constant C, is arbitrary, but will later be 
defined by the normalization (22a). 

Two other solutions are the longitudinal plane 
waves in the directions X and Y, respectively: 


U2= C2 cos kx, v2=0, (36) 
and 
us=0, v3=C3 cos ky. (37) 


* A. Lissiitin, Zeits. f. Physik 59, 265 (1930). 
7V. Petrzilka, reference 5, and Hoch: tech. u. Elek: akus. 
50, 1 (1937). 


These three modes are degenerate, i.e., they 
have the same frequency : 


pw? =k*y11. (38) 


We now abandon the simplifying assumptions 
(34). We then will find the familiar phenomenon 
of splitting of a degenerate eigenvalue by a per- 
turbation. 

For the more general case, we let two unequal 
edges of lengths a and 6 be parallel to X and Y, 
respectively. The zero-order functions can be 
written: 

C; cos sin key, 


= —C, sin cos key, (39) 
u2=C2 cos v2=0, (40) 
us=0, v3=C;3 cos key, (41) 


h 


The normalization (22a) requires that: 
Ci? = C2? = C3= 2/ab. (42) 
The non-orthogonality integrals Nx are, ac- 
cording to (21a): 
Nis=—2/r, Nos=0. (43) 


The equations (10a), (11a), (12a), (13a), and 
(25a) yield after some elementary calculations: 


y= — (Re — Ri) 
Hy: = Hy = (2/4) (hi? yu — 


Hy yu ( 44) 
Hy = Hy = (2/ J, 
= 


Fos = = (8/*) 


These are the elements of the secular deter- 
minant (28). 

We first discuss the case where the major face 
is normal to a crystallographic threefold axis. In 
this case, according to Voigt, the plate behaves 
as if it were isotropic, i.e., 


(45) 
If the plate is also square (ki=k2=k) the ex- 
pressions (44) reduce to: 
1), 
kyu, Hy = — 112), (46) 
= k*y11, = k?(8/?)y12. 


| 
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The secular determinant (28) has the form: 


One root of (47) is obvious from inspection: 
= k*(y11— (48) 


With this value, the ratio of the coefficients A; 
in (24) is: 

A,™:Ae™:A3@ =1:0:0. (49) 
The first solution is, according to (8) and (39): 


Ui=C; cos kx sin ky, 
(39) 
Vi=—C, sin kx cos ky; 
ie., the result is the unchanged function (39). 
This was to be expected, for (39) is known to be 
a rigorous solution for the isotropic case with 
a=b. 
To find the two other roots of (48), only a 
quadratic equation has to be solved. The result 


is: 
(50) 
and with these values 


A, =0:1:1. (51) 


TaBLe I. Frequencies of square Z-cut quartz plates with 
edge a=1 mm (in kilocycles). 


Mode Neale. Nobs. 
1 2.537 108 
2.564 108 2.551 108 
2 2.721 108 
2.734 108 2.724 108 
3 2.861 X 108 
2.863 < 108 2.8962 x 108 


TABLE II. Frequencies of a square Z-cut tourmaline plate 
with edge a= 1.836 cm (in kilocycles). 


Mode Cale. Obs. 
1 218 215.1 
2 253 247.6 
3 278 291.0 
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2 2 2 2 
T 
8 
—pQ?——k*(yu— 712) =0. (47) 
us 
2 2 8 
——pQ?+—k? (yu — 12) k®—y12 pQ?— 


The two equations (51) refer to the roots Q,? 
and Q;*, respectively. The solutions U2, V2 
and U3, V3 are then: 


U.=C(sin ky — (x/4)) cos kx, 
V2=C(—sin kx+(2/4)) cos ky 
Us=Ccoskx, V3;=Ccosky. 


R. Bechmann® has observed some resonant 
frequencies of square quartz plates of the Z-cut 
type. He finds three frequencies for longitudinal 
modes and expresses them in terms of the fre- 
quency constant N, i.e., the frequency (in kilo- 
cycles) of a plate with @=1 mm so that 
N=10~va. Table I gives a comparison between 
these measurements and the constants N as 
computed by Eqs. (48) and (50). 

In this case, the parameters yx are given by: 


(52) 


and 


Sir Sig 
» 


Sir? — S12? 


For s;; and S$. the values 


Si Sig 

4 12.948 x — 1.690 10-8 
an 

12.79X10-"8 —1.535X10-" 


are given by Bechmann® (according to Voigt’s 
measurements) and W. P. Mason,® respectively. 

The calculated values have been obtained with 
these two sets of elastic constants of Voigt (upper 
figure) and by W. P. Mason (lower figure), 
respectively. Two of the observed frequencies 
agree numerically with the formulas given by 
Bechmann although the expression in terms of 
the sj is quite different. 

V. Petrzilka’ has measured the frequencies of 
a square tourmaline plate of edge a=1.836 cm. 


®R. Bechmann, Zeits. f. Physik 118, 515 (1941). 
. *W. P. Mason, Bell Sys. Tech. J. 22, 178 (1943). 
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Fic. 1. Longitudinal vibration of a thin square plate. 
Mode 1: U,=C cos (x/a)x sin (x/a)y, Vi=Csin (x/a) 
x cos (x/a)y. Dotted lines—crystal at rest. Full lines— 
crystal vibrating. 


The frequencies calculated by Eqs. (48) and (50) 
with Voigt’s values for the elastic moduli of 


.tourmaline are shown in Table II. 


Two of these calculated frequencies are equal 
to those given by Petrzilka. The first frequency 


2 2 
—p2?——k*[yn— v2] 
T 


2 2 
— 712] 


The three eigenvalues of (55) are: pQ®=k? 
X73.9X 10"; k?X89.6X 10"; k? 106.8 X 10". 

The corresponding frequencies expressed in 
terms of the frequency constant N=10~*va are 
shown in Table III. As a comparison the values 
observed by Bechmann® and Wright and 


- Stuart" are given. The second calculated fre- 


quency was not observed by these authors. 
Wright and Stuart studied the mode No. 1 


”R. B. Wright and D. M. Stuart, Bur. Stand. J. 


Research 7, 519 (1931). 
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2 2 2 2 


p2?— k*y22 


is the rigorous value and is, of course, in agree. 
ment with our value. 

The vibrational modes calculated by Eqs. (39) 
and (52) are represented by Figs. 1-3. It can be 
seen that none of them has nodal lines, but only 
nodal points. Petrzilka attributed to the fre. 
quency (2) a mode equal to our zero-approxima. 
tion modes (40) or (41), but did not find the 


-expected nodal line by lycopodium powder ob- 


servations. 

The calculated mode (3), Fig. 3, involves a 
nearly radial motion. It can clearly be identified 
with the picture of the lycopodium powder 


pattern published by Petrzilka. In this mode, 


according to the author, the powder is displaced 
radially and stays at rest only at the center. This 
is also what the theoretical picture predicts. 

We now apply Eqs. (44) to a square Y-cut 
quartz plate, i.e., a plate whose large face is 
normal to the crystallographic Y axis, while the 
edges are parallel to the optic and electric axis, 
respectively. Using Voigt’s values for the elastic 
moduli, we obtain: 


yu= 78.67 4 
2 = 103.08 X10", yie=v26=0. (53) 
¥12 = 12.180 X 10"°. 

The secular equation takes the form: 


8 
— ky =0. (54) 


(N=2.67X10*) by means of powder patterns. 
While the experiments with powder on the 
crystal surface are difficult to interpret, the pat- 
tern created by the air currents from the periph- 


TABLE III. Frequencies of square Y-cut quartz plates 
mm, in kilocycles. 


ode Neale. Bechmann 
1 2.64 x 108 2.68 x 108 2.675X 1 
2 2.90 x 108 

3 3.17XK 108 3.19 108 
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ery of the crystal should give a simple idea of 
- the motion of the boundary. 

In attempting to calculate the ratio A;:A2:A3 
corresponding to the mode No. 1, one meets a 
difficulty. While the discrepancy between the 
observed and calculated values of the frequency 
(2.64 and 2.68) is as small as can be expected 
from the uncertainty of the elastic moduli, the 
ratio A;:A2:Az3 is very sensitive to small varia- 
tions of the root, so that the mode cannot be 
described definitely. As a first approximation, 
we can expect it to be similar to the mode No. 2 
described for the Z-cut plate (Fig. 2). In this 


gential direction and create thereby a wind such 
as is shown in Wright and Stuart’s Fig. 13a. 
The agreement between this picture and the air 
currents to be expected from the theoretical 
figure is quite satisfactory. 


(b) Shear Modes 


A type of modes where the shearing strain x, 
predominates is usually called “shear mode.” 

The zero-order modes will be found by con- 
sidering two extreme cases: If the plate is in- 
finitely extended along the x direction (a—«) 
standing shear waves are exact solutions: 


Mode 2: U:=C[sin (x/a)y—(#/4)] cos (x/a)x, 
V2=CL—sin (x/a)x+(2/4) ] cos (x/a)y. 
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mode, the corners of the plate move in a tan-— 


Fic. 3. Longitudinal vibration of a thin square plate. 


Mode 3: U;=Ccoskx, V;=Ccosky. Dotted lines— 
crystal at rest. Full lines—crystal vibrating. 


If the dimension b is very long, standing shear 
waves along the y direction are solutions: 


1,2,3-++, (56) 
Vpz = C cos phix, , 
where again 
ki=x/a, ke=x/b. 


The modes (55) and (56) will be taken as zero- 
order modes for the study of almost square plates. 
The normalization (22a) requires that 


for all these modes. 
The definition (25a) yields, after elementary 
calculations, the elements of the matrix H: 


There is no interaction between the modes 
with » or g even and the odd modes. Therefore, 
the even modes can be dropped, and in the fol- 
lowing, p and g will be only odd integers. 


Fic. 2. Longitudinal vibration of a thin square plate. 


Dotted lines—crystal at rest. Full lines—crystal vibrating. 


C?=2/ab (57) 


py = Yoo? 


(58) 


H 0 if p or g is even 
(8/ab)-yee if p and are odd, 
Nu=0 for any two different modes. 
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The secular determinant reads: 


As kik, each diagonal element of H is 
quasi-degenerate with one neighboring diagonal 
element. The first step in the solution of (59) is 
to find new linear combinations of the zero- 
order functions so that this degeneracy is 
removed. Considering the case of the square 
plate (ki=k:) and disregarding the coupling 
between non-degenerate functions, we find Eq. 
(24) to be: 


A y2(p2? — yeop?k*) — A py (866/07) =0, 

—A +A py(p2? — yeep?k*) = 0. 
These two equations are compatible if 

Or Apr=—Apy. (61) 


Consequently, the sought linear combinations 
are: 


(60) 


Up+=C cos pkey, up-= —C cos phry, 
where p=1,3,5--- (62) 
Vp+=Ccos pkix, v,-=C cos pkix. 


With these new zero-order functions the 
matrix elements of 7 and N are found to be: 
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8 8 
_ 
ab ab 
8 8 
ab ab 
8 
ab ab 
8 8 


Nx=0 for any two functions (62), 


8 
P+, P+ 2 


8 


The interaction between the modes p+ and 
p— is so small for our case that it can be 
neglected. 


Putting 
b=a(1—e), (64) 


where ¢ is a small number, we find that the 
matrix elements referring to the modes p+ 
become 


H+, p+ = +8) ves, 
H+, a+ = (1/a?)(1+6)8y68 (pq). 


Now, the secular determinant has the form: 


(65) 


a a a 
8 66 8 
po? ——(1+«)(9x? +8) + €) 
a va a (66) 
8 
— yall +6) 
a a a 


In 


or, 
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t 
If we pu r= pQ*a?/8yee(1+€) 


and divide each row of (66) by 8ye(1+.«)/a*, the determinant is simplified: 


The roots of (68) are approximated by the by Hight and Willard." They studied quartz 
Schrédinger pec formula: plates whose edge was parallel to the crystal- 
lographic X axis. In this case 


Sis = 


Equation (71) gives, for p=1, 


The sum in this expression can be readily 
calculated. As p and g are odd numbers, the sub- 
stitutions p= 2n+1 and g=2m+1 give: n=—(—) 1.28, 


a . while Hight and Willard’s observations give 


mao (2m+1)?—(2n+1)? 


1 ( -—(—) 1.25, 


, so that the difference between theory and experi- 
where m takes all integer values except m=n. ment is 2.5 percent. Later observations of 


Inspection of this expression shows that the sum — pahmann?® agree within 1 percent with those of 
Hight and Willard. 

' 1 1 For the resultant modes, the Rayleigh- 

4(2n+1) at Schrédinger perturbation formula yields from 


Thus, (69) yields: Eq. (68): 


Up+ 


V.=9,4+-— 
In view of (67), the frequencies are: q?—p? 


or, in view of (64): Ur=coshiy +=, 
q*l 


k (73) 
cos x 


Systematic observations were made on the “j¢ ¢ Hight and G. W. Willard, Proc. ILR.E. 25, 549 
lowest shear mode of nearly square quartz plates (1937). 
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Fic, 4. Shear vibration of a'thin square plate. Observa- 
tions by R. A. Sykes. Full and dotted lines represent the 
contour deformation in two opposite phases. 


This result can be compared with observations 
on the mode of a square quartz plate vibrating in 
the lowest shear mode. Figure 4 shows the dis- 
torted edges in opposite phases, one a solid curve 
and the other a dotted curve, according to 
microscopic observations by Sykes.” Figure 5 
shows one phase (full line) of the motion cal- 
culated by (73). Only the first terms of (73) have 
been used for the drawing because the supple- 
mentary terms are small. The agreement between 
observation and theory is satisfactory. 

Theoretical considerations of the shear mode 
by Mason! and Bechmann"™ led to a mode in 
which the corners of the plate should be at rest. 


12 R. A. Sykes, Bell Sys. Tech. J. 23, 52 (1944). 
18 R. Bechmann, Zeits. f. Physik 117 and 118 (1941). 
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Fic. 5. Shear vibration of a thin square plate (calculated), 
U=cos (x/a)y, V=cos (x/a)x. Dotted line—crystal at 
rest. Full line—crystal vibrating. 


Their solution apparently bears no resemblance 
to the experimental picture. 

All comparisons with experiment for the longi- 
tudinal as well as shear modes were made in 
cases where the parameters and vanish. 
Where these parameters are not negligibly small, 
the longitudinal and shear modes would be 
coupled and could not be treated separately as is 
done here. 

The author is indebted to Dr. F. Beer for 
advice and help. He wishes to thank Mr. A. H. 
Schlesinger for his help in numerical calculation. 
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The propagation of e-m waves is investigated: (1) in a rectangular metal tube half-filled 
longitudinally with a dielectric; (2) in a cylindrical dielectric guide of radius a, surrounded by 
a coaxial metal tube of radius b. In (1) the waves are of the longitudinal-section type; for high ° 
frequencies they are confined to the medium of higher dielectric constant. The problem is 
discussed also with the model of the criss-cross component waves. In (2) the waves are linear 
combinations of E and H waves. The case of axial symmetry is considered, in which there are 
simple E and H waves. For Epa critical wave-length exists, depending only on the dimensions 
of the inner dielectric, and on the dielectric constants, below which the system behaves more 
or less as a dielectric guide in free space, and above which as an ordinary hollow tube. For Ho 
this critical wave-length depends also on the ratio b/a. The case in which the external medium 
has a higher dielectric constant is also briefly investigated. 


1 


propagation of electromagnetic waves in hollow 


HE purpose of this note is to investigate metal tubes, the interior of which is filled with 
from a purely theoretical point of view the two different dielectrics, the distribution in all 
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cross sections being constant. If one of the media 
is air, the problem may be called that of a dielec- 


- tric guide contained in a metal tube. The features 


of the phenomenon are similar in all problems of 
this kind. Two examples will be treated: (a) a 
rectangular tube with two dielectrics, each half 
filling it longitudinally ; (b) a cylindrical dielectric 
guide of radius a, enclosed by a coaxial cylin- 
sical metal tube of radius 6>a.! 

We consider a hollow metal tube with in- 
finitely conducting walls, of rectangular cross 
section, with axis parallel to the Z axis. The 
dimensions of the tube are a in the X direction 
and b in the’ Y direction. The whole length of the 
tube is filled, from X =0 to X = 4a, with a non- 
absorbing dielectric of constant ¢«,, and from 
X=}a to X=a with a non-absorbing dielectric 
of constant «, (Fig. 1), the magnetic perme- 
ability of both media being taken equal to that 
of empty space. We shall take e2:>,. 

The differential equation satisfied by the 
longitudinal component of the electric field for 
E waves (or by the longitudinal component of 
the magnetic field for H waves) is: 

@Ez 


4 +(h? 1 
aye (h?+euw*) Ez (1) 


Here m.k.s. rationalized units are used; w=2r 
‘frequency, and it is supposed that the de- 


pendence of each field component on Z and ¢ is 
given by a factor e~*#+**, ¢ has the value ¢«; in 


. the first medium’ and « in the second medium. 


We may further assume that the dependence on 
Y is the same as for ordinary tubes, namely, that 
we can write: 


Ez=f(X) sin (= 1,2---. (2) 
b 
f(X) satisfies then the equation 


+ (i+ joo. (3) 


For the bibliogra on the propagation of e-m waves 
in metal tubes, graphy ont L. Lamont, Wave Guides 


(Methuen & Com i niga: articles by E. U. Conten 
in Rev. Mod. Phys. (1943). 


Writing 
m=—, 


ko? 
we obtain the solution 


S(X) sin cos kiX 
in the first medium, 
f(X) sin +B, cos k2X 
in the second medium. 
The boundary conditions at the walls give 


B,=0; B,=—Az tan kya, (5) 
so that 
f(X) sin (0<X <4a), 
f(X) = —(A2/cos sin ke(a—X) (6) 
(4a<X <a). 
xs 


2 
Fic. 1. System of coordinates and dimensions of the tube 
studied in the first problem. 
The expressions of the field components are: 
(f'(X) =df/dX) 


h OE h 
Ex=- sin m +iet, 
h?+euw? k*-+-m? 
h 0Ez 
OY 
h 
cos m 
k?+m? 


Ez=f(X) sin 


we 


h?+-euw® 


iwe OEz 
Hy=-— 
= sin m 
Hz=0, 


(4) 


=———_f(X) m cos m Ye" 
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where we must insert for each medium the cor- 
responding values of e, k, and f(X). 
Similarly for H waves, writing 
Hz=¢(X) cos mY, 


we find: 
¢(X)=C,; cos in the first medium, 


C2 


cos k2(a— X) (6’) 


o(X)= 
cos koa 


in the second medium; 


and the expressions for the field components are: 


Ex= ¢o(X)m sin m Ye! + 
24m? 
Ey= ¢’(X) cos m 
k2+m? 


Ez=0, 
(7) 
Hxy=- ¢’(X) cos m 
k?+m? 
h 
Hy = ~o(X)m sin m 
2 


Hz=¢(X) cos m 


We have now to determine the constants A, Ao, 
Ci, C2, in order to satisfy the boundary conditions 
at the separation of the two dielectrics. This, 
however, is found to be impossible with a simple 
E or H wave (with the exception of the 7 wave 
for m=0). The waves in our tube must then be 
a linear combination of E and H waves. Including 
the coefficients of this linear combination in the 
constants A;, As, Ci, Cs, and introducing the 
(dimensionless) quantities 


x= }k,a, y= koa, (8) 


the continuity of the tangential components of 
E and H at the boundary gives the following 
system of homogeneous equations: 


sin y 


A,sinx=—Ay, 


cos 
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(A ,hm+ Cyiwuk) sin x 


ky?+m? 
(A 2hm + Criwuk,) sin y 


(ko? +m?) cos 2y 
cos 2y 
(—A Cihm) cos x 
(—A + Cohm) cos y 
cos2y 


(9) 


In order to obtain a solution different from 0 for 
Aj, Ao, Ci, Co, the determinant of the coefficients 
must vanish. 

Expansion of this determinant leads, if we 
discard trivial solutions, to the two transcen- 
dental equations 


tan y/y= —tan x/x, (10) 
y tan y= — (€o/e,)x tan x. (11) 


As we shall see, different waves correspond to 
these two equations. Equations (10) and (11) are 
easily solved numerically or graphically, and 
curves can be plotted of y against x. In the case 
of Eq. (11) different curves are obtained for dif- 
ferent values of the ratio ¢2/e:. If we had con- 
sidered media with magnetic permeability dif- 
ferent from that of empty space, Eq. (10) would 
have contained the ratio pe2/m. 

x and y are further connected by the equation, 
deduced from (4): 


€2— €) wa? €9— €) 


y?—x? (12) 


€0 4c? €0 
(€9 =dielectric constant of vacuum), 


_ wa ma (13) 


2¢ 
(Ao being the wave-length in vacuum), which, 
together with either (10) or (11), completely de- 
termines them for a given frequency. When * 
and y are known, the wave-length A in the guide 
is found from (4): 


h?= ew? = ke? +m? — 


show! 
ordin 


Fic. 
the firs 
with 
section. 
on th 
param 

The 
=fe/g 
We se 
consta 
single 
one-di 
freque 
case 
single 


Witl 


ELECTROMAGNETIC WAVES IN METAL TUBES 


or writing 


(14) 
(15) 


g=}a8=na/d 
p=}am=}nal/b 1=0,1,2---, 


and 


we have 


showing that a cut-off frequency exists as in 
ordinary tubes. This, and the ratio A/a, depend 


2 3 4 £ 5 

Fic. 2. Values of x* and y¥* at different frequencies for 
the first two I-s Ey,1, and the first three I-s Hy,1 waves 
with e:/¢:=2. Cut-off line refers to tube of square cross 
section. 


on the shape of the tube, which affects the 
parameter p. 

The phase velocity of the waves is v=w/8 
=fc/g, and the group velocity u=dw/dB =gc/f. 
We see that in this problem x and y are not 
constant for a certain wave, as is the case for the 
single quantity which appears in their place in a 
one-diclectric tube, but are functions of the 
frequency and of the dielectric constants. In the 
case of Eq. (10), however, they depend on the 
single variable 


ra 


3 


With fixed Xo, infinite pairs of values ya(o), - 


‘(No) are found satisfying (12) and either (10) 
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or (11). To each pair of values corresponds a 
particular wave that can be labeled with the 
order of the corresponding root. When e—<«; 
(or Ao), both x, and y, approach nx/2, the 
smallest value of m being 1 for Eq. (10) and 0 for 
Eq. (11), so that we shall label ‘‘one’”’ the first 
wave corresponding to (10) and ‘“‘zero”’ the first 
wave corresponding to (11). 

An examination of the fundamental equations 
shows that when w increases, y, increases, while 
x, decreases, until, for a certain w, x,=0. When 
this happens, y in the case of Eq. (10) has a 
value which is a root of y cot y= —1, namely, 


y=0.645, 1.564, 2.54, (n—4)x, 


and in the case of (11) has a value which is a 
root of the equation tan y=0, or y=nr. If the 
frequency is further increased, x, becomes purely 
imaginary. Equations (10) and (11) become: 


tan y tanh |x| 


(10’) 


€ 
y tan y= +—|x! tanh |x}, (11’) 


€1 


where the modulus of x is indicated by |X}. 
Then |x| increases steadily for increasing w, and, 
for becomes proportional to w: 


for all n, 


while y, approaches the value mz in the case of 
(10’) and the value (n+4)z in the case of (11’). 

These results are shown in Fig. 2. y*/x* is 
plotted as abscissa and x*/x* as ordinate. The 
diagonal lines give 

Go a*w? 

f €0 4c? 
The chain curves give the first and second solu- 
tion of (10), (10’) and are independent of the 
values of the dielectric constants. The broken 
line curves give the first, second, and third solu- 
tion of (11), (11’), for a ratio ¢:/e,;=2. The curves 
for other values of this ratio are easily deduced 
from them. 

The other line drawn on this graph is the cut- 
off line for a particular value of p, namely, 
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p=n/2, corresponding to /=1 for a tube of 
square cross section. Along this line, if equation 


g=0, or 


the wave-length in the guide becomes infinite. 
In the region above it no waves are propagated. 
Lines parallel to the cut-off line, represented by 
(16) for arbitrary values of g, would give the 
values of f, x, y, corresponding to a given wave- 
length in the guide. 

In the region of propagation we have gener- 
ally from (16): 


For small A» we have seen that 
€2— 


and we can neglect p’, so that, at high fre- 
quencies, \=Ao(€o/e2)', the value corresponding 
to propagation in the dielectric of higher con- 
stant, for both kinds of waves. 

To see what happens when the frequency 
decreases, let us consider the case €;=€9; then 


\ remains <A» when x is imaginary and |x| >); 
otherwise \ is >A o, and becomes infinite for a 
value of A» which is higher the higher the dielec- 
tric constant of the other medium. As usual the 
wave-length and the cut-off frequency increase 
for increasing order of the root. 

The general case is analogous. 

It is evident from Fig. 2 that for x imaginary 
the second kind of waves has, at the same fre- 
quency, a smaller value of |x| than the first kind 
of waves of the same order, so that, for x 
imaginary, the first kind of waves has a smaller 
wave-length than the second. For x real we see 
from the figure that the two curves cross several 
times generally, so that the wave-length in the 
first case may be larger than, equal to, or smaller 
than the wave-length in the second. 


4 


We indicate now by x and y a pair of values of 
these quantities satisfying either (10) and (12), 
or (11) and (12). Inserting them in the relations 
(9), we find that the ratios of the constants 4,, 
Ce, to A; are: 

Az sin x 


cos 2 
A, sin y (18) 


Ci €) 
for Eq. (10); 
A, 


(20) 


() for Eq. (11); 


Indicating now by A in the case of the first 
equation, and by B in the case of the second 
equation the only constant which is left, and 
which depends on the experimental conditions in 
which the tube is excited, we obtain the following 
expressions for the field components: 

In the case of Eqs. (10), (10) we have in the 
first medium: 


Ex =0, 


h 
Ey = —A— sin cos m 
m 


Ez=A sin sin 
h?—m? 


Hx=A -sin cos (21) 
TM 


=A— cos k,X sin 


hk 
Hz=A—— cos kiX cos 


In the second medium A _ is_ replaced by 
A sin x/sin y, ki by ks, and by 


Moreo 
Hp. In 


first m 


Ex 
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taste |. Ratio of amplitude in medium of higher ¢ to 
amplitude in medium of lower e. 


For Ey; Ez; Hx For Dx; Hy; Hz 


sin x ysinx 


sin y xsiny 


sin x xsinx 


sin y ysiny 


Moreover, there is a change of sign of Hy and 
Hy. In the case of Eqs. (11), (11) we have in the 
frst medium: 
m? 
cos k,X sin m 
1 


m 
Ey= sin k,X cos m Ye~*+ 


sin sin m 


Hy =0, 


Hy = —B— cos k,X sin m 
1 
meywt Pte 
Hz= —B—— cos k,X cos m 
1 


In the second medium B is replaced by 
Bsin x/sin yy, €1 by €2, ky by ko, and RX by 
bia—X). There is a change of sign of Ex, Hy, 
Hz. In the region in which x (or k,) is imaginary, 
we have 7 sinh |x| instead of sin x, and cosh |x| 
instead of cos x. 

The field described by (21) is characterized by 
the absence of Ex, and that described by (22) by 
the absence of Hx. Such waves are called longi- 
ludinal-section electric and magnetic waves, 
respectively, l-s Z waves and |-s H-waves),? and, 
when ¢=€,, they become those linear combina- 
tions of the ordinary E and H waves of the same 
order, which have a vanishing Ex and Hx, re- 
sectively. (For /=0 or n=0, the wave is simply 
an IT wave.) We label them with the correspond- 
ing values of m and J, namely, I-s Ey, 1, Hn, 
For I-s E waves 1=0 is possible, and the first 
Wave (n=1) reduces, for to Hio0; for 
'sH waves /=1 is the first possible value, and 


"H. Buchholz, E.N.T. 16, 73 (1939). 
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the first wave (n=0) reduces, for to Ho, 1. 
This notation is used to identify the curves of 
Fig. 2, in which the second index is taken equal 
to 1, so that the cut-off line applies to a tube of 
square cross section. 

As ) is generally different for the two kinds of 
waves, we cannot, except in particular cases, 
form linear combinations of them that will 
produce ordinary E or H waves. In both media 
all the field components are sine functions of X, 
but with different amplitudes, so as to satisfy 
the boundary conditions. The ratios of the am- 
plitudes in the two media are given in Table I. 

When «x is real, all these ratios remain finite, 
as x and y reach together values multiple of z, 
the limit of sin x/sin y being finite. When x=0 
the ratios are not immediately given by Table I, 
but it is easily found that the ratio for Dx, Dy, Hz 
remains finite in both case, these components 
being independent of X in the first medium. For 
Ey, Ez, Hx the ratio is also finite for l-s E waves, 
but is infinite for |-s H waves, these components 
being then entirely confined to the second 
medium. 

When x is imaginary, we must write 7 sinh | x| 
instead of sin x and then when |x| increases, all 
the ratios approach infinity, namely, for large 
values of |x| (high frequencies) all components 
are confined to the medium of higher dielectric 
constant; we have already seen that the wave- 
length approaches then the value corresponding 
to this medium. 

Figure 3 shows the dependence on X.of the 
first three components for different values of f’ 
in the case of I-s E;,,,, Fig. 4 in the case of 
l-s H;,, for e2,/e,=2, and Fig. 5 gives the same 
for l-s Ho,; at a fixed f’ (f’ =3.33) and for dif- 
ferent values of €:/e,;=«. All ordinates are in 
arbitrary scales. 

From Eqs. (21), (22) the E and H lines are 
easily drawn: They present a pattern which is 
similar to that of i-s waves in one-dielectric 
tubes, but somewhat distorted. For instance the 
E lines for l-s Ho, ; (coinciding with Ho,; in a one- 
dielectric tube) are not straight lines in the 
X-Y plane. 


It is interesting to examine the above results 
by use of the picture of the component criss-cross 


es of 
(12), 
As, |-s E waves 

(18) |-s H waves 
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waves.’ They are plane transverse waves, criss- ponent waves in the first medium, but only waves 
crossing along the tube, and moving in each propagated parallel to the surface of separation, 
medium with the velocity that they would have and with an amplitude falling off exponentially 
in the same medium unrestricted by the tube. (in a region of the order of the wave-length) in a 
The waves in a rectangular tube are generally the — direction normal to this surface. Hence the de- 
pendence on X of all the field components in- 
, \ volves the hyperbolic functions. Each of the 


feo / \ 
i \ component waves falls off in the first medium 
° 3 according to a factor e~®/®!='X, This is in agree- 
A ment with the well-known expression for the 
: ay \ light entering the medium of lower refractive 
aig index in total reflection; in fact, this has a factor! 


exp [—(w/c)(n? sin? i—1)4X ] 


443 \ es 
(where 7 is our 02, m?=€2/¢o, and and we 
have 


Fic. 3. Ey, (Ez, Hx) as function of X at different wthie 2 
frequencies for £;,:. | 1| 


resultant of 4 such waves, in some cases only of 2. feo / 
/ 


The angle @ that the propagation vector of the f \ 

component waves makes with the X axis is (for ee 

x real) given by A 

cos 0; = (€0/€1)!(cki/w) = (€0/€1)*(x/f) 


in the first medium, 
COS = (€0/€2) = (€0 /€2) 


in the second medium. 


= 


6, and 6 increase with increasing frequency, 
always satisfying the law of refraction 


sin 4; 


sin 


(as until, for x=0, 
reaches the value 1/2. 6:(e=0) is then the Fic. 4. Ey, (Ez), as function of X at different 
critical angle. This explains why x becomes frequencies for I-s Hi: and ¢2:/e= 

imaginary when the frequency is further in- 
crease: There can be no more ordinary co1n- 


When the frequency is high and the falling off 
is rapid, multiple reflections are negligible, and 


3L. Brillouin, Rev. gén. de l’Elec. 40, 227 (1936). Lees 
Page and N. I. Adams, Phys. Rev. 52, 647 (1937). L. 4 See, for example, the article by W. Kénig in Handbuch 


Pincherle, Phil. Mag. (in print). der Physik, Vol. XX, p. 231. 
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hen also the ratio of the amplitudes of the com- 
gonent Waves at the surface of separation coin- 
sides with the usual expression. In fact, this ratio 
i, our tube for Ey and |-s E waves is 


r=(sin ye'*!)/(sinh |x!). 


1G. 5. Ey, (Ez), as function of X for l-s Ho,1, with 
‘= 3.33, and different values of «= e2/e;. 


\ow for high frequencies |x| is large, and the 
ratio Is 


2 tan y 2y 


y? 


tan’? y 


r=2 sin 


Using (10°) with tanh |x| =1, we have then 
2y 2y 2n cost 

r=— = = 
fl (e2/eo)—1]! (m?—1)! 


the well-known expression. 

\Ve consider now the propagation of e-m waves 
ina cylindrical dielectric guide of radius a and 
constant €;, surrounded by a medium of dielectric 
constant €2, bounded by a coaxial metal tube of 
radius b (Fig. 6). We use cylindrical coordinates 
'.¢, 2. We assume the conductivity of the two 
dielectrics to be zero, and that of the metal to be 
infinite; we take the magnetic permeability of 
both dielectrics to be that of empty space. This 
problem reduces to the ordinary hollow tube for 
=e (or a=0, or b=a), and reduces to the 


problem of the dielectric guide’ for b= ©. We 
try to write solutions of Maxwell's equations of 
the E and H wave type, namely, with either an 
E; or an H,. Writing, as before, =k 
= ko? — jw", we consider the solution 


N (kr) 


where J, and N, are the Bessel and Neumann 
functions of order n,® and the field components 
are determined, in the case of E waves, by: 

1 au eu 
=—; E,=-——; E,=—+ew’ 

r 

iwe OU 
H,=——-; H,=—-we—; H,=0; 

r de or 


cos y=0,1,2--- (25) 


E, 
(26) 


and in the case of HT waves by: 


H,= ; H,=- ; H,=—+ eww? 
Ordz r 


We have two different solutions in the two 
media. In the first (inner) medium we cannot 
have the Neumann function, as it becomes in- 
finite for r=0; the solution in this medium will 
be: 

U,= CJ a(kir) cos (28) 


In the second (outer) medium the solution is: 
U2=[AJa(kor) + BN, (kor) ] cos nge™*t*, (29) 


At the walls of the tube EZ, and E, must vanish; 


Fic. 6. The tube studied in the second problem. 


5D. Hondros and P. Debye, Ann. d. Physik 32, 465 
(1910). J. R. Carson, S. P. Mead, and S. A. Schelkunoff, 
Bell Sys. Teck. J. 15, 310 (1937). 

6 Definition of Bessel, Neumann, and Hankel functions 
as in E. Jahnke and F. Emde, Tables of Functions (Teubner, 
Leipzig, 1933), pp. 194-200. 
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this gives, in the case of E waves, 


= (30) 
and in the case of H waves 
(31) 


the dash indicating differentiation with respect 
to the argument. We have further to satisfy the 
boundary conditions at r=a, namely, the con- 
tinuity of E,, E., H,, H.. These, as in Section 2, 
cannot generally be satisfied by a single E or H 
wave, and a linear combination must be intro- 
duced. The calculation can be carried out, the 
equation for the proper values of k; and ke 
deduced, and solutions found in rather a similar 
way as that used for dielectric guides by Carson, 
Mead, and Schelkunoff; but, as the solutions for 
different n’s are not fundamentally different, it 
will be sufficient to discuss the problem for the 
simpler case = 0 (circularly symmetrical waves), 
in which all the boundary conditions can be 
satisfied by a simple E or H wave. The funda- 
mental equation is in this case 


Ji(x) _@ Jo(my)Ni(y) — Ji(y) No(my) (32) 


xJo(x) Jo(my)No(y) —Jo(y) No(my) 


for the Ey wave, and 


Ni(my) — Ji(my) Ni(y) 1 
xJo(x) Jo(y)Ni(my) —Ji(my) No(y) 


for the Hy wave, where x=k,a, y=kea, and 
m=b/a. If we write 


€omowa = (24a) /Xo 


(Ao=wave-length in vacuum), (13’) 
g=(ah)/t=(24a)/d 
(\=wave-length in the guide), (14’) 


x and y are connected, as in the first problem, by 
the additional relation 


y?=x?— (€;—€2/€0) =x? — 


which, together with either (32) or (33), com- 
pletely determines them for all frequencies; and 


(12’) 
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the wave-length in the guide is found from the 
relation 


g°= (€2/e0) f?—x*, (35 


Equations (32) and (33) have to be solved ny 
merically. Tables for the left-hand side 


Ji (x) 
¢(x) = 
xJo(x) 


Each ¢ 
the tuk 


and for the right-hand side (for E waves) 
Jo(my) Ni(y) — Jily) No(my) 
y Jo(my) Noy) — Jo(y) No(my) 


have to be constructed for various values of m. 
Then the fundamental equation becomes 


¥(y) 


€2 
¢(x) = 


and the relation between x and y is immediately 
found from the tables. We have a series of x-y 
curves, each of which gives a solution. In the 
case of H waves the procedure is simpler, as then} 


o(x)=H(y), 


and the x—y curves depend only on m and not 


Fic. 8. 
guide ani 
values of 


Hy. 1, Eo, 2 
Ev, Ho, 
reduce te 
and Deb 
previous 

sponding 
constant 
and vanis 
frequenci 


| 


+6 


negative 
— convenier 
‘ Hankel { 
ms! mel \ 
ae Equation 
-6 i\ same wit 
either m- 
Fic. 7. Values of x* and », at different frequencies and 4 islarge, J 
for different m, for Eo,:. Dotted curves «:/e:=1.1; chain 
curves ¢:/¢:=5; broken line curves ¢:/e2:= 81. instead o 


on the ratio €:/e;. (In media with these 
would depend on 

If the hyperbola given by (12) for the required § 
value of the frequency is also ‘plotted, the inter- § 
cepts of this hyperbola with the other curves @ 


ive the values of x and y at the given frequency 
nd for each solution of the fundamental equa- 


non. 


7 


Each curve corresponds to a different wave 
athe tube, which can be labeled Eo,1, Eo,2, **-, 


6 6 


Fic. 8. The relation between the wave-length in the 
wide A and the wave-length in vacuum Xo, for different 
values of €;/e9 and m, in the case of Eg; with e:= 0. 


Hy ++ +, and which reduces to the ordinary 
for €eg=e,. When b= the E waves 
reduce to the waves investigated by Hondros 
aud Debye in a dielectric guide. As in the 
previous problem the parameter, say y, corre- 
sponding to the medium with smaller dielectric 
constant decreases when the frequency increases, 
and vanishes for a certain value of it. For higher 
frequencies the curve y—x is continued for 
iegative values of y*. When y is imaginary it is 
convenient to substitute in (29) one of the 
Hankel functions, say H,,, for Neumann's. 
Equations (32) and (33) remain formally the 
ame with H, written instead of N,. When 
tither m—+0, or the frequency is high, and |y| 
slarge, H, (my) can be neglected, and we have 
instead of (32) and (33): 
Ji(x) 

xJo(x) ey Ho (y) 
Ji(x) 


xJo(x) Ho (y) 


for Eo, (36) 


for Ho. (37) 
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Equation (36) is identical with Hondros and 
Debye’s equation for dielectric guides. 

We see that it is then more convenient to plot 
curves of y’ against x*, instead of curves of y 
against x, so that both cases, y real and y 
imaginary, may be inciuded in the same diagram. 
Also in such a diagram the curves representing 
(12) are straight lines and can be conveniently 
drawn. Figure 7 gives such a diagram, in which 
the y?—x* curves are shown for Eo,:, for the 
values 1,1 (dotted curves), 5 (chain curves), 81 
(broken line curves) of the ratio €:/e2, and for 
m=1 (ordinary hollow tube), m=1,1, and m= © 
(dielectric guide). 

The wave-length in the guide is found from the 
diagram by means of (35). The equations of the 
cut-off lines are y?=(e2/e,)x*. Figure 8 gives the 
relation between A and Xo in the case of Eo,;, for 
€:/eo=2 and 5, being equal to and for 
m=1 and 2. The curves for m>2 differ little 
from the curves for m=2 below the critical 
point. Above this point they remain nearer to 
the line A=Xo, representing propagation in 
vacuum, than the curves for m=2, and they 
actually coincide with this line for m= «. We 
have 


> €0 > 
r=) (x2. if the second medium is air), 
€2 < 


according as 


do > 
= 2.613 


and independently of the radius of the metal 
tube. This wave-length 


(38) 


€0 


can be called the critical wave-length; it cor- 
responds to y=0, and to the component waves 
(of which there is an infinite number in a circular 
guide) falling on the surface of separation under 
the critical angle. 

For \o<A, our tube behaves more as a dielec- 
tric guide than as a hollow tube: The curves 
A—Ao for all values of m do not differ much from 
the curve referring to m= ©; \ increases when 
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TABLE II. Values of x, y, \ at different frequencies for Eo,; and €; > 2. 


CHERLE 


x 


First root of the equation:* 
No(mx) e:Ji(x)No(x) — o(x) Ni(x) 
J o(mx) (€1— €2) Jo(x) Ji (x) 


Static field 


Ao>Ac Real and increasing Real and decreasing A> 

decreasing Decreasing with increasing m 
Ao=Ac 2.4048 0 = (No/a) (€o/€2)* = 2.613(4/e—1)! 
do <A- de- Real and increasing, approaching the Imaginary and ap-  <Aol€o/e2)# 

creasing value 3.8317 proaching the value Approaching the value \ = o(€/«:)# 

€1— €2 increasing with increasing m 
14.682 — — 
€ 

Limit case 2.4048/m at all frequencies =x a\? a\?« 

€2 -) -(<) ——0.1465 

Nos €0 


the radius of the metal tube increases, in con- 
tradiction to the case of one-dielectric tubes. 

For \»o>A-, on the other hand, our tube 
behaves as an ordinary tube. A decreases when 
m increases, and increases rapidly with Xo be- 
coming infinite at a value of XA» which is the 
larger, the larger the value of «. When m—>~, 
\ remains always equal to Xo in this region, and 
the propagation is accompanied by radiation 
from the sides of the dielectric guide. 

For H waves the critical wave-length is not 
independent of m. For instance, for Ho,1, ./a is, 
as before, equal to 2.613 [(€1:—€2)/eo | for m= ~, 
but decreases with decreasing m, and reaches the 
value 


\./a=1.640[ (€:—€2) 


for m=1. For H waves X decreases with in- 
creasing m even when \o<A,. Table II for Eo,1, 
and Table III for 7,;, summarize the variations 
that x, y, \ undergo when the frequency in- 
creases from a static field to very high values. 
No experimental results are known for the 
wave-length in such a guide other than the ex- 
periments of Schriever’ and Ka§Spar* on water 
guides in air, without metal tube (m= ~); they 
confirm the theory for \»<A, and the fact that 


70. Schriever, Ann. d. Physik 63, 645 (1920). 
8 E. KaSpar, Ann. d. Physik 32, 353 (1938). 


* This root (for e2=e0) is 2.4048 when e1 =e and m =1; it decreases both for increasing e:1 and increasing m, approaching the value 0, both for 


A=Xpo at the critical frequency. In this case we 
expect that for \)>A., A should always be equal 
to A». Schriever’s experiments show instead that 
A>Ao, and his results are consistent with what 
would be expected, had there been a metal tube 
with m~5. These results, however, have been 
criticized by KaSpar, whose investigation does 
not extend to the region Ao >A-. The expressions 
for the field components are: Eo wave; first 
medium: 


hx r 
E,=A 
a a 
r 
a? a 


x r 
H,=A 
a a 
E» wave; second medium: 


hy r 
Nolmy)( 
a a 


—Jo(my)N 


Static 


ho>A 
cre 


Limit 


with 


(A, d 
The « 
nary 
instea 
medit 
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TABLE III. Values of x, y, A at different frequencies, for Ho,; and €, >. 


x » 
Static field  3.8317/m 3.8317/m 


ho> Ac de- Real and increasing Real and decreasing 
creasing increasing with increasing m 


First root of the equation* 0 d/a = (Ao/a) (€o/€2)# = (24 /x)(e1/e2—1)4 
Ji(x)xJ o(x) = 4(1—m?*) 


ho de- Real and increasing, approaching Imaginary, approach- <Aol(eo/e2)# 
creasing the value 3.8317 ing the value approaching the value \ =)o(¢o/e:)# 
€1:— 
14.682 — — 
Ao 


Limit case 3.8317/m at all frequencies =x 


= €3 


* This root is 3.8317 for m =1, decreases for increasing m, and is 2.4048 for m= ©. 


r\ Ho wave; second medium: 
E,=A o(my) J y-) 
a y r 
(40) E,= 9") 
a a 


Jotmy) Nol 9) 


r 
a a 


—Ji(my)N (27) 


h r 
H, = 
a a 


x? Jo(x) 
(41) 
y? Jo(y) No(my) — Jo(my) No(y) H,=B (my) 


4,=A 


(A, depending on the experimental conditions). ‘ 
The other components vanish. When y is imagi- —J;(my) Ni) foie 
nary we have in (40) one of the Hankel functions a 

instead of the Neumann function. H» wave; first x? To(x) 

medium: B.=B 


Jo(y)Ni(my) — Ji(my) No(y) 
x) As in the previous problem, when y is imaginary 
the waves are confined to the inner medium, the 
hx r ; more so, the higher the frequency. For E waves, 
H,=B,—J moreover, and £, are confined to the medium 
” ” of higher dielectric constant also at the critical 
x? ‘ point. The way in which the dependence on r 
H,=B,;— of the field components varies with the frequency 
a’ a is similar to that of the other problem. 
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TABLE IV. Values of x, y, \ at different frequencies for Eo,; and ¢:> «1. 


First root of the equation:* 

No(my) 0(¥)Ni(y) — No(y) 
Tomy) (exe) Jo(9) 
First root of the equation:** 

Jo(my) No(y) — Jo(y) No(my) "2 


First root of the equation:*** 


-=f — 

a €1 a y¥\a ) 
No(my) No(y) 


= A—(€0/€2)4Xo 
Jo(my) Jily) 


* This root, for «: =e, is 2.4048 when ez =e, m=1, and decreases for increasing m approaching the value 0. For a given m it increases with 


increasing e: approaching, for e:—+ ©, a finite limit, always <2.4048. 
** This root (for ¢1=e0) is infini te for m=1 for all values of es, and decreases when m on to the value 0 for m= ©. For a given m it 


increases with increasing ¢2, approaching a finite limit for es ©. It has a finite limit also for exe 


*** The roots of this equation are in 
=0.514 for m=5; =2.4048/m for m— 


lependent of the dielectric constants. The first root is 0 for m=1; itis = 3.4 for m =1.5; =1.80 for m=); 


TABLE V. Values of x, y, \ at different frequencies for Ho; and ¢:> «1. 


x 


3.8317/m 3.8317/m 


First root of the equation:t 
1 Ji(y) Ni(my) — Ji(my) Ni(y) 1 


Joly) Ni(my) —Ji(my) No(y) “2 


First root of the equation:ftf 


J\(my) 


Ji(y) 


—0 


Ni(my) Ni(y) 


A—>(€0/€2) ro 


+ This root is independent of the dielectric constants. It is infinite for m =1; =6.32 for m=1.5; =3.17 for m=2; =0.847 for m=5; =3.8317/m 


for ©. 


tt This root, independent of the dielectric constants, is © for m=1; 


In the frequency region in which hollow tubes 
are used, almost all dielectrics are absorbing. 
The problem should then be treated by intro- 


ducing a complex o1 being 
the conductivity of the guide. The solution of the 
fundamental equations would then require the 
use of Bessel functions of complex argument, a 
cumbersome procedure, as complete tables of 
these functions are not available. The introduc- 
tion of a conductivity does not, however, lead to 
any important change in the features of the 
phenomenon, except those already known for 
ordinary hollow tubes. 


=3.95 for m=1.5; =2.275 for m=2; =0.710 for m=5; =0 form=~, 


9 


If the outer medium has the higher dielectric 
constant, for high frequencies the waves are con- 
fined to the space between the two cylinders of 
radii a and b. This represents a new kind of 
guide, limited on the outside by the metal of the 
tube, and on the inside by the medium of smaller 
dielectric constant. 

Tables IV and V, analogous to Tables I! and 
III, give for these waves the values of x, y, \ at 
different frequencies for Eo, ,and Hp, 1, respectively. 

The relation between \ and Xo is similar to the 
Case > €2. 
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The equation of state of monolayers is worked out from exactly the same point of view as the 
statistical theory of liquids previously given by the author. The formulae obtained should be 
applicable as well to gaseous films as to condensed films, though not to the domain of transition 
between the two, since a change in the potential must be supposed to occur in the transition. 
The resulting formulae are compared with two groups of observations of gaseous films. The 
first group concerns diesters which form insoluble films where the surface pressure has been 
observed directly. The second group concerns adsorbed films of fatty acids where the surface 
pressure has been deduced indirectly from measurements on the surface tension. The agreement 
is satisfactory in both groups. For condensed films the potential can be determined in such a 
way that the density and the surface vapor pressure (where this is known) come out correctly. 
The compressibility is then of the right order of magnitude. 


I. INTRODUCTION 


T is well known that in some cases surface 

films, both the insoluble and the adsorbed 
sind, approach with a reasonable accuracy the 
ehavior of “‘two-dimensional gases.”’ ' However, 
most measurements of surface pressure, direct or 
indirect, refer to conditions where the actual area 
of the molecules and their lateral adhesion cannot 
‘« disregarded. Therefore, two-dimensional equa- 
tions of state of the van der Waals type, with 
wo corresponding constants a and 6, have been 
ntroduced for the study of adsorbed films. ‘“‘The 
usefulness of all such equations is limited to 
moderate pressures and, as in three dimensions, 
neither @ nor 6 are really constant.”’ ? 

lt was pointed out by Fowler*® that the dif- 
culties in going beyond the existing treatments 
ire of exactly the same nature as those en- 
countered in the establishment of equations of 
state for compressed gases and liquids. This sug- 
yests that any progress made in the method of 
evaluating the partition sum in the three-dimen- 
sional case will, almost automatically, result in 
applicability to the two-dimensional case of 
monolayers. 

Now the author has developed recently a 
statistical theory of liquids which derives the 


‘N. K. Adam, The Physics and Chemistry of Surfaces 
Clarendon Press, Oxford, 1938). Cf. pp. 40 and 116. 


*N. K. Adam, reference 1, p. 120. 

'R.H. Fowler, Statistical Mechanics (University Press, 
Cambridge, 1936), second edition, cf. p. 828. 

‘G. Jaffé, Phys. Rev. 62, 463 (1942) ; 63, 313 (1943). We 
eler to these two papers in the text as I, II, respectively. 


essential features of a liquid from the knowledge 
of the mutual potential energy of the molecules, 
and which, after the choice of a special form of 
potential, leads to a reasonable agreement with 
experimental data for liquids of widely different 
character. It seemed worth while to the author 
to attempt an extension of his method to the 
equation of state of monolayers. This is being 
done in what follows. Since the calculations 
involved are so very similar to those of the 
previous problem, the reader is referred for details 
to the two papers quoted. 


Il. THEORY 


Here again we are going to assume that the 
molecules are extensionless points, of mass m, 
and that the total potential energy of the mole- 
cules in the monolayer is the sum of binary con- 
tributions. These assumptions appear to be even 
more ebjectionable in the present case than in 
the idea'ization of liquids. The molecules forming 
monolayers consist mostly of long chains, and 
it is known that these chains lie flat on the surface 
of the carrier liquid in gaseous films whereas they 
stand up eventually under some angle, and are 
closely packed in condensed films. Consequently, 
there must be a domain of transition which can- 
not possibly be described by the simple model 
of central forces. 

A more complete treatment would have to 
take into account at least the dependence of the 
potential upon the orientations of two interacting 
molecules. This can be done and has been at- 
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tempted by Mitchell,5 who used a different 
approach to the problem. However, so as not to 
make the calculations too involved, it seemed 
desirable to the author to start from the simplest 
assumptions and to investigate how far they 
would lead. Naturally, the results can be ex- 
pected to hold, with a fair approximation, only 
for films which are either gaseous or liquid. In 
the domain of transition, “expanded films,” new 
phenomena must be anticipated. 

Furthermore, we are going to restrict the 
treatment to one-phase systems. The phenomena 
of ‘‘surface vapor pressure’ and of “‘boundary 
line tension’’ could be included in exactly the 
same way as vapor pressure and surface tension 
are included in the former theory. However, too 
few empirical data are known to warrant such a 
generalization. 

Let there be N molecules in a total surface of 
area S so that n= N/S represents the melecular 
surface density, and A = S/N the area disposable 
for each molecule. Let, furthermore, ¢(r;;) be 
the mutual potential energy of two molecules 
separated by the distance r;;. Then the total 
potential energy can be written in the form 
(compare I, Section 1) 


N 
(2.01) 


ful 


N 
U=> ui, 
i=1 


The gist of our method consists in transforming 
the sum in the expression for u; into an integral 
which is to be extended over the entire surface 
S with the exception of a little circle surrounding 
the molecule i. The radius of this circle ¢ is the 
characteristic parameter of our theory and will 
be determined separately. The surface S may be 
replaced by a circle of very large radius, provided 
boundary effects are to be neglected. Thus we 
obtain, dropping the index 1, 


u=nrx(c), (2.02) 
where x(a) is defined by 
x(o)= f g(r)rdr. (2.03) 


Next, the partition function has to be deter- 
mined (1, Section III). The normalized phase 


5]. S. Mitchell, Trans. Faraday Soc. 31, 980 (1935). 
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integral can be written as 
Q=0%Q,/N!}, 
= 2xmkT/h’, 


and 
a= ff: exp (— U/RT)dx,dy,- -dxydyy. 
(2.06) 


Since the potential energy, under our assump. 
tions, separates into N terms not depending on 
any coordinates, Q, becomes the product of 
identical factors: 


(2.04) 
(2.05) 


with 


(2.07) 
where q, is defined as 


exp (—u/kT). 


(2.08) 


Now the work function A becomes (with the 
usual approximation for N !) 


A=—kTInQ 
= —kTN{In q,—In N+1}, 


—A=kTN {In (OS/N)+1}—U. 


(2.09) 


or 


(2.10) 


Here we have reintroduced the total energy (I, 
Section V1) 


U=Nu=r(N/S) (2.11) 


Calling F the surface pressure, we obtain from 
(2.10) the equation of state of the monolayer in 
the form 


0A NkT 8U 


as aS’ 
which corresponds exactly to the three-dimen- 
sional case [I, form. (6.04) ]. 

The dependence on n and T of o, the average 
value for the distance between nearest neighbors, 
still remains to be determined. By calculations 
which are entirely analogous to the _ three- 
dimensional case (1, Section V and Appendix !]), 
we find the following results. Let the probable 
number of mass points in a small fractional area 
s at the distance r from a given mass point be 
pi(r)ns; then the probability p2(r)dr that the 
nearest neighbor to the given mass point is 0 
be found at distances between r and r+dr will 
become 


p2(r)dr=nP’(r) exp (—nP(r))dr. 


(2.12) 


(2.13) 
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Here P(r) is given by 


(2.14) 


P(r)=2r f p:(p)pdp. 


For Pi(r) we assume the form (1, p. 470) 
pi(r) =exp (— ¢(r)/kT), 


and replace the average value 7, 


(2.15) 


f= f (2.16) 


), the most probable value, i.e., the value which 
makes the integrand in (2.16) a maximum. 


Thus we obtain for the determination of our 
characteristic parameter o the equation 


exp (— =0, (2.17) 


which corresponds to our former equation I, 
5.05) in the three-dimensional case. 

In the limiting case ¢/kT-—>0, i.e., either for 
sufficiently small potentials of interaction, or for 
sufficiently high temperatures, (2.17) reduces to 


o=(1/2xn)! (2.18) 


in agreement with Hertz’s geometrical formula 
|, p. 469) for the two-dimensional case. 

For the calculation of 8U/AS in (2.12) the 
knowledge of d0/0S is required. From (2.17) we 


derive 


do 


with 


D=2—2ag’/kT +079" 


(2.19) 


(2.20) 


ln order to obtain numerical results we intro- 
duce, as in the case of liquids (I, Section VII), 
the following form of potential 


¢(r) = fr-*—gr-*. (2.21) 
Then the total energy becomes 
U=n(N?/S) 
X | O-*) /(X— 2) — /(u—2)}, (2.22) 
and the surface pressure 
F= NkT/S+a(N/S)*{ fo-°-?[1/(A—2) 


(2.23) 
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The relation between o and the area per 
molecule A is given, from (2.17), by 


(2.24) 
(2.24a) 


A =2ro*/5, 
exp (¢/kT). 


with 


Making use of this relation we can bring (2.23) 
into the form 


FA/kT =1+ (6/2) 
(2.25) 


It should be remembered that in this equation 6 
and « are still functions of « and JT given by 
(2.24a) and (2.19). Since the relation between ¢ 
and A cannot be solved with regard to o in an 
explicit way, it is necessary, for numerical cal- 
culations, to start from given values of o and |. 
to determine A and F as functions of ¢. Thus 
(2.24) to (2.25) represent the equation of state 
in the parameter form. 

A considerable simplification of our result can 
be realized if we define o by the purely geo- 
metrical formula (2.18). This amounts to taking 
the limit ¢/k7-—+0, which makes 6=1 and e=}. 
Then the equation of state takes the form 


FA/kT =1+(1/4) { f/RT)(A/(A—2))o™ 
— (g/RT)(u/(u—2))o-*}, 
o=(A/2r)!. 


(2.26) 
(2.27) 


with 


We shall require the expression for the “bulk 
coefficient B,”’ i.e., the reciprocal of the com- 
pressibility b= —S(@F/dS). It can be derived 
from (2.23), if (2.19) and (2.20) are taken into 
account. For the condensed state a suitable 
approximation can be obtained by treating 6 (for 
isothermal changes) as a constant (I, Section 
VII, A). Thus we obtain 


—g[u?/(u—2) Jo-*}. 


In the same approximation (which amounts 
to setting «= 4) we obtain for the equation of state 


(2.28) 


F=kT/A+(6/4A) 


X 1 f(A/(A—2))o* (2.29) 
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Fic. 1. FA/kT is represented as function of F (in dynes/cm) for insoluble films of the diethyl esters of the two 
acids COOH(CH:),COOH with »=10 and n=11. The curves are theoretical, the points observations by Adam and 


Jessop. 


Finally, the heat of (surface) vaporization per 
molecule is given by 
—u=—(8/2){f[1/(A—2) 
—g[1/(u—2)]o~*}. 


This follows from (2.22) and (2.24). 


(2.30) 


Ill COMPARISON WITH EXPERIMENTAL DATA 
A. Gaseous Films 


From what has been said above it follows that 
gaseous and condensed films must be treated 
separately. For gaseous films we are going to use 
the formulae (2.26), (2.27). It is true that the 
condition ¢/kT<1 is not realized except for the 
highest dilutions. We have, therefore, carried 

6 The numerical calculations based on the exact formulae 
show that, as ia the case of liquids, ¢/kT comes out of the 


order 1, and og’/kT and o*y"’/kT correspondingly larger 
(I, p. 472). 


through a considerable number of calculations 


based on the exact formulae (2.24) to (2.25 
These calculations have convinced us that, in th 
cases tested, the two sets of equations may bx 
made to represent the data about equally well 
though with different values of the exponents \ 
and 

It should be pointed out that the two sets ol 
formulae are different only in the calculation o 
the average distance of nearest neighbors. Th 
simpler geometrical expression (2.18) disregards 
the influence of the force field about a molecu 
on that average distance. On the other hand, the 
expression (2.17) probably exaggerates the influ- 
ence of this force field as was pointed out 0! 
several occasions (I, p. 470, p. 475; II, p. 316 
We shall see that the study of condensed films 
leads to the same conclusion, and it seems that 
the geometrical formula represents, for is* 
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thermal changes, at least as good an approxima- 
‘ion as the dynamical formula. 

Under these circumstances we have thought 
that it was not worth while to carry through the 
cumbersome calculations based on (2.24) to 
?.25), and we have used the simpler expressions 
2.26) and (2.27). With these formulae, com- 
paratively low values for the exponents A and u 
have to be chosen. However, the choice is by no 
means very definite if observations are available 
only on the equation of state. We have adopted 
the values A=5 and A=6, and w=3. With the 
eact formulae larger values have to be chosen, 
namely, A~10, and p=6. It follows that not too 
much significance can be attached to the specific 
values of X and yw. However, that is evident from 
the fact that the simple notion of central forces 
can, at best, represent a rough approximation to 
the interaction of molecules in the monolayer. 

\Ve have carried through the calculations for 
wo sets of empirical data on gaseous films. The 
first set deals with insoluble films of dibasic 
esters having a carboxyl group at each end. The 
measurements are due to Adam and Jessop’ and 
iyproach, at low pressures, the theoretical value 
lor gaseous films with an accuracy of about 10 
percent. The experimental points for the diethyl 
eters of the two acids COOH(CH;),COOH, 
with n»n=10 and n=11, are represented in Fig. 1 
where FA/kT is plotted as function of F.* The 
curves are theoretical curves calculated with the 
values of the constants AX, yu, f, and g given in the 
irst two rows of Table I. Though there seem to 
be systematic deviations at low pressures, the 
general trend of the observations and the almost 
linear increase of FA /kT with higher values of F 
is well rendered.® 

The second set of observations concerns ad- 
vorbed films of fatty acids, from 4 to 12 carbons 
long, on dilute solutions. In this case the surface 
pressures have not been measured directly but 


*\. K. Adam and G. Jessop, Proc. Roy. Soc. London 
All2, 376 (1926). 

*We have taken the observed data from N. K. Adam, 
relerence 1, Fig. 12. 

"Not much importance should be attached to the fact 
‘nat for m= 10, the highest point is considerably off the 
‘neoretical curve. The authors state (reference 7):“*. . . the 
‘act that the curve for 10 appears below that for 11 is 
probably due to the former being so volatile that the 
molecules leave the film too rapidly for exact measure- 
ments,”” 


have been calculated by Schofield and Rideal'® 
from earlier measurements of the surface tension 
of solutions," by the aid of Gibbs’ adsorption 
equation. However, it has been pointed out by 
Adam" that the series of curves thus deduced 
shows a continuous transition without break to 
the curves obtained by direct measurements for 
longer chains and that in the case of the 12 
carbon acid the two methods yield concordant 
results. Therefore, this combined set of observa- 
tions represents a particularly valuable series 
since it leads from films which are distinctly 
gaseous through conditions in the neighborhood 
of the “critical point’’ to films which are dis- 
tinctly condensed. 

In Fig. 2 the drawn-out curves are calculated 
from Eqs. (2.26) and (2.27) with the constants 


TABLE I. Constants used in the calculation of Figs. 1 and 2. 


t=8°C 
Substance X10" 


Diester Cio 
Diester Ci, 


Substance 


Butyric acid 
Valeric acid 
Caproic acid 
Caprylic acid 
Capric acid 
Lauric acid 


tabulated in Table I in lines 3 to 8. The data of 
Schofield and Rideal are inserted, and it will be 
seen that the agreement is about as good as 
might be anticipated from the nature of the data, 
Schofield and Rideal have pointed out the close 
resemblance of their curves with the well-known 
isothermals of CO, above and below the critical 
temperature. It is this very point which is 
brought out clearly by the theoretical curves. 
The constants in Table I have been found by 
trial and error, without any reference to the 
chemical constitution of the substances studied. 
It will be noticed that both f and g are almost 


”R. K. Schofield and E. K. Rideal, Proc. Roy. Soc. 
London A109, 57 (1925); A110, 167 (1926). 

" B. von Szyszkowski, Zeits. f. physik. Chemie 64, 385 
(1908); A. Frumkin, Zeits. f. physik. Chemie 116, 466 
(1925). 

K. Adam, reference 1, p. 118. 
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Fic. 2. FA/kT is represented as function of F (in dynes/cm) for adsorbed films of fatty acids with 4 to 12 carbons. 
The curves are theoretical, the points represent determinations by Schofield and Rideal. 


linear functions of the number m of carbons in 
the chain. 


B. Condensed Films 


From the last-mentioned fact it is possible to 
extrapolate the values of f and g to higher values 
of n. If that is done the resulting isothermals, 
i.e., F plotted as function of A, begin to descend 
under the A axis. For »=13 the isothermal 
becomes almost exactly tangent to the A axis 
(for Ayo=46 sq.A). For higher values of » there 
are two intersections. Of these only the inter- 
section on the descending branch can represent 
a physically possible state since the other inter- 
section would involve a negative compressibility. 

The physically possible intersections indicate 
states in which the condensed film could exist 
under a surface pressure zero. This fact should 
not be considered a rigorous consequence of our 
theory. The consideration of boundary phe- 
nomena would introduce a finite surface vapor 
pressure, exactly as in our theory of liquids the 
tension necessitates a finite vapor 


surface 


pressure (I, Section IV). However, the inter- 
sections in question represent the “‘zero approxi- 
mation” to the condensed state. 

Thus it follows, even from the extrapolated 
values, that above n=12 the films become con- 
densed, as they really do. The values of Ao thus 
obtained by extrapolation (A»= 36, 33, 27 sq.A 
for n= 14, 15, 18, respectively) are of the proper 
order of magnitude but change with m, whereas 
it is characteristic of the condensed films in this 
(and other) series that A» is constant, having 
the value A»= 20.5 sq.A for the fatty acids. 

This fact, of course, cannot be brought out by 
potentials obtained from the study of gaseous 
films. As pointed out in the beginning of Section 
II the potential must change in form when the 
films pass from the gaseous to the condensed 
state, and extrapolated values can, at best, give 
qualitative results. Therefore the liquid films 
must be studied separately in the present form 
of our theory. 

If sufficient data were available, the four con- 
stants involved in our potential (2.21) could be 
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determined from our equations (2.28), (2.29), and 
(2.30) to the same extent that this was done for 
the liquid state previously (I, Section VIII). 
Unfortunately the available data are very 
meager. 

For the fatty acids from n=13 to n=16 the 
surface vapor pressure has been determined 
directly by Adam and Jessop."* This is equivalent 
10a knowledge of the heat of vaporization. How- 
ever, the compressibilities are known only 
roughly and cannot be used for the determination 
of \ and yw. Under these circumstances the ex- 
ponents A and yw have to be assumed from other 
considerations. 

lf the surface vapor pressure Fo is known for 
any temperature and if it is assumed that the 
vapor film obeys the gas law 


Fo=kT/Av’, (3.01) 


A, ie., the area per molecule in the surface 
vapor, can be calculated. Then Boltzmann's 
theorem yields (1, p. 467) 


Ao/Ao' =exp (—I/kT), (3.02) 


where /=—w is the heat of vaporization per 
molecule. Hence / can be calculated from the 
knowledge of Ag and Fo. 

Now from (2.29) (with A =A» and F= Fo) and 
from (2.30) the two combinations of constants" 


x=bfoo, y=dgou* (3.03) 


can be calculated provided \ and yu are assumed. 
The value ¢9 corresponds to the condensed state. 

Thus our potential, even with \ and y given 
within certain limits, is flexible enough to give 
the correct values for Ay and Fo. An independent 
test can then be obtained by calculating the bulk 
cocthcient B from (2.28) and comparing it with 
experimental data. 

We have carried through these calculations for 
ilms of tridecylic, myristic, pentadecylic, and 
palmitic acids, and for two sets of values for \ 
and «. The first set is given by A=6, u=3, similar 
‘o that used for gaseous films above. The second 
set, \=9 and 4 =6, was chosen because it repre- 
sents the empirical values of density, compres- 
sibility, and heat of vaporization for liquid 


"\. K, Adam and G. Jessop, Proc. Roy. Soc. London 
All0, 423 (1926). 
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decane without error. Decane data were used 
because it is the hydrocarbon with the highest 
n for which the three data in question are 
available. 

The results of our calculations are tabulated 
in Table II, together with the corresponding 
constants for decane (as a liquid). It will be seen 
that the constants (3.03) can be determined, for 
both sets of exponents, in such a way that Ao 
and Fy come out without error; however, the 
bulk coefficients are different in the ratio of about 
1:3.4. 

From indications of Adam" it can be derived 
that the bulk coefficient for condensed films of 
fatty acids is ‘‘usually about’’ B = 1250 dynes/cm. 
‘In some instances, on chains of 18 carbon atoms 
and upward in length”’ values as high as B = 2500 
dynes/cm are observed. It would follow from 
this that our theoretical values with A=9 and 
u=6 are just about correct. If necessary, they 
could be increased by choosing larger values of 
and 

Furthermore it will be seen that the values of 
the ‘energy terms’’ (3.03) are reasonable inas- 
much as they are very nearly the same as those 
obtained for liquid decane in bulk by a similar 
calculation. 

The calculation of the bulk coefficient can be 
made without an explicit knowledge of 6 and apo. 
However, the value of 6 is involved implicitly 
[by its definition (2.24a)] and can be deter- 
mined in a way similar to that for liquids (1, p. 
475).'5 When 6 is known, oo can be determined 
from (2.24). 

The values of 6 and oo thus determined are 
added in Table I]. It will be seen that they dre 
well in line with the corresponding values ob- 
tained for liquid decane. However, the values of 
ao appear to be too large, even more so for the 
monolayers than for the liquid. As stated pre- 
viously, this should be due to the fact that the 
probability function p,(r) [see (2.15)] exag- 
gerates the influence of the force field about a 
molecule. 

“™ N. K. Adam, Proc. Roy. Soc. London A101, 452 (1922). 


“ As there is a misprint in the formula I, (8.17) we give 
here the corrected form 


6=(1—c2/(25kT)) exp [c,/(5kT) (8.17) 


It should be pointed out that the definition of 3 is slightly 
different in the two-dimensional case from the a. 
dimensional one. 
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TaBLE II. Constants referring to films of tridecylic, 
myristic, pentadecylic, palmitic acid, and to liquid 
decane. 


Cie Decane 


Substance Cis Cis 
14.5°C 23°C 


Cis 
Temperat. 14.5°C 14.5°C 


foo 108 
1018 
B dynes/cm 
6 


oo X 108 

108 
bgoo ®X 108 
B dynes/cm 


oo X 108 
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If the geometrical formula for @ is accepted, 
this amounts to setting 6=1 in Eqs. (2.28) to 
(2.30). The theoretical values of the bulk coef. 
ficient remain the same, but the mean distances 
between nearest neighbors become, uniformly 
for all films and with both sets of exponents, 
o9=1.8A. This appears to be a little too small, 
and a more correct value of the function p,(r) 
would probably lead to a value intermediate 
between the two values given here. 

On the whole it can be said that the agreement 
between empirical data and our theory is as good 
as might have been anticipated, given the pre- 
liminary nature of the potential of interaction. 
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It is pointed out in the first place: (1) in Birkhoff's 
gravitational theory based on ‘“‘flat’’ space-time, the “‘red 
shift’’ is accounted for by the energy change of the photon 
as it travels from the emitting body, whereas the photon 
plays no especial role in the Einstein theory; (2) the 
solution of the problem of two or more bodies is feasible 
in the new theory because of its simpler character. Four 
comments of H. Weyl concerning the Birkhoff theory are 
discussed, and it is concluded that these are to be taken 
with much reserve. In regard to the third of these com- 
ments it is pointed out that the “‘perfect fluid’ used by 
Birkhoff as‘the ultimate carrier of mass and electric charge 
is to be characterized as the simplest fluid with disturbance 
velocity that of light (c). It is affirmed to be a glaring 
defect of earlier relativistic theories that the disturbance 


I. PRELIMINARY OBSERVATIONS 


N a recent number of a journal of wide cir- 

culation,! Hermann Wey! has given expression 
to several critical remarks on G. D. Birkhoff's 
theory of gravitation of 1942.* In this note we 
intend first to analyze briefly the substance of 
Weyl’s comments, secondly to consider the 
structure of the new theory from the physical 


1H. Weyl, Math. Rev. 4, 285 (1943). 
2G. D. Birkhoff, Proc. Nat. Acad. Sci. 29, 231 (1943). 


velocity in matter has been taken as arbitrary, although 
that of gravitation and of the electromagnetic field have 
been equal to c. The differential equations of the theory are 
then set up. An additional cosmological term in the gravita- 
tional potentials A,; is suggested, namely, 


hy = 22) gay, 


where x, y, 2, fare Lorentz coordinates and K is the (small 
cosmological constant. The explicit formula for the rate 
of advance of periastron P of two bodies (mass points) of 
masses m and mz is given, as obtained from the solution of 
the two-body problem in the theory, and its possible appli- 
cation to double stars is referred to. The authors propose 
to give a detailed development of the theory and its appli- 
cations in papers to be published shortly elsewhere. 


point of view, and lastly to refer briefly to its 
physical applications. 

Before doing this, however, we would like to 
make certain general observations. 

The explanation of the “red shift” is funda- 
mentally different in the gravitational theories 
of Einstein and Birkhoff. In the new theory the 
red shift has turned out to be accounted for by 
the energy change of the photon as it travels 
from the emitting body to the Earth; this 
explanation fully takes account of the role 0 
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A =6, p=3 
11.2 12.6 14.2 17.4 
6.12 6.82 7.67 9.20 
279 321 355 447 
17.2 18.6 20.0 23.1 
7.5 7.8 8.1 8.7 
A=9, u=6 
42.9 47.7 53.8 64.3 64.1 
37.7 42.0 47.0 56.3 48.7 
958 1067 1206 1450 
6 26.5 28.1 30.8 34.1 34.3 
| 9.3 9.6 10.0 10.6 12.1 
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the light signal (photon). In Ejinstein’s inter- 
pretation, on the other hand, the frequency of 
the light wave emitted by an atom on a distant 
body is compared with the frequency of the 
light wave emitted by the same atom on the 
Earth, and any phenomena taking place while 
the light signal (photon) travels from the emitting 
body to the Earth play no part whatsoever. 

Again, the difficulties accompanying the two- 
body problem and other questions in the general 
theory of relativity are too well known to require 
mention, while in Birkhoff’s theory the solution 
of the problem of two or more bodies is quite 
within reach and will soon be available. This is 
due to the essentially simpler character of the 
latter theory. 

An objection which may be made to Birk- 
hoft's theory is the introduction of an absolute 
reference system, which runs counter to Ein- 
stein’s general relativity principle that matter 
determines space, and analogous philosophical 
ideas previously developed by Ernst Mach. But, 
ina sense, this objection to Birkhoff might also 
be urged against Einstein because in the latter's 
theory a single rotating body can still be sup- 


posed alone in the universe, which is absurd 
from Mach’s point of view. It would be difficult 
indeed to set up any physical theory to which 
objections of this general nature could not be 


raised. 


Il. ANALYSIS OF WEYL’S COMMENTS 


To begin with, Weyl contends that Birkhoff's 
theory is much the same as Einstein's theory of 
1916, for the case of weak gravitational fields. 
ln a note appearing elsewhere, Barajas* has 
shown that the factual consequences of Einstein's 
theory in this case differ from Birkhoff's ; further, 
that the choice of the gravitational potentials 
suggested by Wey] is unsatisfactory ; and, lastly, 
that the trajectories of a test particle in Birk- 
hoff's theory are not geodesics ‘n any four- 
(imensional space-time. 

Weyl further raises four fundamental objec- 
tions to Birkhoff’s theory: (1) ‘“‘The connection 
between metric and gravitation is dissolved."’ 
2) The proportionality between inertial and 
gravitational mass “has again become as mys- 


*\. Barajas, Proc, Nat. Acad. Sci. 30, 54 (1944). 
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terious as it was before Einstein.’’ (3) Birkhoff's 
perfect fluid “‘appears as a primitive irreducible 
physical entity.”’ (4) ‘There seems to be no 
indication that the mechanical equations spring 
from a universal law of conservation of energy 
and momentum.” We now propose to take up 
these points in detail. 

(1) With regard to the first, it should be 
recalled that the connection between matter and 
geometry, as developed by Einstein, is purchased 
at the expense of giving up a fundamental 
reference system. This implies abandoning the 
description of nature in terms of four funda- 
mental independent variables, essentially unique 
except for the arbitrariness involved in position 
and velocity of a single point (Lorentz group). 
In the abandonment of such coordinates may be 
discerned a sufficient reason for the early ex- 
haustion of all observational tests of the general 
theory of relativity, and also for the fundamental 
difficulty of assigning actual physical meaning to 
the coordinates introduced in problems of this 
theory. Indeed, almost thirty years of intensive 
research have failed to provide another test 
besides the three classical ones, or to apply the 
theory to other fields of physics. In spite of a 
great deal of work, there does not exist any satis- 
factory solution of the two-body problem, and 
the n-body problem has not even been formu- 
lated. Even in the simple case of Schwarzschild’s 
solution of the one-body problem, no clear-cut 
physical interpretation of the Schwarzschild 
coordinates seems to be available. This last dif- 
ficulty is so serious that it led Milne to give up 
the use of general coordinates in relativity. 

(2) Weyl’s second objection seems to be based 
largely on a misinterpretation. In Birkhoff's 
theory, just as well as in Einstein’s, the equality 
between inertial and gravitational mass comes 
out of the fact that in the gravitational equations 
the energy tensor is linear in the mass density. 
While nothing similar to Einstein's ‘‘equivalence 
principle” has been explicitly formulated in 
Birkhoff's theory, it must be stated that the 
exact significance of the principle has yet to be 
found. As far as can be seen at the present time 
it merely asserts that bodies moving freely in 
empty space near attracting matter behave, 
relatively to a hypothetical attached reference 
system, in the same way under all circumstances 
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—insofar as they so behave! From the point of 
view of differential geometry, it amounts to 
saying that in the infinitesimal neighborhood of 
any point in a Riemannian curved space there 
exists always a tangent flat space, and this is 
just as true in Einstein's theory as in any other. 

According to Einstein, the proper language for 
the description of nature is that of tensors, to- 
gether with the underlying group of general 
transformations; and the equivalence principle 
expresses the theorem just stated. According to 
Birkhoff, the proper language for this description 
is that of four vectors, and the underlying space- 
time is everywhere flat. In the latter case, the 
basic group involves only ten arbitrary con- 
stants, while in the former it involves four com- 
pletely arbitrary functions of four completely 
arbitrary variables. Whichever choice is made, 
the equality between gravitational and inertial 
mass follows as a consequence both of Birkhoff’s 
and Einstein's theories. 

(3) Birkhoff has chosen the name “perfect 
fluid” for perhaps the simplest substance in which 
all disturbances are propagated with the velocity 
of light, on account of the similarity between the 
equations governing such a substance and those 
of a perfect fluid. The particular type of perfect 
fluid considered by Birkhoff is characterized by 
only one scalar (the mass density) and vector 
(the velocity). Fundamental difficulties arise if 
the disturbance velocity is different from that 
of light. From our point of view, the fact that in 
Einstein’s theory no attention has been paid to 
this requirement constitutes a glaring defect, 
which by itself explains the possibility of Birk- 
hoff’s theory. Indeed, the postulation of a pri- 
mordial substance in which all disturbances are 
propagated with the velocity of light is funda- 
mentally a consequence of the assumption that 
the basic space-time is everywhere that of 
Minkowski. This is actually the only physical 
assumption made in Birkhoff’s theory. 

(4) Weyl’s fourth objection also seems to 
spring from a misunderstanding. In Birkhoff's 
theory the law of conservation of energy and 
momentum at material velocities small compared 
with that of light is as much a consequence of the 
mechanical equations of motion (and conversely) 
as in Einstein's. It is true that in Birkhoff’s case 
the conservation of energy and momentum is not 


connected with a fundamental geometrical 
theorem (Ricci’s theorem), as in Einstein's 
theory. But, although Einstein obtains the 
formal result that the divergence of the energy- 
momentum tensor must vanish, this does not 
imply the conservation of energy and momentum 
in an exact sense because the four-dimensional 
integrals of a covariant partial derivative jn 
curved space-time cannot be transformed into 
three-dimensional integrals. Consequently the 
conservation theorems of Birkhoff's theory are 
at least much more precise. 

These remarks may serve to point out that 
Weyl’s assertions are to be taken with much 
reserve and that additional research is required 
before the usefulness of Birkhoff’s theory for 
physics can be adequately assessed. 


Ill. STRUCTURE OF THE BIRKHOFF THEORY 


In the light of the preceding remarks, it is 
possible to recapitulate as follows the considera- 
tions which lead to Birkhoff's theory. The funda- 
mental concept of electromagnetic space-time, 
associated with the names of Fitzgerald, Larmor, 
Lorentz, Einstein, and Minkowski, has been of 
the first importance for physics and has played 
an ever increasing role in the developments of the 
last fifty years. The revolutionary change which 
this concept has brought about in physical theory 
is reflected in the mathematical apparatus em- 
ployed. Four homogeneous variables of space and 
time replace the three homogeneous variables of 
space and the single disparate variable of time; 
the underlying Lorentz group replaces the 
Galilean group, and the language of 4 vectors 
replaces the language of 3 vectors characteristic 
of classical physics. 

The early attempt of Nordstrém (1912) and 
others to incorporate gravitation in this frame- 
work failed to explain certain delicate gravita- 
tional phenomena which alone provide a crucial 
test. Thus, without an intensive study, electro- 
magnetic space-time was abandoned for a curved 
or Riemannian space-time, latent in the ideas of 
Minkowski and realized in Einstein's brilliant 
gravitational theory of 1916. 

It is exceedingly easy to exaggerate the sig- 
nificance of the three so-called critical confirma 
tions of this theory, of which by far the most 
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certain is afforded by the excessive perihelial 
advance of the planet Mercury. In fact, the 
building of theories from the aesthetical-mathe- 
matical point of view has shown that dimensional 
considerations enter which lead always to the 
sume result, aside from simple numerical factors. 
More definitely, the three formulas for perihelial 
advance, bending of light, and red shift take 
always the respective forms 


qm/a(1—e*), rm/p, sm(1/r1—1/re) 


where q, 7, S, are simple numerical constants. 
Consequently, the basic test of any such theory 
really reduces to the single requirement that the 
first constant r has a value not very different from 
Einstein’s value 64! Thus the theory of Einstein, 
stripped of all mystical trappings, is seen in its 
proper perspective. It becomes obvious that the 
question ‘‘What is the simplest theory of gravita- 
tion and other physical phenomena, based on 


Table I. Table of disturbance velocities. 


Nordstrém, 
Einstein 


Newton, 


Maxwell Birkhoff 


Matter 
Gravitation 
Electromagnetic field 


Arbitrary Arbitrary 
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clectromagnetic space-time, which explains the 
known facts?” deserves the most careful con- 
sideration. This question becomes all the more 
urgent since the Einstein theory, with its enor- 
mous mathematical complication and its lack of 
proper independent variables, seems to be essen- 
tially unworkable. 

Birkhoff’s theory of 1942 provided perhaps the 
first thoroughgoing attempt to answer this ques- 
tion. Its point of departure arises from the valid 
criticism of earlier theories that the forms of 
matter therein employed are inconsistent with 
the actual requirement that the disturbance 
velocity be that of light. This situation may be 
roughly set forth in the comparative Table I, 
in which the entries indicate disturbance velocity. 
The physical necessity for a disturbance velocity 
c of matter appears from the fact that, with any 
other velocity, the equations of motion break 
down at the collision of two atomic constituents. 
Here the matter referred to is not gross matter, 
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but the ultimate carrier of mass and electric 
charge. 

If this requirement is admitted, it appears that 
the first ein of any new theory based on 
electromagnetic space-time should be the con- 
dition that matter has a disturbance velocity 
equal to that of light under all circumstances. 
The theory of Birkhoff is based on one such type, 
the “perfect fluid,’’ which seems to be the sim- 
plest conceivable from a conceptual point of view. 
But it appears almost certain that any other 
type of medium obeying this same fundamental 
requirement will lead to essentially the same 
gravitational theory. The “perfect fluid’’ may be 
approached in the following manner: The state 
of matter is regarded as characterized by a single 
scalar density p and vector velocity u‘=dx‘/ds 
in the sense of local causation; the equations of 
motion are to be linear in the rates of change of 
these variables; free equilibrium is possible at a 
certain density po; the disturbance velocity is to 
be that of light c, which becomes 1 if the light- 
second is the unit of distance. 

We have established the mathematical theorem 
that the equations of motion of the perfect fluid 
may be written in the 4-vector form 


div T=0T**/dx*=0, T= p(u‘ui—}g‘/), 


where gii=1, —1, —1, —1 for i=1, 2, 3, 4, re- 
spectively, and g;;=0 for i#j, with appropriate 
choice of the scalar density p, unique up to a 
choice of a unit. Furthermore, the perfect fluid 
so obtained satisfies the conservation principle 
that the integral is conserved (dv, 
element of volume referred to a rest system). 

The equation of motion under arbitrary force 
densities f‘ may be written correctly 


OT */dx* = f'. 


There is then the essential further requirement 
to be imposed on the forces f‘ that if the particles 
of the fluid return to an initial state of position 
and velocities, the density p must also return. 
For this, it seems necessary and is certainly suf- 
ficient that the condition of orthogonality 


fau* =0 


be identically satisfied. Thus the force vectors 
are always required to be identically orthogonal 
to the velocity vectors. A primary force vector of 
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this type is evidently the acceleration vector 
itself. 

Now all force vectors which have been used in 
previous relativistic theories are rational and 
integral in the components of the velocity 
vector u‘. It is therefore natural to set: 


fi=A'+ 


where the coefficients are functions of position 
alone, But the electromagnetic force is known to 
be linear and homogeneous in the velocities and 
also in the first partial derivatives of the electro- 
magnetic potential g,, which itself is of dimension 
0 in the unit of length and time. Similarly, in the 
theory of Einstein, the gravitational forces are 
homogeneous and quadratic in the velocities, and 
linear and homogeneous in the first partial de- 
rivatives of the gravitational potentials g;; which 
again are of dimension 0. It is therefore natural 
to assume that B and C pertain to electromag- 
netic and gravitational forces, respectively. 
Birkhoff imposes the analogous requirements 
upon the vector A and thus arrives at the first 
(covariant) form for the forces‘ 


f oy re ( 8 
= u*+ 
Pax Ox® Ox® = ax 


where y is his ‘‘atomic potential,’’ constant along 
the world-lines of the fluid, ¢; is the electro- 
magnetic potential satisfying the Maxwell- 
Lorentz equations (¢, density of electricity) 


0 
—{ — ——— = —4ron, 
Axe 


and where /;; is his gravitational potential satis- 
fying 


which is the generalized form of Poisson's equa- 
tion of the theory. 

Hence the gravitational theory involved may 
be singled out as follows: If the equations of 
motion for a small total amount of matter 
(absence of gravitational force) are written 


OT ‘*/dx*= f'; 


4 Birkhoff has changed the form of A; from d¥/dx* to 
pdy/dx* because of the dimensionality requirement. 
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then for the general case, these equations are 


ax? 


where f** is the gravitational force and where the 
gravitational potentials 4;; are defined by 


8x7, ;. 


Such is the structure of the Birkhoff theory as 
hitherto formulated. However, to account for the 
phenomena of nebular recession on this basis, it 
is necessary to suppose that at some time in the 
past there was a nuclear distribution of matter, 
with a wide range of velocities. But from the 
physical point of view it appears to be more 
natural to suppose an initial nuclear distribution 
at low relative velocities. In our recent inves- 
tigations, we have been led to the following 
extension of the original theory involving the 
introduction of a cosmological constant. 

From the formal point of view the following 
more general type of Poisson equation needs 
to be especially examined: 


With proper choice of a unit of density, it is then 
possible to specialize further the above equation 
to the form 


847, ;+Kgi;, 


where K is the ‘“‘cosmological constant,”’ which 
has the dimensions of a density and is supposed 
to be very small. 

Now, with such an extended Poisson equation, 
it is no longer possible to demand that the 
gravitational potentials h;; approach 0 regularly 
at «, nor even that these functions are linearly 
infinite. However, the condition may be imposed 
that hj; are regularly infinite to at most the 
second order, with a boundary distribution at 
infinity which is spherically symmetric in 4 
spatial sense. Obviously, the most natural pos- 
sibility from the electromagnetic point of view 
is to take the cosmological term in the gravita- 
tional potentials to be 


2.5, 


where x, y, 2, ¢ are any Lorentz coordinates and 
x=y=z=1t=0 is the origin in space-time. 
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lt is our purpose to publish shortly detailed 
studies of the theory outlined above, its cos- 
mological extension, and its applications. 


IV. APPLICATIONS AND TESTS 
OF THE THEORY 


lt has not yet been possible for us to make an 
intensive and thoroughgoing study of the Birk- 
hoff theory in relation to known gravitational 
phenomena ; and still less has it been possible for 
us to consider other possible applications. 

However, the basic two-body problem has 
already been solved by us to the requisite order 
of approximation. Thus, for instance, the formula 
for rate of advance of periastron P of two bodies 
mass points) of masses m, and mz has been found 
to be (in absolute units) 


Pp 2x 


which in case of an infinitesimal mass (m,=0, 
m,=m) reduces to that made familiar by the 
Einstein theory, 6mm/a(1—e®). Here, of course, 
the application to the rotation of the line of the 
apsides for double stars is at once suggested. But 
such a comparison with observation must await 
the further analysis of the observations them- 
vlves, so that appropriate data may be available 
in cases Where the rotation of the line of the 
apsides due to tidal forces is relatively small. 
We here naturally hoped also that the new 
theory would account for some of the hitherto 
unexplained features of lunar motion. However, 
this does not seem to be a likely outcome 
although it may turn out that the solution of the 
three-body problem (Earth, Sun, and Moon) 
will yield an explanation of the desired type. 
Instead, at the present writing, we are inclined 
to think that these supposed features are 
attributable to the large mass ratio (1/83) of 
Moon to Earth, which makes the convergency 
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of the mathematical calculations involved slow 
and uncertain. 

It seems not unlikely that the astronomical 
consequences of this theory for the interpretation 
of the red shift from extragalactic sources will be 
of special interest. 

As yet no means appear to be available for a 
decisive experimental test between the gravita- 
tional theory of Birkhoff and that of Einstein. 
This theoretical uncertainty is likely to continue 
for some time, especially as it appears to be very 
difficult to carry the Einstein theory to specific 
conclusions other than those found by him at 
the outset. 

There remains for later consideration the study 
of atomic phenomena on the basis of the electro- 
magnetic and atomic potentials. It is to be hoped 
that, by proper assignment of the atomic poten- 
tial y to the constituents of the atom, the ex- 
planation of atomic phenomena may be advanced 
and possibly the fundamental Schroedinger wave 
equation may be obtained in a_ conceptual 
manner.’ Here the proton and electron consti- 
tuting the hydrogen atom, for instance, are con- 
ceived of as freely interpenetrable and super- 
imposed in case of equilibrium and as oscillating 
under disturbances. 

In conclusion we would like to point out that, 
for the physicist, all mathematics is funda- 
mentally a form of abstract model building, of 
more or less general aspects of nature; and that 
no experiments which the physicist may perform 
in his laboratory can advance very far without 
free access to a variety of abstract models which 
are not to be thought of as final. The new theory 
of matter, electricity, and gravitation in flat 
space-time proposed by Birkhoff, would seem to 
afford a model of unusually fundamental, simple, 
and complete type. 


5 See G. D. Birkhoff's two notes in the Proc. Nat. Acad. 
Sci. 13, 160, 165 (1927). 
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The density of potassium chloride has been determined 
by the method of “‘crystal suspension” in a liquid mixture 
of bromoform, n-hexanol, and n-pentanol. The density of 
the liquid mixture at the average suspension temperature 
of fifty-three different crystals was taken as the density 
of the potassium chloride. The value determined was: 


P28.098°C = 1.98651 +0.00002 g/ml 
which when corrected to 25°C gave: 
p2s°c = 1.98721 +0.00002 g/ml 
= 1,98715+0.00002 g/cm’. 
It was found that six precipitations of potassium chloride 
from aqueous solutions with hydrogen chloride were neces- 


sary to effect purification to constant density (+4 10~° 
g/ml). Exposure of crystals to the atmosphere was found 


to produce some manner of surface contamination which 
resulted in a change of density by as much as 7 X 10™ g/m| 
within two hours. Storage of crystals over phosphorus 
pentoxide under vacuum eliminated this effect. An an- 
nealing of crystals for four hours at 50° below the melting 
point with gradual cooling over a period of four hours was 
found necessary to obtain agreement of suspension tem- 
peratures among different crystals. By combination of 
density and x-ray data, atomic weights have been cal- 
culated. The values of the atomic weights of fluorine and 
calcium were found to be 18.9967+0.0010 and 40.0851 
+0.0011, respectively. In the crystals of calcite, lithium 
fluoride, and potassium chloride which were used to obtain 
the data underlying these calculations, there is no evidence 
for a secondary structure as proposed by Zwicky. 


I. INTRODUCTION 


PRECISION method for the determination 

of the densities of solid crystalline sub- 
stances has been developed by C. A. Hutchison 
and H. L. Johnston.! 

The method consists of suspending?’ crystals of 
the substance whose density is to be determined 
in a liquid of the proper density. The above 
workers applied this method to the determination 
of the density of very pure lithium fluoride. 
Later, H. L. Johnston and the writer*® were able, 
with certain refinements of the method, to obtain 
a value for the density of very pure sodium 
chloride with reduced probable error over that 
of the lithium fluoride determination. The present 
paper deals with the application of the ‘‘sus- 
pension method” to the determination of the 
absolute density of very pure potassium chloride. 

The purposes of the present determination are : 
(1) to obtain density data which, in conjunction 


*Faculty Associate, Engineering Experiment Station, 
and Instructor in Chemistry. 

1 (a) C. A. Hutchison and H. L. Johnston, J. Am. Chem. 
Soc. 62, 3165 (1940); (b) H. L. Johnston and C. A. 
Hutchison, J. Chem. Phys. 8, 869 (1940). 

?In previous papers the words “‘float’’ and ‘‘flotation” 
are used. These words do not clearly express what is done. 
Since the temperature is found at which a crystal would 
be suspended by a given suspension medium, the writer has 
used the words ‘‘suspend”’ and “‘suspension”’ in place of 
the above words. 

37H. L. Johnston and D. A. Hutchison, Phys. Rev. 62, 
32 (1942). 


with x-ray data, can be used in the calculation of 
various physical constants, e.g., atomic weights," 
and Avogadro's number ;° (2) to obtain density 
data which can be used in detecting changes in iso- 
topic abundance ratios of the elements potassium 
and chlorine due to a given isotopic separation 
process.! 6 

The most accurate values for the density of 
potassium chloride, reported to the present time 
in the literature, appear to be those of Kohler’ 
and of Tu. 

Kohler determined the density to be 1.990! 
+0.00015 g/cm’ at 15°C while Tu obtained the 
value 1.98930+0.00014 g/cm*® at 18°C. When 
correction is made for the coefficient of thermal 
expansion of potassium chloride,® Kohler’s value 
reduces to 


pos = 1.9879+0.00015 g/cm? 
and Tu’s value reduces to 


pos? = 1.98773+0.00014 g/cm’. 


4 (a) C. A. Hutchison and H. L. Johnston, J. Am. Chem 
Soc. 63, 1580 (1941); (6) C. A. Hutchison, J. Chem. Phys 
10, 489 (1942). 

®*R. T. Birge, Phys. Rev. 61, 206A (1942). 

*H. L. Johnston and D. A. Hutchison, J. Chem. Phys. 
10, 469 (1942). 

7]. Kohler, Zeits. f. Physik 78, 375 (1932). 

5 Y. Tu, Phys. Rev. 40, 662 (1932). f 

a(cubical) = 113 10-*. Klemm, Tilk, and Mullenhei™, 
Zeits. f. anorg. allgem. Chemie 176, 1 (1928). 
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KXohler’s value exceeds the result of the present 
determination by 0.038 percent while Tu’s value 
is greater by 0.029 percent. Neither of these 
workers attempted purification of his potassium 
chloride, but used the best commercial product 
available. Their larger values may be due in part 
to the presence of impurities and to a greater 
extent to the effect of prolonged exposure to the 
atmosphere, which, as found in the present work, 
generally increases the density value (see 
Table 11). 

The result of this investigation on very pure 
potassium chloride when corrected to 25°C gives 


pos’ = 1.98721+0.00002 g/ml 
= 1.98715+0.00002 g/cm?’ 


with reduced probable error when compared with 
other determinations. 


II. EXPERIMENTAL DESCRIPTION 


A. Preparation of Chemically Pure 
Potassium Chloride 


Chemically pure potassium chloride was pre- 
pared by a method similar to that used by H. L. 
Johnston and the writer® in their work on the 
density of sodium chloride. In essence this is the 
method employed by Richards and Wells" in 
their determinations of the atomic weights of 
sodium and chlorine. 

Pure hydrochloric acid was prepared by distil- 
lation of C. P. Baker’s Analyzed reagent hydro- 
chloric acid to which had been added a small 
quantity of C.P. potassium permanganate. The 
distillation was carried out in a still composed of 
a flask and condenser made from pure fused 
quartz. The middle third of the distillate was 
saved from each distillate of two successive dis- 
tillations. The pure hydrochloric acid was stored 
for future use in clean Pyrex bottles, the inner 
surface of which had been exposed to concen- 
trated C.P. hydrochloric acid for about four 
months prior to storing. 

Pure sulphuric acid was prepared in a similar 
manner by double distillation of C. P. Baker’s 
Analyzed reagent sulphuric acid. This acid was 
also stored in Pyrex which had been exposed to 


 T. W. Richards and R. C. Wells, “A revision of the 
atomic weights of sodium and chlorine,” Carnegie Institute 
of Washington, Publication No. 28 (1905). 
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C.P. sulphuric acid for about four months prior 
to storing. 

Pure sodium chloride which was used in the 
generation of pure hydrogen chloride was frac- 
tionally precipitated three times from saturated 
aqueous solutions of sodium chloride by the 
addition of hydrogen chloride. The beginning 
material was Baker’s Analyzed C.P. sodium 
chloride. The hydrogen chloride used in these 
precipitations was generated by the reaction of 
pure sulphuric acid (see above) with an aqueous 
solution of Baker’s Analyzed C.P. sodium 
chloride and pure hydrochloric acid (see above). 

Pure water which was used in all reagent 
preparations was obtained by double distillation 
from an alkaline permanganate water solution in 
a pure fused quartz still. Water used in the deter- 
mination of the volume of the glass bob (see 
below) was triply distilled and had a specific 
conductivity of about 10-7 ohm~. 

Pure hydrogen chloride which was used in the 
purification of potassium chloride was generated 
by the reaction of pure sulphuric acid (see above) 
with an aqueous solution of pure sodium chloride 
(see above) and pure hydrochloric acid (see 
above). The generated gas was passed through 
two tubes containing ground glass covered with 
pure sulphuric acid (see above) and finally 
through a tube containing Baker's Analyzed C.P. 
phosphorus pentoxide which had been sublimed 
into this tube. 

Purification of the potassium chloride used in 
the density determination was effected by suc- 
cessively precipitating potassium chloride from a 
saturated aqueous solution by the addition of 
pure hydrogen chloride (see above). After each 
precipitation the potassium chloride was centri- 
fuged to remove the mother liquor. The beginning 
material was Baker’s Analyzed C.P. potassium 
chloride. Precise density determinations showed 
that after four precipitations the densities of the 
potassium chloride from different samples did 
not agree within the experimental limits of the 
density determination (see Table III). However, 
constant density values were obtained from dif- 
ferent potassium chloride samples after six pre- 
cipitations. As a precautionary measure eight 
precipitations were performed. 

After the eight successive precipitations, the 
hydrochloric acid was expelled from the purified 
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salt by fusion in a platinum crucible. The salt on 
cooling was dissolved in pure, triple-distilled 
water (see above) and given three crystallizations. 


B. Preparation and Calibration of the 
Suspension Liquid 


For the crystal suspensions a liquid had to be 
found which had the proper density, the proper 
change of density with temperature, and which 
would not decompose or react with the potassium 
chloride crystals. A suitable liquid mixture, 
similar to that used by H. L. Johnston and C. A. 
Hutchison’ in their determination of the 
density of lithium fluoride, was found to have the 
following composition: bromoform, 40.00 ml; 
n-hexanol, 881 drops; m-pentanol, 832 drops (52 
drops equal 1 ml). The suspension liquid was 
prepared as follows: Baker’s bromoform (boiling 
point range, 150°-151°C) U.S.P. IX, was washed 
with water, dried with calcium chloride, and 
distilled under vacuum. The middle third of the 
distillate was collected in a receiver containing a 
small amount of n-hexanol and n-pentanol (25 
drops of each) which acted as a negative catalyst 
preventing decomposition. The composition of 
the bromoform was adjusted to that of the above 
suspension liquid with Eastman practical grade 
alcohols. The isothermal density of this suspen- 
sion liquid was found to remain constant within 
the experimental limits of error (+5 X10~-® g/ml) 
of the density determination for a period of about 
seven days. Suspension temperatures after the 
seven-day period were corrected by means of the 
standard crystals suspended within the seven- 
day period. 

The density of the suspension liquid at dif- 
ferent temperatures was determined by means 
of hydrostatic weighings. A Pyrex glass bob 
about 9.4 ml in volume which contained mercury 
for additional weight was weighed in air, in 
water, and in the suspension liquid at different 
temperatures. As in the sodium chloride work,’ 
to eliminate surface tension effects! at the point 
of entrance of the 0.01” platinum suspension wire 
into the liquid, an auxiliary mercury-weighted 
Pyrex bob about 1.1 ml in volume was attached 
to the lower end of the suspension wire. The large 
bob was detachable from the auxiliary bob. The 


"Osborne, McKelvey, and Bearce, Bull. U. S. Bur. 
Stand., Reprint No. 197 (1912). 


weighing procedure has been described _preyi- 
ously.’ Two series of weighings were made in the 
suspension liquid, those of the wire suspension 
with the auxiliary bob immersed and those of the 
wire suspension with the auxiliary bob and large 
bob immersed. 

TABLE I. Density of the bromoform mixture as a func- 
tion of temperature. The bob volume at 25.000°C was taken 
as 9.41336 ml and corrected to the temperature of each 


weighing by means of the cubical coefficient of thermal 
expansion of Pyrex. 


Density (g/ml) 
(at 28.098°C) 


1.98652 
1.98651 
1.98651 
1.98652 
1.98650 
1.98652 
1.98651 
1.98651 
1.98651 
1.98651 
1.98650 
1.98650 
1.98651 
1.98651 
1.98650 
1.98649 
1.98649 
1.98651 


Temperature Density (g/ml) 
(Centigrade) (at T) 


29.320 1.98422 
29.068 1.98468 
28.819 1.98515 
28.547 1.98567 
28.304 1.98611 
28.048 1.98661 
27.806 1.98706 
27.543 1.98756 
27.310 1.98799 
27.060 1.98846 
26.802 1.98894 
26.556 1.98940 
26.301 1.98989 
26.072 1.99032 
25.794 1.99084 
25.553 1.99128 
25.289 1.99178 
25.027 1.99229 


The precise volume of the bob was determined 
by weighings in triple-distilled water (see above) 
which had been heated to about 80°C before use 
to expel dissolved gases. Twenty sets of weighings 
were made in water at temperatures ranging from 
26.174°C to 27.769°C. These twenty inde- 
pendent volume determinations when corrected 
to 25.000°C by use of the cubical coefficient of 
thermal expansion of Pyrex glass" yielded an 
average bob volume of 9.41336+0.00002 ml. The 
densities of the water at the various temperatures 
were obtained from the International Critical 
Tables."* The average deviation from the mean 
of the twenty determinations was +0.00001 ml, 
and the extreme deviation was — 0.00009 ml. Al 
weighings were corrected for the buoyant effect 
of air. 

The temperatures at which the weighings were 
made were determined with a Beckmann ther- 


R. M. Buffington and W. M. Latimer, J. Am. Chem 
Soc. 48, 2305 (1926). The cubical expansion coefficient 
was taken as three times the linear coefficient given 4 
3.6X at 300°K. 

13 International Critical Tables (McGraw-Hill! Book Con 
pany, Inc., New York, 1928), Vol. III, p. 25. 
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mometer graduated in 0.01° intervals. The scale 
of this thermometer was shown to be quite 
uniform by the U. S. Bureau of Standards cali- 
bration of the scale. A standard mercury in glass 
thermometer with a temperature range of 18° to 
28°C and with graduations in 0.01° intervals, 
the scale of which had been calibrated by the 
U. S. Bureau of Standards, was compared with 
the Beckmann thermometer at nine different 
temperatures at 0.5° intervals over a temperature 
range of 1° to 5° Beckmann. The comparison at 
the extremes of the temperature range gave 


1.000° Beckmann = 29.320°C 
and 
5.000° Beckmann = 25.289°C. 


The densities of the bromoform suspension 
liquid are given in column two of Table I at 
temperatures covering the range 25.027°C to 
29.320°C. 


las_e II. Suspension temperatures of crystals grown from 
Baker's Analyzed C. P. potassium chloride.* 


Crystal Suspension temperatures in ° Beckmann 

desig Newly After After 

nation prepared 4 days 7 days See footnote® 

A2-1 1.998 2.010 2.007 

A2-2 1.991 1.999 2.005 

A2-3 2.009 2.004 2.007 

Al 1.986 2.013 2.018 

V2-1 1.999 1.992 1.990 

V2-2 1.995 1.987 1.984 

\2-3 1.887 1.896 1.899 

\2-4 2.006 1.994 1.997 

\V4-1 1.984 1.986 1.984 1.985 1.998 1.972 
\4-2 1.995 1.995 1.996 1.998 2.017 1.991 
V4-3 1.988 1.990 1.989 1.987 2.024 1.999 
V4i44 1.976 1.977 1.975 1.975 1.974 1.975 


* Explanation: All crystal segments whose suspension temperatures 
are given in this table were taken from the same melt. The crystal 
designation A2 denotes crystals which were annealed for two hours at 
50° below their melting point and stored exposed to the atmosphere. 
The designation V2 denotes a two-hour annealing period and storage 
over phosphorus pentoxide under vacuum. The designation V4 denotes 
a four-hour annealing period with gradual! cooling to room temperature 
over a period of four hours together with storage over phosphorus pen- 
toxide under vacuum. The Beckmann thermometer graduations are 
such that increases in numerical values are decreases in temperature. 

* After the seventh day, crystals V4—1, V4—2, and V4-3 were ex 
to the atmosphere for twenty minutes and then resuspended. The sus- 
pension temperatures are given in the first column. The second column 
gives the suspension temperatures after exposure to the atmosphere 
for two hours. Crystal V4—4 was not exposed. The third column gives 
the suspension temperatures after a four-hour annealing at 50° below 
the melting point. 


C. Preparation of Suspersion Crystals 


The suspension crystals were grown from a 
melt of purified potassium chloride (see Section 
\) by the method developed by Kyropolous" 


san Kyropolous, Zeits. f. anorg. allgem. Chemie 154, 308 
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which has been described previously. Similar to 
the procedure of the density determination for 
sodium chloride,’ a four-hour annealing period of 
about 50° below the melting point was found 
necessary as shown in Table II. After annealing, 
the crystals were cooled gradually over a period 
of four hours. Phosphorus pentoxide was intro- 
duced at about 150°C into the furnace around 
the platinum crucibles containing the crystals. 
At 60° to 70°C the crystals were removed from 
the furnace and, with platinum-tipped forceps, 
were placed in small glass bottles which previ- 
ously had been heated to about 65°C. The bottles 
were stored immediately in a desiccator over 
phosphorus pentoxide under vacuum. 

The above procedure was found necessary as 
evidenced from the observations recorded in 
Table II. The time required to transfer a crystal 
from its desiccator to the suspension liquid is 
short enough that the crystal is not contam- 
inated. The length of time that a crystal may 
stay in the suspension liquid without change of 
the suspension temperature is about one week, 
and then any change is largely due to decom- 
position of the suspension liquid. 

It will be seen from the data of Table II that 
generally crystals which were exposed to the 
atmosphere for about two hours or more were 
suspended at lower temperatures indicating an 
apparently greater density. This may to a certain 
degree account for the discrepancy between the 
value of the density in the present determination 
and those of Tu® and Kohler.’ The remainder of 
the discrepancy may be due to the lesser purity 
of the potassium chloride used by the ather 
workers and to the lesser precision inherent in 


their methods. (See Tables II and III.) 


D. Data for the Absolute Density 


Fifty-eight crystal segments were suspended 
to determine the suspension temperature of the 
purified potassium chloride. The suspension tem- 
peratures are summarized in Table IV. 

The average suspension temperature of the 
fifty-three unstarred crystals in Table IV was 
1.989° Beckmann with an average deviation of 
+0.000 and an extreme deviation of +0.003 
which is equivalent to —5.6X10-* g/ml in the 
density of potassium chloride. This extreme 
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TaBLE III. Suspension temperatures of crystals taken average suspension temperature, 28.098°C. The 
“ yYarious stages of the purification as outlined in Section corrected values are given in column 3 of Table L 


The average of these values, 1.98651 g/ml, is 


y Crystal Suspension temperatures taken as the best value for the density of the 
— > ae suspension liquid at 28.098°C. The value for the 
| + eo} absolute density of potassium chloride can be 
4-3 1.978 stated as | 
4-4 1.984 
6-3 1.988 When this value is corrected to 25°C by use of 
6-4 1.987 the cubical coefficient of thermal expansion of 
8-1 1.989 potassium chloride,’ the value obtained is 
8-2 1.989 

= pos’c = 1.98721 +0.00002 g/ml 

= 1.98715-+0.00002 g/cm. 


* Explanation: The crystals designated at 4— were grown from 
potassium chloride that had been ipitated four times with hydrogen Ill. MOLECULAR WEIGHT COMPARISONS. THE 


ueous solutions, followed by fusion in platinum, three 


hloride from 
fram tripe dailed water. anda fourhour ATOMIC WEIGHTS OF FLUORINE 
as 6- differed in treatment from the 4- in that they were given six AND CALCIUM 
precipitations with hydrogen chloride from water. Those designated 
a It has been shown that combination of crystal , 
TaBLE IV. Suspension temperatures of fifty-eight densities and x-ray data can be used in the com- ¢ 
crystals which were prepared and preserved as described parison of molecular weights.** With the as- 
sumption of certain atomic weights, other atomic 
Number of crystals suspended at Suspension temperatures weights may be calculated. Combination of the : 
density of potassium chloride obtained in this 
1 1.999% determination and the density of lithium fluoride 
obtained by C. A. Hutchison and H. L. John. 
| 1 1.992 ston!) with the x-ray data of Tu'® and that of 
Straumanis, Ievins, and Karlsons'* makes pos- 
37 1.989 sible the calculation of the atomic weight of 
: 4 . fluorine. The resulting atomic weight is found to i 
1 1.984* be 18.9967 +0.0010. Combination of the data on 
. 1.900" potassium chloride with the density of calcite 
obtained by Bearden" and the x-ray data on 
dase they may be on the calcite of Tu® and Bearden'’ gives 40.0851 
the mean of the fifty-eight crystals suspended. +0.0011 as the atomic weight of calcium. 2 
A detailed account of the calculations of these 
deviation is within the experimental limits of and other atomic weights will appear elsewhere." P 
the determination. The results of these atomic weight calculations *t 
When the experimental densities of the do not support the hypothesis of Zwicky,'® who “ 
bromoform mixture (see column 2 of Table 1) proposed a secondary structure for perfect Ps 
are plotted against temperature, a straight line crystals which requires periodic variations in the Pt 
is obtained. The slope of this line is —1.883 grating spaces in a direction normal to the . 
+0.004X10-* g/ml/degree. Comparison of the crystal planes. th 
Beckmann thermometer with the standard 
thermometer, calibrated by the U. S. Bureau of 4 Y, Tu, Phys. Rev. 40, 662 (1932). é . 
Standards, yielded 1.989° Beckmann equal to on and Karlsons, Zeits. f. physik. 
28.098°C. By employing the value of the thermal ~17J. A. Bearden, Phys. Rev. 38, 2089 (1931). - 
coefficient of density change, the entries in A. Hutchison, obe published 
column 2 of Table I have been corrected to the 211 (1930). 24 
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Statistical Mechanics at Extremely High Temperatures 


G. WATAGHIN 
University of Sio Paulo, Sio Paulo, Brazil 
(Received July 18, 1944) 


Equilibrium conditions between elementary particles and nuclei at temperatures 710° 
degrees are studied. Three temperature intervals below the upper limit ~10" degrees 
[kT~10 ev] are considered. It is found that because of the behavior of some high energy par- 
ticles, known experimentally from observations on cosmic rays and nuclei, some limitations 
arise to the validity of the laws of quantum statistics, in accordance with the idea of the 
existence of a supplementary indeterminacy for high energy particles and a lower limit for 
measurable lengths. Some astrophysical aspects of the phenomena of pair production and of 
the gravitational effect of light particles are discussed. The results concerning pair production 
at temperatures ~10° degrees are summarized in Table I. General formulae for thermal equi- 
librium between nuclei and light particles are given. 


1 


N some previous papers! the statistical phe- 
nomena taking place at extremely high 
temperatures were considered. The purpose of 
the present research is to give a fuller account of 
some of the results obtained. Although the idea 
of the supplementary indeterminacy in the high 
energy regions served as a guiding principle in 
all this research, we tried to base our conclusions 
mainly on the known experiments on cosmic rays 
and nuclei, and on the obvious irreversibility of 
some of the processes in which particles are 
created. 

We shall distinguish three temperature inter- 
vals below the upper limit 7, 137 

1. Temperatures well below the critical tem- 
perature e.g., T<k“mc*/10, for 
which the usual quantum statistical laws hold. 

2. Temperatures character- 
ized by the appearance of positrons and electrons 
produced by thermal photons, at which an 
approximate treatment of the assembly, by 
means of elementary statistical formulae, is still 
possible. In this temperature interval it is 
possible to neglect, as a first approximation, the 
production of few high energy particles with 
E>10% ev. However, the increase of density, 
thermal capacity, and pressure owing to the 
created electron pairs must be taken into 
account. 

3. Temperatures 7 <137 which can be 

1G. Wataghin, Phil. Mag. 17, 910 (1934); Comptes 


rendus 203, (1935); Phys. Rev. 63, 137 (1943); 64, 
248 (1943); 65, 205 (1944). 


introduced only for assemblages of heavy par- 
ticles, e.g., for nuclei. 

4. Conditions [kT>10%ev] at which the 
average amount of energy per particle is so 
great that the majority of the processes become 
irreversible and no equilibrium is possible. The 
lack of the equilibrium is due to the simultaneous 
or successive creation of many particles (espe- 
cially of unstable mesotrons) and to the processes 
involving emission of neutrinos, in which an 
appreciable fraction of the energy and mo- 
mentum escapes all observation. 

We shall not discuss here the problem of 
relativistic invariance examined briefly in a 
preceding paper.? 

2 


Let us consider the equilibrium between 
photons, electrons, positrons, neutrons, protons, 
and nuclei at TAT. Let Mee, 


Nz be the numbers (per cm*) of photons, elec- 
trons, positrons, neutrons, protons, and nuclei of 
charge Z, and atomic weight A, which belong 
to the momentum-interval p, p.+dp,. We assume 
the validity of the laws of conservation of charge: 


conservation of energy: 

pe 
+nz.Eq]+&=E=const. (2) 


[where & indicates the interaction energy, 
E..=mc(1—6*)-* is the total energy of an 


*G. Wataghin, Phys. Rev. 65, 205 (1944). 
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electron, etc. ], and of conservation of the total 
number of neutrons and protons: 


(3) 


We assume also that it is possible to specify 
the states of weakly interacting particles in the 
usual way, namely, by means of eigenstates 
having appropriate symmetry. Then, indicating 


with 
82 Vh*p2dp, = 82 dE, 
the number of quantum states or eigenvalues of 
the energy belonging to the sth momentum- 
interval, we obtain, in the usual way, the fol- 
lowing expression for the number of different 
complexions : 
8 N,-1 ! 
nf + ) g 
(ge—1)!No! mes!(Zs— Mes) ! 
Neglecting the interaction &, and calculating 
the maximum of 


log W—an—BE—yN, 


(4) 


we have, in the usual way: 


log W Nes 
Ones Nes 


+a—BE,,=0, etc., 


and 


N,= gs gs 


exp (Bhy,)—1 exp (—a+BE.)+1 
gs 
Nps = 
exp (a+BE,,)+1 
gs 
nHs= ’ 
exp 
Zs 
NNs = ’ 
exp (y+8Eyn.)+1 
A 
Nzs = 
exp (Za+Ay+BE* )+ 1 


The constant 6 is =1/kT, at least within the 
limits of applicability of the usual concept of 
thermodynamical temperature; the ether con- 
stants a, y can be determined when we introduce 
(5) in (1) and (3). In order to take approxi- 
mately into account the existence of a supple- 
mentary indeterminacy in the region of high 
energy and momenta, the author suggested a 
more general statistics in which the number g, 


(5) 


of quantum states is: g.=8xVh“G(p,)p.dp,, 
where G(p,) is a cut-off factor, which decreases 
more rapidly than p,* for values of p,>p, 
= 137 me [e.g., G(p.) ~exp (— J. 

For temperatures T’<k~'mc’, the following 
approximation is valid : We can treat the protons, 
neutrons, and nuclei by the usual statistical] 
methods, neglect the 8-ray processes and the 
concentration of the mesotrons, put G(p,)=1 
and V=1, and, denoting by m*, n~ the numbers 
of positrons and electrons per cm’, assume: 


Dd =n- — nt 


3x10%., (1 
where N+ is the total number of protons per cm? 
and p, is the density of the nuclear matter (in 
g cm~*) calculated with exclusion of the masses 
of the electron pairs. 

Indicating with x9 = mc?/kT, x=E,/kT, where 
E,=c(p2+m'c)! is the total energy of the 
electron, and putting in (5) 

(kT )3 (x? — x9?) ixdx, 


we obtain easily the following formulae: 


8x (x? — x9?) Ixndx 
n = e-atz4 ’ (6) 

Br 
n= (kT) f 


Substituting (6), (6’), and (1’) and remember- 
ing that: 


sinh a 


cosh x+ cosh a’ 
we can calculate a from: 


mc 
(x? — x9") 


0 cosh x+cosh a’ 


(7) 


or from: 


(7’) 
sin =1. 
zo cosh x+cosh a 


For values of x»<1 and a<i, Eq. (7’) be- 
comes: 


3.29a+0.3308+ (8) 
a~5.2X 10-8 


or 


\\ 


hs 


| | 
1. 
1. 
1: 
1. 
| si 
| 
1 
t 1 
: 
| 
1 
1 
1 
1 
1 
| 
mé 
|| ne 
co! 
; ha 
fic 
ing 
gre 
| 
Ob 
ass 
| | 
ma 


(6) 


STATISTICAL MECHANICS AT HIGH TEMPERATURES 151 


TaBLeE I. Compilation of value calculated for various temperatures. 


Thermal Thermal Radiation Density of 
capacity of capacity of density and energy of 
electrons radiation pressure electrons pairs 
per cm? U =3pr E.* +E.~ 
T —pn(g/cm*) a m [Cr] pair aU /aT =4aT* (erg/cm*) (erg/cm') 
1.2x10° 5000 10-* 4983.0 
12x10? 10-8 479.8 
12x10" 500 10? 498.2 
10" 500 506.8* 
142x108 50 10-¢ 23.0 
1.2 108 50 10 44.6 
1.2108 50 107 93.0* 
5.9x 108 10 10-6 4.8x 10-8 5.81073 5.0 10" 6.2 10" 10% 7.2X 
5.9 108 10 10 2.3 5§.0X10 «9.3K 10 7.2K 108 
1.2X 10° 5 10-¢ 8.7X10- 3.1 1.3X10% 5.010" 1.510% 3.710" 
12X10° 5 1 8.7 10-5 3.1 1.3X10% §.010" 1.5X10® 3.7K 10% 
1.2X 10° 5 102 8.7 X 10-3 3.1 1.3X10% 5.0X10% 1.5X10® 3.710" 
3.0X 10° 2 10-6 550.0 8.0 10" 10" 5.8X 10% 10% 
3.0X 10° 2 10? 6.2xX10-% 550.0 8.0 10" 10" 5.8 x 10% 10% 
3.0x 10° 2 105 6.2107? 550.0 8.0 10" 8.0 10" 5.8x 10% 5.9X 10% 
1.210" 05 10-8 10-% 108 9.0X10% 5.0X10% 1.510% 10% 
12X10" 0.5 10 6.5 10-7 4.2 104 9.0X10% 5.010% 1.510% 2.510% 
1.2 0.5 106 6.5X 10-4 4.2 10 9.0 10'* 5.0 1016 1.5 106 
1.2 10'° 0.5 3x10? 1.9107 4.2 10* 1016 5.0 2.5X 10%¢ 
1.210" 0.05 10° 6.51077 4.2107 10" 10!* 108° 2.6X 10% 
1.210" 0.05 3xX107 19X10] [4.2107 9.0 1019 1.5 108° .6X 10” 
1.210" 0.005 3X10? 1.9X 10-7 4.2x10%] [9.0x10%] [1.5x10"] 10 
* Degeneracy. 


In Table I some values of a, for different 
temperatures and densities, are given. 

The numerical calculations were made on the 
assumption N+=}p,/my, that is, by assuming 
matter contains an equal number of protons and 
neutrons. In the case when the stellar matter 
contains a large fraction of hydrogen atoms, one 
has N+~p,/my and one obtains for the coef- 
ficient in (7’) a value 3.4 10-7. The correspond- 
ing values of a in (8) become nearly twice as 
great. 

From (6) we have: 


n-~—n* 
n-+nt 
(x? — x9") 


o 
e*+e-* (9) 


Observing that the factor [1+(e-*/e+e~¢)] 


assumes values between 1 and 2, one has: 


we, (10) 
~ or —~e~**, 
n-~+nt n~ 


If a<1 the total number of pairs is approxi- 
mately given by 


1/2(n-+n*)~1/2a Nt. 


In the case when a<1, x9<10 the density of 
the electron pairs ppsir Can assume values com- 
parable to or greater than p, (see Table I). These 
values can be calculated from (1’), (8), and (10): 


Ppair= (11) 
n-~~n*~3X 


The question of degeneracy can be treated in 
a manner similar to that of the Fermi statistics. 
We note that the graphics of the two charac- 
teristic functions (e***+1)- can be obtained 
one from the other by means of a translation of 
modulus 2a parallel to the x axis (or by changing 
a in —a@) and that only values of x 2x9>0 enter 
into consideration. Then from (6) follows that 
no degeneracy can occur for the positron gas 
because the values of x 2x» correspond always to 
the quasi-maxwellian region of the distribution 
curve. 

In the low temperature region (large values 
of xo) the negatrons are in a degenerated state 
when a>xXpo. Introducing into (6’) the critical 
value of a: a=xo>1, we obtain in the non- 
relativistic approximation : 


oman 


h 3 
~xo (12’) 


*dp,, 
> Py 
wing 
tons, 
tical 

the 
)=1 

cm? 
(in 
iSSes 
here j 
the 

| 
| 

( | 


152 GLEB WATAGHIN 


One can easily verify that this relation coincides 
with the usual Fermi criterion for degeneracy. 

At extremely high temperatures, x»<1, the 
values of a 2x correspond to values of 


Pn 2Nt+my~ > per, 


where Ag=h/2xmc and X10? 
g cm“ is the critical value of the nuclear density. 
Indeed, the approximate value, for x»—0, of 
n-~nt, is: 


mc \* © 
f 

h o e*+1 

~ = 1.736X 10*]. 


For such values of p, the average distances 
between the nuclei (and also between the elec- 
trons) become <Ao, and the assumptions of the 
classical or the quantum electrodynamics, on 
which the preceding calculations are based, are 
no longer valid. Thus we exclude such cases from 
our consideration. At x»<1 degeneracy of nega- 
trons does not occur at densities pra < per. 

The last columns and lines in Table I are 
reproduced in order to illustrate some interesting 
features of pair production and radiation proc- 
esses in regions where we cannot be sure of the 
validity of the preceding simple statistical cal- 
culations. 


3 


Obviously we must expect the failure of the 
preceding simple calculations at temperatures 
and densities above certain limits for several 
reasons. Let us examine some of them. 

(1) We neglected the nuclear reactions, which 
can give rise to 8-ray processes connected with 
the emission of neutrinos. In these processes the 
laws of conservation of energy and momentum 
are out of our direct experimental control as 
long as we have no means to measure the energy 
of the neutrinos or ‘at least to observe them. In 
some cases, when we can neglect 8-ray processes, 
an approximate treatment of nuclear reactions 
is possible on the basis of the formulae (5). One 
can calculate the equilibrium concentration of 
the protons, neutrons, and nuclei determining 
the constants a, 8, y from (1), (2), (3). The 
calculation of y leads to especially interesting 
results, which we shall mention briefly later. 


(2) At nuclear densities pra per~3 X10" g 
the average distance between the nuclei becomes 
of the order of the Compton wave-length. At 

10" degrees [x050.05] the density of the 
electron pairs n-+n+ becomes 52X10"~A 
In both cases the linear equations of Maxwell 
are no longer valid. The field strengths become 
of the order of or greater than | E| .,=e/p*, where 
l~e?/mce?~h/2xmc is the universal length. The 
first non-linear correction terms in the La. 
grangian are known to be of the form: 


a, (E*— B*)?+6,(E-B)?, 


where a, 8; are constants calculated by several 
authors. The representation of a field by means 
of stationary waves or quantum states of the 
linear field theories also becomes impossible at 
least in the region of high frequencies. The writer 
has suggested* that in such conditions our pos- 
sibilities to distinguish eigenvalues, e.g., measure 
“observables,”’ is limited in a new way, which 
depends essentially on a new kind of phenomena 
appearing at high energy collisions, such as 
pair production, 8-ray processes, and mesotron 
showers. These phenomena are unavoidable in 
high energy collisions and thus appear neces- 
sarily in the interaction between the measuring 
devices (which are at rest in the respective refer- 
ence frames) and the objects of observation. The 
easiest way to obtain converging results in the 
usual calculations of cdllision processes, problems 
of self-energy, and interaction energy, is to 
abandon the point source models and to use the 
so-called cut-off prescriptions. In the foregoing 
formulae (5) these cut-off factors are introduced 
by means of the factors G(p,) appearing in g,. 

(3) Recently the writer pointed out the irre- 
versible character of the collision processes which 
give rise to the production of several particles, 
and which are accompanied by 8-ray processes. 
This irreversibility (due to the fact that the 
inverse processes involve ternary or multiple 
collisions and emission of neutrinos) is in con- 
trast with the validity of laws of the type (5) at 
T>10" degrees, but agrees well with the 
modified formulae obtained from (5) introducing 
the cut-off factors G(p,). 

For instance, Planck’s formula could not be 


3G. Wataghin, Nature 142, 393 (1938) ; Comptes rendus 
207, 358, 421 (1938). 
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STATISTICAL MECHANICS 
valid at temperatures of the order of 10" degrees 
[kT ~~ 10° ev ]. At these temperatures the number 
of photons of energy >10° ev should be 5X 10° 
times the number of those having hv<10° ev. 
From cosmic-ray observations we know that 
each of the photons or electrons having E>10%ev 
gives rise to a multiplication shower in which 
the initial energy is subdivided among a great 
number of created particles. The rate of the 
inverse processes from many low energy par- 
ticles to a high energy one is negligible. Indeed, 
the annihilation of an electron pair can occur 
with the creation of two photons or of one 
photon. In the first case the number of light 
particles is not changed (besides, the cross 
section for annihilation with two-photon pro- 
duction is maximum at E~10* ev and decreases 
rapidly with the energy). The probability of a 
three-particle-collision (two electrons with a 
nucleus), necessary for a one-photon anni- 
hilation, is much smaller than that of a collision 
of one electron with one nucleus. Only these 
three-particle collisions give rise to a reduction 
of the number of particles, and the rate of these 
is necessarily very low. We conclude 
that the equilibrium is impossible because the 
processes of multiplication of particles cannot be 
balanced by processes which reduce their 
number. 
(4) Let us consider a region near the center of 
a star during the contractive evolution which can 
give rise to extremely high temperatures. From 
Table I (last two columns) one can see that at 
temperatures 7 ~10" degrees, the density of the 
radiation energy and the density of the electron 
pairs are much greater than p,, so that one can 
neglect the nuclear density. If such conditions 
are verified for a sufficiently extensive region 
[R>=10'° cm], a very peculiar situation arises. 
The mass of the radiation and of the electron 
pairs increases ~7*. The gravitational energy 
of this mass varies as — 3(44/3)'®GM®*%p!/8, At 
low densities the gravitational energy is neg- 
ligible in comparison with the total relativistic 
energy of the nuclei, radiation, and electrons. 
But if praatppair is >10’ in a region of linear 
dimensions ~10'° cm, the gravitational energy 
becomes ~Mc*. The gravitational pressure 
varies as $(42/3)'/*GM?/%p*/8, whereas the pres- 
sure of the relativistic electron gas and the 
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radiation pressure vary ~ p. Thus at sufficiently 
high values of T and (praa+Ppair) the gravita- 
tional pressure becomes and remains greater than 
the radiation pressure (in the outer layers of the 
high density region) and determines a con- 
traction of the gas. These conditions are verified 
when the negative gravitational energy 


—§GM?/R= — 


becomes comparable to the total energy of the 
star Mc*. The increase of the gravitational mass 
of the star at expense of the potential gravita- 
tional energy implies a change of the gravitational 
field due to the star, and thus the variation of the 
mass must be accompanied by a gravitational 
wave which produces this change of the field. A 
similar variation but of opposite sign is pro- 
duced when energy is carried away by neutrinos 
in 8-disintegration processes. Obviously the con- 
traction considered above cannot continue indef- 
initely and must be counterbalanced by losses 
of energy (and mass) by radiation, electrons, and 
neutrinos. 

There exists also the effect of the gravitational 
red-shift, which produces an apparent decrease 
of the energy of the electrons and photons and 
which becomes appreciable in the considered 
case: 


The paradoxical aspect of the phenomenon in 
which gravitational mass of electron pairs is 
created by a gravitational field induces us to 
admit that on the condition of such strong inter- 
action between the gravitational field and other 
fields some new physical process counteracts the 
transformation and makes it impossible. The 
most plausible assumption seems to be the one 
that the neutrino emission prevents the increase 
of the mass and of the gravitational contraction 
by means of induced 8-ray processes. 

In this connection it seems noteworthy to us 
that accepting Friedman's idea of the expanding 
universe, one can calculate the increase of the 
gravitational energy of the receding nebulae by 
means of the Newtonian potential. Assuming 
that some 10° years ago the whole mass of the 
observed galaxies was concentrated in a reduced 
volume of linear dimensions ~1 parsec, one can 
make a very rough calculation of how much 
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gravitational energy the observed galaxies have 
acquired, and one finds that the average energy 
per proton is ~10!’ ev. Thus, if we accept the 
idea of Weizsacker, Chandrasekhar, and Henrich‘ 
on the prestellar stage of the universe and assume 
that the gravitational energy of expansion was 
produced at expense of the particle-energy, we 
must say that no equilibrium could possibly 
exist at an earlier epoch with such an amount of 
energy per particle.® 

But the analysis of the abundance and dis- 
tribution of the isotopes of the chemical elements 
in the universe strongly suggests the idea that 
nuclei were formed in the whole universe in 
similar conditions. We assume that such uni- 
formity of distribution could be attained only 
in conditions which approximate a statistical 
equilibrium. In order to study this possibility, 
we can apply the formulae (5) to the equilibrium 
among nuclei, electrons, and photons. The cal- 
culations are similar to those used in the con- 
struction of Table I. Here we want to mention 
only some general results. 

At temperatures T>10!, the constant y is 
determined essentially by the concentration of 
neutrons and protons. The concentration of 
heavier nuclei is entirely negligible and only the 
concentration of ‘He2 is appreciable. The con- 
centration of the heavier nuclei increases gradu- 
ally in the temperature region between 3X 10° 
and 10° degrees (x»~2 or 5). The analysis of the 
nuclear structure (mass defect as function of Z 
and A-—Z) shows that the stability of the 
nuclei is closely related to the probability of the 
B-disintegration processes, and the formation of 
stable nuclei follows closely the evolution of the 
concentration of electron pairs, photons, and 
neutrons. The formation of heavier nuclei by the 
decrease of the temperature from 6X 10° to 10° 
degrees follows in the order of increasing atomic 
weights from lightest to heaviest nuclei. An easy 
calculation shows that at temperatures of a few 
billion degrees the nuclear photo-effect becomes 
of great importance and gives rise to the emission 
of many neutrons, so that equilibrium is possible 


4 Weizsacker, Chandrasekhar, and Henrich, Astrophys. 
J. 95, 288 (1942); Physik. Zeits. 39, 633 (1938). 
5G. Wataghin, Phys. Rev. 63, 137 (1943). 
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with a high concentration of neutrons. At tem. 
peratures below 10° degrees the neutron concen. 
tration vanishes and the usual cyclical thermo. 
nuclear reactions take place. All these concly. 
sions derive from the formulae (5) corrected by 
the introduction of the cut-off factors G(p,). 

The limits of applicability of these formulae 
depend on the induced #-ray processes. In ac. 
cordance with a fundamental idea of G. Gamow 
and M. Schoenberg related to the so-called 
‘‘urca-processes,”’ at temperatures > 10"° degrees, 
we can expect the 8-ray processes to acquire such 
an intensity that the energy losses due to the 
neutrino emission must invalidate the possibility 
of equilibrium, and thus the law of conservation 
of energy (2) becomes out of our control. We 
think that this fast loss of energy by neutrinos 
was at a certain epoch in competition with the 
expansion of stellar matter in the non-equi- 
librium process which increased the gravitational 
energy of the universe. 

Leaving the highly speculative consideration 
on the prestellar stage of the universe, we want 
to call attention to the fact that 8-ray processes 
introduce a new limitation to the concept of 
temperature in an assemblage of nuclei and light 
particles because the non-applicability of the con- 
servation law (2) implies the impossibility of 
equilibrium and of identification of the constant 
8B with (k7)—. Of course, one can introduce the 
notion of temperature also for a nucleus as is 
done usually in order to measure the average 
energy per nucleon. But even in nuclear matter 
the equilibrium becomes impossible when the 
average energy per particle is greater than the 
rest energy of the meson yuc?~10* ev because at 
higher temperatures irreversible collision proc- 
esses take place with production of unstable 
mesons. Thus the temperature 7,,~10" degrees 
must be an upper limit also for heavy particles. 

In this paper we shall not discuss the problems 
of compatibility of the supplementary indeter- 
minacy with the relativistic invariance. Also the 
production of mesotron showers by collision of 
high energy nuclei and correlated phenomena 
were not considered. An account of these 
problems will be published elsewhere in a short 


time. 
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Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


An Attempt at Finding the Relationship Between 
the Nuclear Force and the Gravitational 
Force 


K. C. WANG AND H. L. Tsao 


Department of Physics, National University of Chekiang, 
7” Meitan, Kweichow, China 


June 2, 1944 


HE best we know about the nuclear force is that it is 

of short range type. Its analytic form isstill unknown. 
The meson theory of nuclear force is, up to the present at 
least, inconclusive.t On the other hand, the gravitational 
force is explained by the theory of general relativity. If we 
temporarily disregard the theory of general relativity and 
pay attention to the important fact that both the nuclear 
force and the gravitational force are attractive, we may 
assume that these two forces are of the same origin. Let 
us tentatively take the function 


V=—Aek!r (1) 


as the mutual potential between any two nuclear particles, 
i.e., between two protons, two neutrons, or one proton and 
one neutron, with A and K as constants. Then, as the 
distance r between the particles becomes much larger than 
K, Eq. (1) turns out to be 


(2) 


In this expression, the first term on the right-hand side 
is a constant, and therefore is of no consequence. Thus 


V(r) =—AK/r. (3) 


On the other hand, the gravitational potential between 
the two nuclear particles of masses M and M’ should be 
given by 


V=—G(MM'/r). (4) 


NUMBERS 5 AND 6 


SEPTEMBER 1 AND 15, 1944 


If the nuclear force and gravitational force are really of the 
same origin, we must have, from (3) and (4), 


AK=GMM'=1.83X10-* erg cm, (5) 


since G=6.66X10-§ and M=M’=1.66X10™. Recently, 
Ragan, Kanne, and Taschek,? from the proton-proton 
scattering experiment, found that the potential between 
two protons can be represented by a rectangular well of 
width 2e*/mc* and depth 10.5 Mev. Now if from our equation 
(5), we put K as the Compton wave-length of electron, i.e., 
K=h/mc =3.83 cm, we get A =4.80X10™, and if 
we substitute these values of K and A into (1), we find the 
potential V=10.5 Mev at r=1.5le*/mc* which is only a 
little smaller than the experimental value of 2e*/mc*. In 
fact, if the potential is not a straight rectangular well but 
a slightly curved one, for the same depth of potential, a 
smaller value of r is expected. If, instead of Eq. (1), we take 


V=—(B/r)eX!* = —(AK/r)eX!" (6) 


as the nuclear potential function, we get the result V= 10.5 
Mev at r=1.62 e*/mc* which checks even better with the 
experimental value. Thus we may tentatively state that: 


V=—(B/r)eX'" forr>a (7) 
V=—(B/a)eX for rSa 


for a certain value of a, the exact value of which is of no 
great importance at present. Since the value of K is not 
purely arbitrary, but is chosen as the Compton wave- 
length of electron, the agreement between the experi- 
mental and the calculated value is perhaps not accidental, 
and may indicate that the electron, as has long been sug- 
gested, plays an important role in nuclear interactions. If 
it is required to avoid the rather artificial way of “cutting 
off" of the expression (7) at r=a, we give an alternative 


expression 
K /{K’'\* 
V=—A exp | (8) 


with n lying between 1 and 2, and K’ a new constant which 
may tentatively be put as the Compton wave-length of 
meson, i.e., K’=h/m'c ~h/180mc, but this makes the width 
of the potential well smaller than 1.51 e/mc*. It will be 
seen in a paper by M. S. Wang? that good results for the 
binding energy of deutron and the scattering cross section 
of protons by fast neutrons can actually be obtained by 
calculations in employing either the potential function (1) 
or (6). In conclusion, we express our gratitude to discussions 
with our colleague, Professor H. P. Soh. 


1H. A. Bethe, Phys. Rev. 57, 260, 390 (1940). 

2G.L. , W. R. Kanne, and R. F. Taschek. Phys. Rev. 60, 628 
(1941). See also T. H. Osgood, J. App. Phys. 14, 53 (1943). 

3M. S. Wang, Phys. Rev. 66, 103 (1944). 
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Proton Primaries of Atom-Annihilation 
Cosmic Rays 


Dana T. WARREN 
Department of University of Omaha, 
August 7, 1944 
ILLIKAN, Neher, and Pickering have presented the 
atom-annihilation hypothesis of cosmic-ray origin 
in two recent papers.' In I, footnote 4, they point out that 
the magnetic effects depend only on the energy of the 
primary particle, and consequently their hypothesis would 
not distinguish one kind from another. In II, page 245, they 
insist that their results would be valid only for electron 
primaries, not for mesotrons or protons. Their reasoning 
in this case depends on a determination of the kinetic 
energy available when a helium atom decomposes into two 
electrons or into two protons, respectively. The latter, 
of course, gives about half as much kinetic energy, from 
which they conclude that the latitude of appearance 
would be appreciably higher. 

In view of the recent experimental indications? and 
theoretical calculations’ concerning the protonic nature of 
the primaries, it seems important to point out that I is 
correct, and not II. This is because the magnetic effects 
depend, not on kinetic energy, but on total energy, which 
is of course the same for both kinds of disintegrations. 
This can readily be seen, for instance, from the table of 
values that relates the magnetic deflection variable (called 
r- in I, Fig. 2) in stoermers with the energy in electron 
volts.*5 A given value of the kinetic energy corresponds to 
a larger value of r, for protons than for electrons, and this 
effect very closely neutralizes the decrease in mass annihi- 
lated. A detailed calculation can be made on the basis of Eq. 
(12) of reference 4. For the case of electrons, this reduces 
to Eq. (66) of reference 5. On the basis of this relation are 
drawn the values scaled in I, Fig. 2 (copied from Fig. 14 
of reference 5), specifically labeled electron energies. For 
protons, however, the equation takes the form 


when rf, is in stoermers and 6 is the kinetic energy in Bev. 

When these equations are applied to the values given 
in II, page 245, for electrons (1.88 Bev) and protons (0.94 
Bev), we obtain r,=0.178 stoermer for electrons and 0.165 
stoermer for protons. This difference of 7 percent in the 
value of r. would correspond to a latitude difference of 
about 1°, instead of the 6° stated in II. Similar calculations 
for the other atom-annihilation bands give differences of 
1 percent or less, which would be quite unobservable. 

We conclude, therefore, that the atom-annihilation 
hypothesis is quite independent of the nature of the 
primary cosmic radiation involved. This does not apply, 
however, to alpha-particles, with their double charge. 


V. Saher, ond W. Ph 


. Rev. 6 
397 7 ion) Ch (hereafter ca called 1); Phys. Rev. 63, 234 (1943) called ip.’ 
ev. Mod. Phys. 11, 208 (1938 


w. and H. W. Phys. 6, 78 (1943). 
Lemaitre and M. S. Vallarta, P ev. 43, 87 (1933). 
5 T. H. Johnson, Rev. Mod. Phys. 10, 193 (1938). 
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Efficiency of Production of Continuous 
X Radiation 
KIRKPATRICK 
Stanford University, Stanford University, California 
August 2, 1944 


OR over ten years we have had in the work of Sommer. 
feld and others the basic structure of a proper theory 
of continuous x-ray production, but little has been done 
to put the predictions of theory into numerical forms 
which could be compared with experiment. We have now 
completed computation of the elementary electric dipole 
moments for radiation between the quantum limit and 
1.431A resulting from collisions of electrons in the energy 
range from 15 to 90 kev with atoms in the range from 
Z=13 to Z=92. The results of these computations, 
together with a description of the computational procedure, 
will shortly be submitted for publication. 

From these computations one readily obtains W, the 
theoretical total energy of continuous spectrum x radiation 
in unit frequency interval per bombarding electron per 
atom per cm? of target material. This is a numerical result 
known for a variety of bombarding potentials V, radiation 
frequencies, and target atomic numbers Z. By graphical 


integration of W over all frequencies from zero to the © 


quantum limit, the total continuous x-ray energy per 
electron per atom per cm? is obtained. This quantity is, 
of course, a function of Z and V and it has been evaluated 
for a variety of values of these parameters. 

To deduce the efficiency of such x-ray production one 
needs to know the total energy lost per atom per cm? by 


TABLE I. Theoretical efficiencies expressed in percent. 


Z=13 Z=47 Z=70 Z=92 
Efficiencies 


0.20 10.53 0.042 0.137 0 0.555 
0.30 24.66 0.077 0.270 0.474 0.685 
0.40 46.52 0.132 0.438 0.675 0.895 
0.50 79.00 0.179 0.645 0.971 1.25 
0.526 90.00 0.182 0.693 1.06 1.34 


an electron of velocity ratio 8 while penetrating the target 
material. The formula of Williams' as modified by Webster 
may readily be put in the form: 
Energy loss per electron per atom per cm? 
= 5.62 
By simple division we get the efficiencies reported at the 
Berkeley meeting of the American Physical Society.* 
The total continuous spectrum energy R emitted during 
the complete stoppage of an electron may be evaluated by 
a further graphical integration of the form 
where E is the electron kinetic energy and 8» is the maxi 
mum or initial value of 8. The first derivative in the 
integrand is, in absolute value, the calculated efficiency 
aforementioned; it is, of course, a function of 8. The 
second derivative in the integrand is a function of 8 
obtained by differentiating the relativistic kinetic energy 
expression. The integration was carried out for some 
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twenty cases and the over-all efficiency of continuous 

x-ray production obtained by dividing R by the initial 

electron energy. These theoretical efficiencies, expressed in 
t, are shown in Table I. 

These values are to be compared with observed effi- 
ciencies. Compton and Allison have summarized existing 
work in the equation 

Efficiency = 1.1 10~Z V, 


in which V is in volts and the efficiency is a fraction 
rather than a percentage..The calculated efficiencies all 
increase with increasing Z and V but not in all cases 
lly as the experimenters have found. The 
increase with Z for the three higher voltages is an accurate 
proportionality in each case but at lower voltages an 
upward concavity appears. As to the increase of calculated 
efficiency with V the curve is linear for the higher values 
of Z but misses the origin; for lower Z the curves go 
through the origin but are slightly concave downward. 
On the whole, the agreement between theory and experi- 
ment as to the functional relation between efficiency an& 
Zand V is as good as’could have been wished, in view of 
the difficulties of such observations and the approximations 
(neglect of relativity and screening) of the theory. 
The theoretical efficiencies are from 50 to 100 percent 
ter than those given by the summarizing equation of 
Compton and Allison, but they lie within the spread of the 
results of the several observers. A discrepancy in this 
sense should be expected since the theoretical emission 
extends to zero frequency and thus includes a minor, 
though not negligible, contribution from a spectral region 
unavailable to measurement because of heavy absorption. 
This agreement encourages the hope that all the observable 
features of thick-target x radiation may be found explicable 
by fundamental theory. 
1Compton and Allison, X-Rays in Theory and Experiment (1935), 
>’. Kirkpatrick, Bull. Am. Phys. Soc. 19, 7 (July 22, 1944). 


Further Remarks on the Difficulties of the 
Meson Theory 


D. IWANENKO 
The State University, Sverdlovsk, U.S.S.R. 
August 7, 1944 


OME time ago Oppenheimer and his collaborators! 
arrived at the conclusion that mesons should be 
described by pseudoscalar equations and possess zero 
spin because vector equations for the mesons of the spin 
unity lead to unadmissibly high cross sections at various 
scattering effects. This point of view seems to be supported 
in recent publications of Pauli and Dancoff and Fréhlich.* 
We may remark that pseudoscalar mesons behave satis- 
factorily only for electromagnetic scattering (e forces), but 
in cases of scattering under the influence of specific nuclear 
@ forces, there is no essential difference between behavior 
of vector and pseudoscalar particles. For instance, the 
cross sections for the scattering of pseudoscalar mesons on 
neutrons grow unlimitedly at high energies. The very 
reason for this behavior lies in the presence in the commu- 
tation rules for field quantities of the terms of dipole type 
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having the form »*. This can be deduced on the lines of 
the analysis applied previously by us to the vector case.* 

Most immediately this can be inferred from the fact 
that pseudoscalar mesons, if transferring nuclear inter- 
actions, lead to the forces of the dipole type, as vector 
mesons do, and it is just this coupling of the dipole type 
between nucleons and mesons which leads to the high 
cross sections. 

In our opinion there is no essential ground for preferring 
the pseudoscalar description of mesons as this also needs 
further amendment. 

As to the removing of these dipole difficulties, we remark 
first of all that the reaction of the mesonic field on the 
radiating nucleonic dipole is described in the classical 
picture firstly by a Heisenberg term containing the second 
time derivative of the nucleonic dipole and secondly by a 
term with the third time derivative, just as in the case of 
the radiation damping force acting on a vibrating charge. 
This latter term was introduced by Bhabha and the 
present writer.‘ It is important to emphasize that Heisen- 
berg’s term will be absent for the point source models of 
the nucleonic dipole, so that the whole reaction of damp- 
ing is expressed then by the term with third derivative. 
In the quantum theory the damping was introduced 
simultaneously by Heitler, Wilson, and Sokolow,' and 
these authors have shown that upon taking the damping 
into account the cross sections for the scattering of vector 
mesons on nucleons and of the scattering of light on 
mesons remain limited at high energies. Further calcula- 
tions of Sokolow and myself show that the introduction 
of the damping for pseudoscalar equations also removes 
the difficulty of the growing cross sections for the scattering 
under the influence of specifically nuclear forces. Pre- 
liminary calculations show also that the unpleasant 
divergent part of the polarization of vacuum due to virtual 
mesons remains finite.® 


F. ond 3. Kusaka, Phys. Rev. 59, 

2W. Pauli and S. M. , Phys. Rev. 62, 85 (1942); Frdhlich, 
Phys. Rev. 62, 180 (1942). 

3D. Iwanenko and A. Sokolow, Phys. Rev. 60, 277 (1941); J. Phys. 
Moscow 5 (1941). 

‘D. ong A. Sokolow, J. P 3. 12 (1940). 

tA kolow, J. P’ ys. Moscow 6 agai) Heitler, . Camb. 
Phil soe 37, 291 Wilson, Proc. Soc. 37, 301 (1941). 

* As in Phys. Rev. 60, 277 (1941), the tion of 
vacuum via mesons was calculated in collaboration with M. Rosman, 
now in active service. 


Errata: On the Perturbation of Boundary - 
Conditions 


[Phys. Rev. 65, 307 (1944)] 
HERMAN FESHBACH 
Department of Physics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 
HE abbreviated derivation of Eq. (3.7) and the 
corresponding equations (4.1) and (5.6) in the paper 
“On the perturbation of boundary conditions” has per- 
mitted some misinterpretation. 

Formula (3.7) applies only when ¥, has the éigen- 
value k*. If the eigenvalue of y, is K,*, then the term 
(k8—K,*) f¥.Gidv should be added to right-hand side of 
(3.7). The same term should be added to (4.1). The 
corresponding term of Eq. (5.6) is (K,*—&*) f/A,-Gido. 
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PHYSICAL REVIEW 


HE 262nd meeting of the American Physical 
Society was held in Berkeley, at the 
University of California, on Saturday, July 22, 
1944. In spite of the continually decreasing 
number of graduate students, and the continued 
preoccupation of most physicists in war work, 
the attendance, which varied from 65 to over 
100, was the largest at any Pacific Coast meeting 
since the start of the war. About 35 members of 
the Society and guests had luncheon together at 
the Faculty Club. Brief remarks were made at 
the luncheon by Dr. K. K. Darrow, Dr. E. U. 
Condon, Dr. Otto Stern, and Dr. Paul Kirk- 
patrick. 
As in the case of the preceding Pacific Coas 
meeting, very few contributed papers were 
anticipated, and accordingly arrangements were 
made for four invited papers. These papers 
undoubtedly added greatly to the interest of the 
meeting and to the size of the attendance. But 
the unexpected receipt of 18 contributed papers, 
none of which was read by title, resulted in a 
program definitely too long for a single day. 


VOLUME 66, NUMBERS 5 AND 6 


Proceedings of the American Physical Society 


MEETING AT BERKELEY, CALIFORNIA, JULY 22, 1944 


Abstracts of Contributed Papers 


SEPTEMBER 1 AND 15, 1944 


The morning session opened at 9:30 A.M. and 
adjourned at noon. The afternoon session opened 
at 1:30 p.M. and adjourned at 5:20 P.M. During 
the morning session two invited papers were 
presented, viz., 

Specific Ionization and Energy Loss of Cosmic Rays, 
Wayne E. Hazen, University of California. 


Microwave Oscillation Generators Using Velocity Mody. 
lated Electron Beams. Epwarp U. Connon, Westinghouse 
Research Laboratories. 


During the afternoon session the following 
two invited papers were presented: — 


* Applications of Radioactive Tracers to the Ph 

of Plants. Roy OVERSTREET, Experiment Station, College 
of Agriculture, and Radiation Laboratory, University of 
California. 

’ Some Aspects of Our Present Understanding of Ele. 
mentary Particles. JosepH W. WEINBERG, University of 
California. 


Abstracts of the eighteen contributed papers 


follow. 
R. T. BIRGE 
Local Secretary for the Pacific Coast 


1. Cascade Showers in Lead. Wayne E. Hazen, Uni- 2. Growth of Water Drops in a Wilson Cloud Chamber, 
versity of California.—The calculations for number of E. O. Barrett, University of California.—Water drops 
electrons and positrons in a shower as a function of depth formed in a Wilson cloud chamber were photographed 
have been carried out lately to an accuracy of a few with periodic illumination and their velocity was measured, 
percent.! These latest calculations give curves for number The drop size was determined by means of Stokes’ law and 
of particles versus depth that are practically identical in the square of the radius of the drop was found to vary 
shape with previous curves.? The important difference is linearly with time as soon as the drop was large enough to 
that the corresponding curves are for initiating particles register photographically. A 60-cycle neon light that moved 
differing in energy by a factor of ~1.4. In photographs of | with constant velocity during the exposure was turned on 
a cloud chamber that contained eight lead plates, 50 by the expansion valve. This light was photographed on 
complete cascade showers have been observed. The ob- the same frame and its flashes were used to establish the 
served showers can best be compared with theory by zero of the time scale for the drop images. It was observed 
plotting number of particles at the shower maximum versus that the straight line obtained from the experimental 
total number of particles in the shower. Theory predicts points intersected the time axis about 0.1 sec. after the 
essentially a straight line with slope 0.2. All but one of the zero. This fact seems to indicate that the rate of growth of 
experimental points lie above this straight line. The the drops is smaller during the first 0.1-0.2 sec. and then 
observed showers contain only one-half the predicted total reaches a constant value. The slope of the straight part 
number of particles for a given size of shower maximum. of the curve was found to be 7.8 X 10~* cm?/sec. 
Since a large number of the electrons have energies as low 
as three Mev, it is possible that scattering and increase in 3. The Average Specific Ionization of Cosmic-Ray 
ionization for E<1.5 Mev would account for the dis- Mesotrons in Hydrogen. LesTER L. SKoLIL, University 
crepancy. of California.—Counter-controlled cloud-chamber photo 
1H. J. Bhabha and S. K. Chakrabarty, Proc. Ind. Acad. Sci. AIS, graphs of cosmic-ray mesotron tracks in hydrogen have 


ba yo and K. Greisen, Rev. Mod. Phys. 13, 240 (1941). been taken. The chamber was operated in such a way that 
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the average specific ionization could be obtained from the 
number of drops that condensed on positive ions. A triple- 
coincidence counter telescope was employed, with two 
counters below the cloud chamber and one above. Sufficient 
Jead was placed above and below the top counter so that 
electron tracks could be distinguished from mesotron 
tracks. The average specific ionization, excluding energy 
transfers greater than 800 ev (clusters of more than thirty- 
two drops), is essentially 10.5 ion pairs per cm of hydrogen 
at N.T.P. 


4. Radiofrequency Dielectric Properties of a Cellulosic 
at Low Moisture Contents. W. C. Dun ap, JR., 

anp B. MAKOWER, Western Regional Research Laboratory. 
—As background for a study of electrical moisture testing 
by radiofrequency methods, measurements were made of 
the radiofrequency dielectric properties of dehydrated 
carrots, a material typical of those on which such methods 
are used. Measurements were made of the dielectric 
constant « and specific conductivity @ (a.c. and d.c.) as 
functions of moisture content (1.5-21.6 percent), tempera- 
ture (1.5°C-39.8°C), frequency (18 kilocycles-5 mega- 
cycles), density (0.792-1.45 g/cm!) and particle size (small 
and medium). « increases little with moisture content up 
to 6-8 percent; above this region ¢ increases rapidly. 
Increase of temperature increases both the value of « and 
the rate of change of ¢ with moisture. At 1.5°C there is no 
detectable change of total ¢ with increase in moisture until 
6-8 percent is reached; above this region the variation 
becomes exponential. At higher temperatures this apparent 
discontinuity disappears. Particle size has little effect upon 
«or ¢, for a given bulk density. Increase of bulk density 
increases ¢ and total o uniformly on a logarithmic scale for 
all frequencies. A brief discussion will be given of the theory 


of the properties observed. 


5. Control of Glycol Vapor Concentration in Room Air. 
C. E. NrELSEN, Department of Physics, AnD K. B. DEOMeE, 
Division of Veterinary Science, University of California.— 
The vapors of certain glycols, non-toxic to man, are 
effective in killing air-borne bacteria. The effectiveness 
depends critically both upon the amount of water vapor 
and upon the amount of glycol vapor present. In conse- 
quence, satisfactory clinical experiments cannot be made 
without stabilizing both humidity and glycol vapor pres- 
sure. Although humidity is controllable, no method of 
controlling glycol vapor pressure under room conditions 
has been devised. The Wilson cloud chamber offers an 
approach to this problem. When room air is cooled by 
adiabatic expansion, droplets will grow upon nuclei nor- 
mally present if the temperature falls below the dew point. 
The addition of glycol vapor should permit condensation 
with a lesser expansion, glycol in the droplets lowering 
their vapor pressure and causing them to grow in an 
atmosphere not saturated with water vapor. Reduction in 
the expansion required, water vapor pressure being main- 


tained constant, should thus indicate the vapor pressure - 


of glycol. Experiments indicate that glycol does reduce 
the expansion required for condensation, and that the 
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occurrence of condensation can serve to control the rate of 
addition of glycol to maintain constant glycol vapor 
pressure. 


6. On Electron Microspectroscopy. L. MARTON, Division 
of Electron Optics, Stanford University.—The energy losses 
of fast electrons corresponding to K levels of light elements, 
and observed first by G. Ruthemann,"* were used by J. 
Hillier* in an instrument designed for the investigation of 
the composition of small specimens. His instrument, called 
an “electron microanalyzer,” consists of an inverted elec- 
tron microscope producing a small probe on the surface of 
the specimen, with analysis of the transmitted beam and 
image formation by the shadow method. Our investigations 
show that the energy losses can be observed in a trans- 
mission type of electron microscope. For such a purpose an 
attachment is built on to the Stanford electron microscope 
replacing the transparent fluorescent screen used for end-on 
viewing the image. The attachment consists of a disk, 
coated with fluorescent material, and provided with a 
small slit through which part of the highly magnified 
image can be directed into a magnetic analyzer extending 
belaw the microscope. A simultaneous observation of both 
the highly magnified image and its velocity spectrum is 
therefore achieved. The area used for the velocity analysis 
can be made variable by changing the end magnification 
and keeping the slit width constant. The magnetic analyzer 
is of 135°, direction focusing type. Observation, temporarily 
at least, is on a fluorescent screen inclined toward the 
beam axis and approximately in the correct position for 
the velocity focal plane. Some of the obtained results are 
discussed. 

1G. Ruthemann, Naturwiss. 29, 648 (1941). 


2G. Ruthemann, Naturwiss. 30, 145 (1942). 
3 J. Hillier, Phys. Rev. 64, 318 (1943). 


7. X-Ray Monochromatization by Four Balanced Filters. 
PaAuL KIRKPATRICK AND C. K. CHANG, Stanford University. 
—A limitation of Ross’ balanced filter method is the 
physical impossibility of entirely suppressing a—} outside 
the pass band, a(A) and b(A) being the fractions transmitted 
by filters A and B, respectively. This limitation may be 
surmounted by using two similar filters of element A and 
two of element B. From an x-ray source adjacent beams, 
defined by similar diaphragms, are allowed to pass to a 
common ion chamber. An ion current observation is made 
with the superposed A filters in one beam and the super- 
posed B filters in the other. One A and one B filter are 
then transposed and a second current reading is taken. 
The difference between these currents is then a summation 
to which each wave-length contributes an amount propor- 
tional to (¢—5)*. In Ross’ method the corresponding 
magnitude is (a—b). Thus the contrast between the pass 
band contribution (where a/b may be of the order of 3) 
and the contribution from other spectral regions (where 
a/b=1) is greatly accentuated by the four-filter method 
and a monochromatic effect may be achieved with filters 
which have been only very roughly balanced. 


8. Linearity of d.c. Vacuum-Tube Voltmeters. SEVILLE 
CHAPMAN, Stanford University—In d.c. vacuum-tube 
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voltmeters having a balanced double-triode circuit, the 
meter is usually connected between the plates of the two 
tubes. If, however, the load resistors are inserted in the 
cathode circuits, rather than in the plate circuits, and the 
meter is connected between the two cathodes, the plates 
being connected directly to power supply, it is shown that 
without loss of sensitivity the departure from linearity 
between input voltage and meter reading is reduced by a 
factor of from 2 to 4 in practical cases. Thus an increase 
in sensitivity by a similar factor can be obtained, or a less 
sensitive meter may be used for performance equivalent 
to that of the usual circuit. When degeneration from 
additional cathode resistors is employed in the usual 
circuit, its performance is improved, but it does not equal 
that of the suggested circuit. Practical examples are 
compared with theoretical calculations. With an ordinary 
250-volt unregulated power supply the suggested circuit 
will give a meter deflection of 1 milliampere per volt 
input with a single 6F8-G tube, sensitivity being inde- 
pendent of normal line voltage fluctuations. This high 
sensitivity arises from the low impedance of cathode 
loaded stages. 


9. The Spirotron. LAuRENCE E. Dopp, University of 
California at Los Angeles.—As a means, theoretically 
possible, of stopping very high speed particles by decelera- 
tion using potential gradients applied in successive steps 
as in the cyclotron, but at the same time avoiding corre- 
spondingly high frequencies of field reversal, which may 
be impractical to obtain, the present paper proposes to 
whittle down such high velocity a selected component at 
a time. The name “‘spirotron”’ is given to such a design. 
The type of spiral implied by this term is that of a tapering 
helix rather than a plane curve as at present. A succession 
of several of these in series, of different pitches, with their 
axes at definite angles of inclination with each other, 
should serve to stop particles of higher velocity than could 
be dealt with in the reversed cyclotron because of frequency 
difficulties. 


10. A Simple Adaptation of the Ultraviolet Photographic 
Spectrophotometer to Direct Visual Measurement. Louis 
A. Strait, W. D. KuMLER, AND FRANK M. Goyan, College 
of Pharmacy, University of California.—The conventional 
u.v. photographic spectrophotometer of the rotating sector 
or stationary variable diaphragm type is adapted to direct 
visual measurement by substituting a fluorescent screen for 
the photographic plate. Monochromatic measurement of 
optical density or transmission may be made with speed 
comparable to conventional photoelectric photometers and 
there is an advantage in flexibility with regard to type and 
size of sample which may be measured. A quick survey of 
qualitative absorption may be made conveniently by 
using a Baly tube and observing the absorption spectrum 
on the screen. Where photographic photometry of an 
unknown is to be carried out, the desired concentration 
may be readily found in the same way, saving photo- 
graphic plates as well as time. Potassium chromate solu- 
tions ranging in optical density at 2800A from 0.04 to 1.3 
have been used, and the duplicate-measured transmission 


values deviate from the theoretical by a maximum of 3 
percent with a mean deviation of 1.6 percent. More 
detailed data will be reported elsewhere. A suitable screen 
may be made by depositing a uniformly fine suspension of 
anthracene (in a base such as cellulose acetate or butvar) 
on a glass screen made by removing the emulsion from a 


photographic plate. 


11. The Measurement of Absolute Efficiencies of 
Gamma-Ray Counters. WENDELL PEACOCK, Massachusetts 
Institute of Technology.—By measuring various radioactive 
isotopes, using coincidence counting methods, one may 
determine the absolute efficiency of gamma-counters and 
the absolute rate of disintegration of the radioactive 
isotopes. By the use of the isotopes Na™, Co, Mn*, ja, 
Y*, Fe®*, Mn**, and Co**, the absolute efficiency of screen 
wall copper, platinum, and bismuth counters was deter. 
mined. With the isotopes Br®, Au’, Mn®, Na®, ang 
V*8, the efficiency curves thus determined were confirmed, 
The efficiency curves thus obtained were in substantial 
agreement with the theory proposed by V. Droste. With 
several isotopes, each emitting a single monochromatic 
gamma-ray, it is possible to compare the efficiency of 
gamma-ray detecting instruments at various energies 
rather quickly. Isotopes with conveniently long half-lives 
and emitting gamma-rays of energies above 1.5 Mev are 
not yet available, but Be’, Cu, Mn, and Zn*® serve at 
lower energies. By the use of these isotopes in connection 
with Y* and Na*, the efficiency curve can be extended 
from 0.4 Mev to 2.75 Mev. The absolute efficiency of the 
measuring device can be obtained with a source for which 
the absolute disintegration rate is known. For this purpose 
Co*® (5.3y) has been used. 


12. Direct and Fluorescence Excitation of the Z11: Level 
in Thick Targets of Thorium. Burr G. BURBANK, Stanford 
University (Introduced by Paul Kirkpatrick).—The ratio 
P of the probability of exciting the Lim shell of thorium 
by direct electron bombardment to the probability of 
exciting the same shell by photoelectric ionization in a 
thick target was determined by observation of normally 
emergent x-rays from an oil-cooled thorium-target x-ray 
tube. The intensities at the peak of the La line and in the 
continuous background on each side were measured with 
a Bragg spectrometer. The target was then covered with 
a palladium foil thick enough to stop all bombarding 
electrons, and these intensities were again measured. The 
La; line from the covered target could be produced only 
by the photoelectric ionization process, while the La; line 
from the bare target could be produced partly by each 
process. By taking into consideration the relative intensities 
and the mean depths of production of the continuous rays 
from the palladium and the thorium, P was calculated for 
a range of voltages from 40 to 100 kv, and was found to 
decrease from a value of 3.1 at 40 kv to 1.98 at 100 ky. 


13. A High Voltage D.C. Regulator Entirely A.C. Oper- 
ated. J. F. McGee, Stanford University —The regulator 
is of the fundamental degenerative type, as catalogued by 
Hunt and Hickman' for low voltage applications and 
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suggested for higher voltages. It is used in conjunction 
with a half-wave voltage-doubler power supply in the 
x-ray research laboratory. The present regulator circuit is 
distinguished from others of comparable stabilization in 
that: (a) it is entirely a.c. operated, the usual reference 
battery of several hundred volts being replaced by the 
highly regulated and separate plate supply voltage of the 
dc. amplifier; (b) it makes use of an RC combination 
which for short duration fluctuations of voltage has the 
effect of increasing the fraction of the high voltage that is 
compared with the reference voltage. 


iF. V. Hunt and R. W. Hickman, Rev. Sci. Inst. 10, 6 (1939). 


14. Planck’s Universal Constant of Action. SamueL R. 
Cook, Sacramento, California.—The paper purports to 
find a relation between Planck’s two constants, the unit 
electric charge and what Planck called his universal 
“quantum of action’’ h. Planck considered that the two 
constants were independent, and that although the two 
constants were essential in developing his radiant heat 
formula, they were entirely independent of each other. 
The paper purports to show, first, that Planck's first 
constant, which he considered as a mere formal constant, 
is in reality the real universal constant of action, being 
the unit charge of electricity e. The paper also purports to 
show, second, that the ratio of the values of the two 
constants ¢ and h is a definite constant of proportionality k. 
Since the discovery by Planck of the constant h, which 
when multiplied by the frequency of radiation gives the 
value of the energy of any radiation quantum hy, the 
exact explanation of the nature of the constant h has been 
in question. The paper shows that e, the electric charge 
on an electron, is the real ‘‘constant of action” and that h 
is a constant which acts as a factor to measure the energy 
developed in the radiation quantum developed by e. The 
ratio of e and h is determined as a constant k. 


15. A New Method of Determining the Fraction of 
Electrons in Cosmic Radiation. Ropert R. Brown, 
University of California (Introduced by Wayne E. Hazen). 
—A cloud chamber containing seven lead plates, three of 
0.2-cm thickness and four of 0.7-cm thickness, was operated 
by a three-counter telescope above the chamber. The 
counters were made of thin glass and the usual copper 
cylinder was replaced by a layer of graphite paint. Above 
the counters and the chamber was a thin metal roof. The 
equivalent stopping power of the thin metal roof, the walls 
@ the counters, and the wall of the cloud chamber was 
~10 Mev. The electrons and mesotrons that entered the 
chamber were differentiated by their shower effects, 
penetration, and ionization. The probability that an 
electron would pass through the chamber without initiating 
a shower is extremely small. On the other hand, the 
probability that a mesotron would initiate a shower 
while passing through the absorber is extremely small. The 
advantage of this method over that using a cloud chamber 
and a magnet is that the incoming particle spectrum is not 
distorted by a strong magnetic field. The advantage of 
this direct cloud-chamber method over the indirect counter 
method is that the cloud-chamber results depend on the 
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cascade theory of showers only so far as the determination 
of the lower energy limit is concerned. 


16. Rigorous Calculation of an Electrostatic Immersion 
Lens. Rupotr G. E. Hutter, Division of Electron Optics, 
Stanford University.—Previous investigations have shown 
that the electrostatic immersion lens, given by the axial 
potential distribution, ¢(z) =¢oe* ***** */*, is the simplest 
electrostatic lens for which Newton's image equations are 
satisfied throughout the field. The axial potential distribu- 
tion of this lens is found to approximate that of the 
electrostatic lens consisting of two cylinders of unequal 
radii which are connected to two potentials V;=¢oe~**” 
and V2=¢ee***?, The general solution of the paraxial ray 
differential equation is derived. Based on this solution and 
the fact that Newton's image equations are satisfied, exact 
expressions are derived for all‘important optical quantities 
such as the focal lengths, the location of the four cardinal 
points, the magnification, the object-image relation, and 
the spherical and chromatic aberration for any object 
position. The optimum positions of an object are deter- 
mined so that the chromatic and spherical aberrations are 
reduced to a minimum. All quantities are represented 
graphically as functions of the ratio of the asymptotic 
voltages V2/V;. A plot of the equipotential lines in space 
is given for a special value of the voltage ratio. 


17. Theoretical Efficiencies of Continuous X-Ray Pro- 
duction. PAuL KirKPATRICK, Stanford University—Com- 
putations of continuous x-ray intensities on the basis of 
the Sommerfeld theory, as supplemented by subsequent 
contributors, have been completed for potentials to 90 kv, 
wave-lengths from 0.181A to 1.431A, and representative 
atomic numbers throughout the range. By combining 
these results with Williams’ theory of the electron stopping 
power of matter, values of production efficiency have been 
obtained. The following table gives the ratios of emitted 
x-ray energy (all frequencies and all directions) to total 
energy loss suffered by electrons of the stated energies when 
traversing matter of the given Z. 


Zz 40 kv 90 kv 
13 0.0013 0.0018 0.0028 
47 0.0048 0.0068 0.0115 
70 0.0069 0.0097 0.017 
92 0.0087 0.0125 0.021 


The total continuous spectrum energy emitted during 
complete stoppage of a 90-kv electron is almost directly 
proportional to the atomic number of the stopping element, 
rising to 1.90 10~* erg at Z=92. 


18. The Mean Depth of Production of Continuous 
X-Rays in Thick Targets of Thorium. Burr G. BuRBANK, 
Stanford University (Introduced by Paul Kirkpatrick).— 
A thick target of thorium was bombarded by 40-100 kv 
electrons incident obliquely in a direction making a grazing 
angle of 30° with the target surface. The intensities of 
several selected wave-lengths of the continuous spectrum 
propagated in a direction at right angles to the direction 
of bombardment were observed with a Bragg spectrometer 
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for a variety of orientations of the target with respect to 
an axis of rotation coinciding with the electron beam. By 
Stoddard’s modification of Kulenkampff's method, the 
mean depths of production were deduced from the slope of 
a plot of the logarithm of the observed continuous intensity 
vs. the secant of the angle of rotation of the target. Mean 
depth values were obtained as functions of wave-length 


and bombardment voltage. The average value of the mean 
depth of production within the spectral range between the 
high frequency limit and the Lim absorption limit increases 
from 6.0X 10-5 cm at 40 kv to 22.2X10-* cm at 100 ky, 
The total range along the path of the electrons in the 
thorium, as calculated by Williams’ formula, is about 19 
times the average mean depth of production. 
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